FACILITY FORM 602 


srwrw ' ‘ r ■»' * 


i \ » 1 


ft AS# C.A 


% 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


Space Programs Summary 37-42, Vo/. Ill 

The Deep Space Network 

For the Period September 1 to October 31, 1966 


N6 7 15901 

(ACCESSION NUMBER* 

/4P 

N67 

15929 

(THRU) 

/ 



(CODE) 

(NASA Cft OR^TMX OR AD NUMBER) 

GPO PRICE 

$ 

fCATEGORY) 


CFSTI PRICE(S) $ 

Hard copy (HC) t ? sTT) 
Microfiche (MF> 


653 July 65 

JET PROPULSION LABORATORY 

CALIFORNIA INSTITUTE OF TECHNOLOGY 
PASADENA, CALIFORNIA 


November 30, 1966 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


Space Programs Summary 37-42, Vol. Ill 

The Deep Space Network 

For the Period September 1 to October 31, 1966 


JET PROPULSION LABORATORY 

CALIFORNIA INSTITUTE OF TECHNOLOGY 
PASADENA, CALIFORNIA 


November 30, 1 966 


Space Programs Summary 37-42, Vol. Ill 


Copyright © 1 967 
Jet Propulsion Laboratory 
California Institute of Technology 

Prepared Under Contract No. NAS 7-100 
National Aeronautics & Space Administration 


Preface 


The Space Programs Summary is a six-volume bimonthly publication designed 
to report on JPL space exploration programs and related supporting research and 
advanced development projects. The titles of all volumes of the Space Programs 
Summary are: 

Vol. I. The Lunar Program (Confidential) 

Vol. II. The Planetary-Interplanetary Program (Confidential) 

Vol. III. The Deep Space Network (Unclassified) 

Vol. IV. Supporting Research and Advanced Development (Unclassified) 
Vol. V. Supporting Research and Advanced Development (Confidential) 
Vol. VI. Space Exploration Programs and Space Sciences (Unclassified) 

The Space Programs Summary, Vol. VI, consists of: an unclassified digest of 
appropriate material from Vols. I, II, and III; an original presentation of the 
JPL quality assurance and reliability efforts, and the environmental- and dynamic- 
testing facility-development activities; and a reprint of the space science instru- 
mentation studies of Vols. I and II. 


Approved by: 



Jet Propulsion Laboratory 


JPL SPACE PROGRAMS SUMMARY 37-42, VOL. HI 


Contents 


I. Introduction 


II. Tracking and Navigational Accuracy Analysis 

A. DSN Inherent Accuracy Project 

T. W. Hamilton and D. W. Trask 

B. Effect of Constant Acceleration on the Information Content 
of a Single Pass of Doppler Data 

T. Nishimura 

C. Entry Capsule Orbit Determination Accuracy 

C. Scranage 

D. The JPL Matrix Operations Programming System 

P. M. Muller 

E. Theoretical Basis for the Double Precision Orbit Determination 
Program: V. Doppler and Range Observations 

T. D. Moyer 

F. Pioneer Project Support 

0. W. Curkendall, J. E. Ball, and J. F. Gallagher 

References 



,/ 

*/ 


til. Communications Research and Development 


24 


A. Frequency Generation and Control: Atomic Hydrogen 
Frequency Standard 

W. H. H/ga 

B. Improved RF Calibration Techniques : Daily System Noise 
Temperature Measurements 

C. T. Stelzried 




C. Improved RF Calibration Techniques: Weekly Noise 
Temperature Calibrations 

T. Y. Otoshi 

D. Efficient Antenna Systems: Calculated Gain of the Advanced 
Antenna System 

A. C. Ludwig 

E. Efficient Antenna Systems: X-Band Gain Measurements 

D. A, Bathker, and G. S. Levy 

F. Low Noise Receivers: Microwave Maser Development 

S. M. Petty, and R. C. C/auss 

G. X-Band Cone RF Instrumentation 

C. T. Stelzried 

H. Digital Communication and Tracking : Ranging Measurement 

R. C, fausworfhe 


3^/ 

37X 

49/ 

42X 

4 X 


iv 


JPL SPACE PROGRAMS SUMMARY 37-42, VOL. Ill 


Contents (contd) 


I. Digital Communication and Tracking: Ranging Subsystem 
Acquisition Probabilities 

R. Didday 

J. Digital Communication and Tracking: a Phase-Locked Receiver 
Analysis Program 

R. C. Tausworfhe 

K. Digital Communication and Tracking: Time Synchronization 
Experiment 

W. Martin ' . . 

L. Frequency Generation and Control: VCO Noise Measurement 
and Optimum Tracking Filters 

C. Cooper 

M. Frequency Generation and Control: S-and X-Band Central 
Frequency Synthesizer 

G. U. Barbani . 

References 


56 ^ 



</ 

(j/ 



IV. Communications Engineering Development 

A. S-Band Receiver-Exciter Ranging Subsystem (RER), 
Blocks IIA and MIC 

H. Donnelly , C. E. Johns, R. W. Tappan, and M. R. Wick 

B. 57/221 Frequency Shifter— Block MB 

C. E. Johns 

References 


77 



V. Deep Space Station Engineering and Operations 87 


A. Flight Project Support 

J. OrJb/son • 

B. Facility Construction and Equipment Installation 

J. Orbison 

C. Mariner Venus 67 Project: Ground Instrumentation for 
Occultation 

G. Levy 

D. Venus DSS Operations 

M. A. Gregg , E. 8. Jackson, and A. L. Price 

E. DSIF FR-900 Video Recorders 

W. C. Allen 

F. DSIF Station Control Data and Equipment 

R. Flanders, G. Jenkins, A. Burke, W. Frey, and H. Baugh .... 

G. Conference Nets of the Operational Voice Communications 
Subsystems 

F. E. Bond, Jr., and J. W. Capps 


3 ^ 

3 ^ 

/ 

■/ 


JPL SPACE PROGRAMS SUMMARY 37-4 2, VOL. Ill 


I. Introduction 


The Deep Space Network (DSN), established by the 
NASA Office of Tracking and Data Acquisition, is under 
the system management and technical direction of JPL. 
The DSN is responsible for two-way communications 
with unmanned spacecraft travelling from approximately 
10,000 miles from Earth to interplanetary distances. 
Tracking and data-handling equipment to support these 
missions is provided. Present facilities permit simultane- 
ous control of a newly launched spacecraft and a second 
one already in flight. In preparation for the increased 
number of U.S. activities in space, a capability is being 
developed for simultaneous control of either two newly 
launched spacecraft plus two in flight, or four spacecraft 
in flight. Advanced communications techniques are being 
implemented to make possible obtaining data from, and 
tracking spacecraft to, planets as far out in space as 
Jupiter. 

The DSN is distinct from other NASA networks such 
as the Space Tracking and Data Acquisition Network 
(STADAN), which tracks Earth-orbiting scientific and 
communication satellites, and the Manned Space Flight 
Network (MSFN), which tracks the manned spacecraft 
of the Gemini and Apollo programs. 

The DSN supports, or has supported, the following 
NASA space exploration projects: (1) Ranger , Surveyor , 
Mariner , and Voyager Projects of JPL; (2) Lunar Orbiter 
Project of the Langley Research Center; (3) Pioneer Proj- 
ect of the Ames Research Center, and (4) Apollo Project 
of the Manned Spacecraft Center (as backup to the 
Manned Space Flight Network). The main elements of 
the network are: the Deep Space Instrumentation Fa- 
cility (DSIF), with space communications and tracking 
stations located around the world; the Ground Communi- 
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cations System (GCS), which provides communications 
between all elements of the DSN; and the JPL Space 
Flight Operations Facility (SFOF), the command and 
control center. 

The DSIF tracking stations are situated such that three 
stations may be selected approximately 120 deg apart in 
longitude in order that a spacecraft in or near the ecliptic 
plane is always within the field of view of at least one 
of the selected ground antennas. The DSIF stations are: 


Deep Space 
Communication 
Complex (DSCC) 

Deep Space 
Station (DSS) 

DSS serial 
designation 

Goldstone 

Pioneer 

1 1 


Echo 

1 2 


Venus 

13 


Mars 

14 

Canberra 

Woomera 

41 


Tidbinbilla 

42 


Booroomba* 

43 


Johannesburg 

51 

Madrid 1 

Robledo 

61 


Cebreros ! 

62 


Rio Cofio' 

63 


Cape Kennedy 

71 


(Spacecraft Monitoring) 



Ascension Island 

72 


(Spacecraft Guidance 



and Command) 


1 Planned. 



-Station not yet authorized. 


"Station not yet operational. 
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JPL operates the U.S. stations, and will operate the 
Ascension Island Station. The overseas stations are nor- 
mally staffed and operated by government agencies of 
the respective countries, with the assistance of U.S. sup- 
port personnel. 

The Cape Kennedy Station supports spacecraft final 
checkout prior to lapnch, verifies compatibility between 
the DSN and the flight spacecraft, measures spacecraft 
frequencies during countdown, and provides telemetry 
reception from lift-off to local horizon. The other DSIF 
stations obtain angular position, velocity (doppler), and 
distance (range) data for the spacecraft, and provide 
command control to (up-link), and data reception from 
(down-link), the spacecraft. Large antennas, low noise 
phase-lock receiving systems, and high-power transmit- 
ters are utilized. The 85-ft diam antennas have gains 
of 53 db at 2300 MHz, with a system temperature of 
55 °K, making possible the receipt of significant data 
rates at distances as far as the planet Mars. To improve 
the data rate and distance capability, a 210-ft diam 
antenna has been built at DSS 14, and two additional an- 
tennas of this size are planned for installation at overseas 
stations. 

In their present configuration all stations, with the 
exception of Johannesburg, are full S-band stations. The 
Johannesburg receiver has the capability for L- to S-band 
conversion. 

It is the policy of the DSN to continuously conduct 
research and development of new components and sys- 
tems and to engineer them into the network to maintain 
a state-of-the-art capability. Therefore, the Goldstone sta- 
tions are also used for extensive investigation of space 
tracking and telecommunications techniques, establish- 
ment of DSIF/spacecraft compatibility, and development 
of new DSIF hardware and software. New DSIF-system 
equipment is installed and tested at the Goldstone facili- 
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ties before being accepted for system-wide integration 
into the DSIF. After acceptance for general use, it is 
classed as Goldstone Duplicate Standard (GSDS) equip- 
ment, thus standardizing the design and operation of 
identical items throughout the system. 

The GCS consists of voice, teletype, and high-speed 
data circuits provided by the NASA World-Wide Com- 
munications Network between each overseas station, the 
Cape Kennedy Station, and the SFOF. Voice, teletype, 
high-speed data, and video circuits between the SFOF 
and the Goldstone stations are provided by a DSN micro- 
wave link. The NASA Communications Network is a 
global network consisting of more than 100,000 route mi 
and 450,000 circuit mi, interconnecting 89 stations of 
which 34 are overseas in 18 foreign countries. It is 
entirely operationally oriented and comprises those cir- 
cuits, terminals, and switching equipments interconnect- 
ing tracking and data acquisition stations with, for 
example, mission control, project control, and computing 
centers. Circuits used exclusively for administrative pur- 
poses are not included. 

During the support of a spacecraft, the entire DSN 
operation is controlled by the SFOF. All spacecraft com- 
mand, data processing, and data analysis can be accom- 
plished within this facility. The SFOF, located in a 
three-story building at JPL, utilizes operations control 
consoles, status and operations displays, computers, and 
data-processing equipment for the analysis of spacecraft 
performance and space science experiments, and com- 
munications facilities to control space flight operations. 
This control is accomplished by generating trajectories 
and orbits, command and control data from tracking 
and telemetry data received from the DSIF in near-real 
time. The SFOF also reduces the telemetry, tracking, 
command, and station performance data recorded by the 
DSIF into engineering and scientific information for 
analysis and use by scientific experimenters and space- 
craft engineers. 
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II. Tracking and Navigational Accuracy Analyses 


A. DSN Inherent Accurc$&Jfect, 1.5. 

and D. W. Trask ** 

The DSN Inherent Accuracy Project was formally 
established by the DSN Executive Committee in July 
1965. The objectives of the Project are: 

(1) Determination (and verification) of the inherent 
accuracy of the DSN as a radio navigation instru- 
ment for lunar and planetary missions. 

(2) Formulation of designs and plans for refining this 
accuracy to its practical limits. 

Achievement of these goals is the joint responsibility 
of the Telecommunications Division (3*3) and the Systems 
Division (31) of JPL. To this end, regular monthly meet- 
ings are held to coordinate and initiate relevant activities. 
The Project leader and his assistant (from Divisions 31 
and 33, respectively) report to the DSN Executive Com- 
mittee, and are authorized to task Project members to 
(1) conduct analyses of proposed experiments, (2) pre- 
pare reports on current work, and (3) write descriptions 
of proposed experiments. The Project is further author- 
ized to deal directly with those flight projects using the 
DSN regarding data-gathering procedures that bear on 
inherent accuracy. 
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The various data types and tracking modes provided 
by the DSIF in support of lunar and planetary missions 
are discussed in SPS 37-39, Vol. Ill, pp. 6-8. 

Technical work directly related to the Inherent Ac- 
curacy Project was presented in SPS 37-38, Vol. Ill, and 
in subsequent issues, and is continued in the following 
sections of this volume. 

N67 15903 

B. Effect of Constant Acceleration on the 
Information Content of a Single Pass of 
Doppler Data, T. Nishimura 

1. Introduction 

In SPS 37-33, Vol. IV, pp. 8-16, J. O. Light has pointed 
out that the velocity parallax effect due to the rotation 
of the ground station around the Earth’s spin axis plays 
a dominant role in the orbit redetermination when the 
information is supplied by a single pass of doppler data 
from a distant spacecraft. Then T. W. Hamilton and 
W. G. Melbourne have clarified (SPS 37-39, Vol. Ill, pp. 
18-23) that the mean geocentric radial velocity, its right 
ascension, and the cosine of its declination are three ma- 
jor parameters that can be meaningfully measured from 
a full pass of doppler data. D. Curkendall (SPS 37-41, 

3 


Vol. Ill, pp. 47-50) has extended the analysis to the case 
when the noise in information may be regarded as white 
range noise. 

It is the objective of this article to evaluate the extent 
to which a small constant acceleration acting on the space- 
craft (solar pressure, gas leakage) will affect the esti- 
mates of the three major parameters mentioned above, 
especially the cosine of declination of the spacecraft s 
trajectory. 

A definition of terms used appears at the end of this 
article. 

2. Mathematical Model 

The assumptions made in this article are identical with 
those of Hamilton and Melbourne; namely, the space- 
craft is moving away from the Earth with an almost 
straight line motion while the observer (tracking station) 
is rotating around the z-axis at a distance r s (Fig. 1). 

The range p from the observer to the spacecraft is 
given by 

p 2 = (x - r 8 cos 0) 2 4- (y ~ r a sin S) 2 4 (z - z a ) 2 (1) 

Differentiating Eq. (1) with respect to time yields 
pp = rf — r 8 (x cos 0 4 y sin 0) — zz« 

4 r 8 $ (x sin 0 — y cos 0) (2) 


The speed components x, y, z, are described as follows 
under the constant acceleration 


X = x a + L x (t - 

to) 

(3) 

y =y 0 + Ly(t - 

to) 

(4) 

2 — Zq ( t 

to) 

(5) 

and the angular rate 0 is given by 



0 — o> 


(6) 

However, for the probe at a distance r » r 8y r » z 8r 

Eq. (2) 


can be simplified by dividing by p, expanding it in terms 
of 1/r, then dropping higher-order terms. It is also 
assumed that the positional coordinates of the probe at 
time t do not deviate significantly from the original 
straight line at time t<, so that the same rotational cosines 
are employed throughout a pass 

p = r + r 8 <o cos 8 cos 0 O sin 0 — r g to sin 0 O cos 8 cos 0 (7) 


z 



Fig. 1. Coordinates of probe P and observer O 


Let f 0 and L r , respectively, be the projections of the 
velocity vector and the constant acceleration vector on r 
vector at time t 0 . Then 

f — f 0 4 h r (£ — t 0 ) (8) 

Now, by employing the fact that t 0 = 0, and adding 
the white doppler noise n(t), Eq. (7) can be rewritten 


in the following form: 

p = a4bsin04 c cos 0 4 d 6 4 n (9) 

where 

a = r 0 — L r t 0 = f> (10) 

b — r 8 (o cos 8 cos 0 o r « w cos ^ (H) 

c = —r 8 (i> cos 8 sin B {) 4 — B {) r 8 o cos 8 (12) 

d = L r u (13) 


3. Minimum Variance Filter Design 

When we have almost continuous doppler data (e.g., 
1 sample per minute) ^contaminated by white noise, con- 
tinuous estimators a, b , c, and d can be designed which 
yield minimum variance estimates of constant parameters 
a , b, c , and d, respectively. 
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Following the standard design procedure of the minimum variance filters, the covariance matrix A(^) associated 
with the estimator is derived as 


AM = [A(O)- 1 + 7M] -1 


(14) 


JM = ~ r M8)A(0)de 

J-if; 


(15) 


where A(0) is the observation matrix defined as 


A($) = [1, sin 0, cos $, 0] 


(16) 


and \j/ is the half -width of the pass ( — ip < 0 < \p) and is thepower spectral density of white noise n. When N obser- 
vations are made over a pass at an equally spaced time interval T, the mean squared value of of the discrete data 
noise is related to the power spectral density $$ by the following: 


= *L = jV' 

rrr * T 


"2T 


(17) 


When no a priori information is available, A(0) -1 becomes zero, so that A (</*) is simply reduced to 


am = 


(18) 


assuming J(ip) is nonsingular. 

Carrying out the integration of Eq. (15), A (\fi) is found to be 

1 




1 - p 2 
r 13 

0 

— 2 sin \p 


0 
2 \p 


— 2 sin \p 


0 


(<P - % sin 2\p) (1 - p 2 2i ) 

0 


(if/ + %sin2i p) (1 — p 2 i3 ) 

q — 6 (sin \p — \p cos \p) 

2\p 


(*P + Ms sin 2ip) (1 — p 2 3 ) v (\p 4- Ms sin 2^) (1 — p 2 3 ) 

— 6 (sin \j/ — \p cos 


0 


- % sin 2\j/) (1 - plj 


0 


- %sin2f) (1 - plJ 
0 
3 


f-(i -p 2 J 


where 


= - sin f_ l y \ 

* \ i+ 2^ sin2 V 


sin ip — \p cos ip 


/ w y 

\ip — M> sin 2^ y 


(19) 


( 20 ) 


(21) 


When we compare this with the results of Hamilton and Melbourne (Eq. 3, SPS 37-39, Vol. Ill, p. 19), distinguish- 
ing the components of their matrix by putting a bar at the top, we observe 


— vfr 
cr'e' = V? 


( 22 ) 

(23) 
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CORRELATION FACTORS 


but 


of : 




(i - pir 


> for \p 2i > 0 


(24) 


The estimators a, b, c, and d are also constructed, setting their a priori values equal to zero. 

* = w^hmUl Me) d> ~ m cos H 


A 

b = 


£ ■ 


d = 


- % 


- sin 2t/<) (1 - p* 4 ) ^ S1 


sin — 


3(sin \j/ — \j/ cos i/Q f ^ * 


rx P 


J* P(0)$d9^ 


. . . . , . . . p(0) cos OdO — S * n — f p(0) dO 

(+ + ^ sm 2«p) (1 - p-’ :l ) [_ U * ] 


W (1 - P\ 


ulll M 


m - ^*-*™*) r p W si nM ; 

* ~ /2 sm 2«p rv 


(25) 

(26) 

(27) 

(28) 


The correlation factors p u and p 24 are plotted in Fig. 2, 
and square roots of the diagonal elements of the A matrix 
are depicted in Fig. 3. 

Suppose that the doppler data are processed as if there 
were no acceleration, whereas the spacecraft is actually 
subject to an acceleration L r . Then the actual estimate 
of b, which is denoted b ay would be different from %, the 
optimal estimate for the no-accelcration case. According 
to Hamilton and Melbourne, b is given by 

% = 77—: — 7- / p(0) sin OdO (29) 

~ % sin 2i// v 7 



6 



half-width of pass \p 
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Then b a is computed using p of Eq. (9) for nonzero d as 
the actual data in the above equation. 

The variance associated with this actual estimate % a , 
denoted as A£, is an ensemble average of — b) 2 , which 
is found to be 


— Jr! + 


4 d 2 (sin if/ — \p cos i j /) 2 


(i// — V 2 sin %p) 
When \j/ = 90 deg, or 2 AAft becomes 
16L; 


— o-| + AAft (30) 


ar 2 AA$ = 


4 0.64 X 10 17 X L 2 (m 2 ) (31) 


f or xp = 90 deg. We may therefore safely conclude that 
the deviation term may be neglected when L r is smaller 
than 10~ 10 g, but cannot be ignored if the acceleration is 
larger than this value. Hence, the formulation developed 
in this article becomes necessary. 

The concept of the observability which has been first 
introduced by Kalman (Ref. 1) is closely related to the 
preceding analysis. According to Kalman, the param- 
eters are observable if the matrix / in Eq. (15) which is 
called Fishers information matrix, is positive definite, or 
equivalently if all the elements of the observation matrix A 
are linearly independent of each other (Ref. 2). 


4. Discussion 

The variances <r| and <r|, respectively associated with 
the estimate of the geocentric radial velocity and the right 
ascension, are not affected by the constant acceleration 
because the even functions 1 and cos 0 over — if/ < 0 < \j/ 
are orthogonal to the odd function 6 . Therefore, these two 
quantities remain identical with those obtained by Ham- 
ilton and Melbourne. 

The variance <r,| associated with the estimate of the 
cosine of declination is severely degraded by the presence 
of constant acceleration because the odd function sin 0 is 
nonorthogonal to the odd function 0 over — \p< 0 < if/. 
These two quantities are strongly correlated over the 
intervals 0 < ^ < 90 deg, as observed from Fig. 2, where 
the correlation factor p 2 4 lies in the range — 1 < p 2 4 < 
— 0.99. This indicates that a small acceleration acting on 
the spacecraft in the radial direction may seriously dam- 
age the precision of the cosine of declination. For example, 
for if/ = 90 deg, or 1 o-ft jumps from 0.725 m without accel- 
eration, to 6.027 m with acceleration. Care must therefore 
be exercised in determining cos 8 when the spacecraft is 
subject to an acceleration. 

If the radial acceleration L r is so small that the term 
associated with it becomes the same order of magnitude 
as higher-order terms in the expansion of Eq. (7), which 
have been neglected as small, the preceding discussion 
does not apply. However, the original model in Hamilton- 
Melbourne's article, as well as their conclusion, remains 
valid. 

For example, when 

Lr = 10- 8 5 (m/sec 2 ) = KT 9 5 g 

the deviation term or 2 AAft of Eq. (30) becomes compar- 
able to <o~ 2 o-| which is (0.725 m) 2 , according to Fig. 3 


Although 0 and sin 0 are linearly independent of each 
other, and hence d and b are observable according to 
the above definitions, they are so closely correlated over 
— ^ < 0 < $ and 0 deg < \p < 90 deg, as indicated by 
their correlation p 2 4 , that they behave almost like unob- 
servable (or undistinguishable) parameters. This accounts 
for the poor performance revealed by the estimators ^ 
and d (Fig. 3), especially for small 

However, if some a priori information is available con- 
cerning these parameters, a simple analysis using Eq. 
(14) will indicate that the performance of these estimators 
can be significantly improved, and the curve for ar'of 
will approach that of orV# in Fig. 3. 


Definition of terms 

x, y , z 

Ts 

rectangular coordinates of the spacecraft 

distance of the station from the Earth's 
spin axis 

distance of the station from the Earth's 
equatorial plane 

O) 

Earth's rotation rate 

r, f 

geocentric range and range-rate of the 
spacecraft 

8 

declination of the probe 

t 

time measured from the expected time 
that the probe will cross the meridian of 
the station 

to 

actual time of the meridian crossing 

X 0 , t)o, Z 0 

speed of the spacecraft at time t 0 

Lf, Ly , L z 

components of acceleration 

II 

0 

1 

CO ( t to) 
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varying drift in one-way doppler (this study is now 
underway). 


C. Entry Capsule Orbit Determination 
Accuracy, C. Scranage 


1. Introduction 

Future planetary missions will employ an entry-lander 
capsule to collect data pertaining to the climatological, 
biological, and surface conditions existing on and near a 
planet. Advanced mission planning such as that required 
for Voyager 1973 specifies that these data be transmitted 
by the DSN tracking system via two links: (1) the capsule- 
orbiting bus link (COL), and (2) the bus-Earth tracking 
station link (BEL). The presence of two vehicles in the 
vicinity of a planet presents the interesting possibility 
of using an on-board coherent two-way doppler system 
to determine range rate between the vehicles. Not only 
can this data type be used to determine the orbit of the 
capsule, but, as in the case where one spacecraft is an 
entry capsule, it can be used to estimate the atmospheric 
parameters of the planet. This observable should be 
highly accurate because of the large position and velocity 
parallax that is possible between the two moving space- 
craft. This article presents some of the results obtained 
from a preliminary investigation on the effect of high 
frequency noise of two-way coherent doppler on the esti- 
mates of the orbit of the capsule and the atmospheric 
constants. 

The material contained in this introductory article will 
encompass (1) a description of the computer programs 
and procedures, (2) a brief discussion regarding entry and 
impact statistics of the orbital parameters, and (3) a 
presentation of the estimated uncertainties on the at- 
mospheric constants. To be discussed in future articles 
will be an estimate of the atmospheric winds (most of 
this work is completed), and an investigation of time- 


The Martian atmosphere was simulated as piecewise 
continuous with a constant lapse rate. Since this is a 
feasibility investigation, a trajectory of simple geometry 
was chosen. This trajectory constrained both the bus and 
capsule to orbits lying in the equatorial plane of the 
planet; therefore, it was not possible to solve for latitude 
and flight path azimuth of the capsule with respect to 
Mars. These restrictions did not limit the generality of 
the method or the results. 

Most significant of the results is that this new observ- 
able determined the stratosphere scale height, surface 
temperature, and stratosphere base temperature to a 
greater accuracy than was attained by the Mariner IV 
Mars occupation experiment. 

2. Mathematical Simulation 

a. Trajectory program. Two-body equations of motion 
employing three degrees of freedom were used to simu- 
late the orbits of the bus and capsule. Both vehicles were 
represented as point masses in the program. Atmospheric 
drag was assumed to be the only nongravitational force 
acting on the capsule; the orbit of the bus was beyond 
the specified limits of the atmosphere. Normal trajectory 
conditions are given in Table 1. The bus-capsule geom- 
etry is illustrated in Fig. 4. 

The atmosphere was assumed to be composed entirely 
of CCL. Since a temperature maximum does not occur 
for Mars (due to lack of ozone), the adiabatic lapse rate 
for temperature was applied immediately above the sur- 
face. The stratosphere is defined as the region extending 


Table 1. Nominal trajectory conditions (parameters referenced to coordinate system centered on Mars) 


Injection To 

Entry 

Capsule 

Bus 

Capsule 

T t/ inertial longitude, 235 deg 

X/, inertial latitude, 0 deg 

yvt, inertial flight-path angle, —21.381 deg 

/3 » 7 , flight-path azimuth, 90 deg 

Ri, inertial radial distance from geocenter, 7118.1966 km 

V/, inertial velocity, 2.1224864 km /sec 

GMM, gravitational constant of Mars, 42977.8 km'V*ec 2 

a, 6478 km 
e, 0.2933004 
■,0 

Q, 270 deg 
w, 90 deg 

T p , -11681.825 sec 

Ri, altitude, 243.8 km 
time of entry 3061 sec 
7r/, inertial flight-path angle, — 15.65 deg 
V/, inertial velocity, 4.014 km/sec 
77 , inertial longitude, —33.61 deg 
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EQUATOR PLANE 


The scale height of the stratosphere H ' is given as 


INJECTION 



»J = 


T b R 

m v g b 


( 3 ) 


where m 0 is the molecular weight of the atmosphere at 
the surface and g b is the acceleration of gravity at h b . 

If m 0 gh is assumed constant for small changes in h. 


8H; =K8T b 


(4) 


where 


Fig. 4. Capsule-bus geometry 


from tropopause to mesopause. Values of atmospheric 
density and temperature at the surface and at the base 
altitude of the stratosphere (namely, p 0 , T,„ p b , T b ) and 
the temperature gradient T = + dT/d(h — h b ) are given. 
The temperature and density as a function of altitude 
are computed as follows: 


T — T 0 + r (h — h* J, h„_, <h < h* 


(2-AD 


(T nXf+Ar.) / r \ 

' *(t) (t) 

X exp j^— A (f 0 — r 6 r)r 6 ^l 


where 



( 1 ) 


( 2 ) 


K = — — 

m og b 


The ratio (standard deviation of H' b )/H' b , called coeffi- 
cient variation, becomes 



(5) 


b. Statistical program . An IBM 1620 statistical program 
(MOPOS) employing Bayes’ estimation theory was used 
to compute the uncertainties on the parameters at injec- 
tion. The effects of data noise and a priori uncertainties 
were included in the statistical model. The covariance on 
the estimates of the parameter set X are contained in the 
matrix 


A = g°*o 
R'(T b - rv 6 ) 2 

M = — — — 

(yR'TY* 

R = universal gas constant 
FP = gas constant 
R 0 = equatorial radius 
h ~ geometric altitude 

h* = base altitude for atmospheric calculations 
h b = base altitude of stratosphere 
r — R 0 + h 
r b = R 0 + h* ^ 

V A — velocity with respect to atmosphere 


T = [A T WA + T- 1 ]- 1 (6) 

where A is a matrix of partials of range rate with respect 
to X, and F x is the a priori covariance matrix. In the W 
matrix the data is weighted as the inverse of the square 
of standard deviations on the observables. In this pro- 
gram it is assumed that the observations are not corre- 
lated with each other. 

3. Procedures 

a . A priori uncertainties. The sigmas on the estimated 
parameters are affected by inherent inaccuracies of the 
BEL which limit the uncertainties in the a priori injec- 
tion conditions. The present S-band BEL configuration 
estimates position at the distance of Mars to within 
10 km (3<r). 
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For this study, the position of the bus is assumed to 
be known perfectly, and the 10 km uncertainty is assumed 
to be simulated as an error in the position of the capsule 
only. A retro-thrust of 500 m/sec and l-o- uncertainty in 
the pointing angle of % deg was used. The uncertainty 
of 1 kmVsec 2 in the gravitational constant of Mars was 
determined from Mariner IV data. From these data the 
l-o- a priori uncertainties of Table 2 were computed. 
A priori uncertainties assumed for the atmospheric con- 
stants are tabulated in Table 3. 

b. Data weighting. The high frequency noise observed 
on p in the BEL is generated by the hardware associated 
with the coherent two-way doppler. The level of this 
noise is directly related to oscillator instability (which 
varies with the type of oscillator) and oscillator tempera- 
ture. Tracking data noise of 10~' 5 km/sec (60-sec count 


Table 2. A priori uncertainties, lo 


Capsule parameter 

A priori values 

Longitude 

0.05 deg 

Latitude 

0.05 deg 

Flight-path angle 

0.1 deg 

Azimuth 

0.1 deg 

Radial distance from planet 

3.0 km 

Velocity 

0.002 km/sec 

GM Mars 

1 .0 km 3 /sec 2 


time) is considered a reasonable level for an Earth-based 
S-band system. 

For this study it is assumed that the doppler is counted 
in the COL. Since it is highly improbable that capsule 
design will permit a rigid on-board environmental control 
to maintain oscillator stability, it is equally improbable 
that the doppler counted in the capsule would be of the 
quality associated with a high frequency noise level of 
10 _r> km/sec. A knowledgeable guess of <rj = 10 3 km/sec 
is more likely. 

c. Statistics . Finite difference partials to be used in the 
A matrix were computed by the trajectory program. The 
partials were generated to simulate observations taken 
at 1-min intervals before entry and at 2-sec intervals 
through the atmosphere. The injection statistics com- 
puted by the program were mapped to entry and impact 
altitudes. Elements of the mapping matrices were also 
computed by finite differences. 

4. Discussion of Results 

a. Orbital parameters 

Statistics at entry . As mentioned previously, the geom- 
etry of the trajectory permits only four of the six orbital 
parameters to be estimated; namely, inertial longitude 
r/, inertial flight path angle y Vh radial distance from the 
geocenter R h and inertial velocity V/. Uncertainties in 
position were mapped into uncertainties in time of entry 
and time of impact. Impact statistics on all four parameters 
are given. In this article, the discussion of entry statistics 


Table 3. Coefficient of variation on atmospheric parameters 




Coefficient of variation, % 

Values obtained from Mariner 

Parameter 

Nominal 


Data noise, km /sec 

occultation experiment 

value 

A priori 

io * 2 

10 “ 

Nominal 

value 

> 

o 

7 

1,37 

0.7 

0.7 

0.2 



km 2 

R\ 5 

««e*- • K 

0.1 889 X 10“ 

15 

1 1 

2 



hi,, km 

18.6 

5 

5 

0.06 



To, °K 

200 

12 

1 1 

2 

180 

1 1 

Tt, °K 

100 

25 

13 

2 



kg 

P °' km s 

0.132 X 10 s 

70 

27 

0.3 

0.143 X 10 * 

7 

H\, km 

5 


13 

2 

8-10 

20 

I', k/km 

-5.37 
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is limited to the uncertainties in V/ and y Vh the two pa- 
rameters which are of interest to the aerodynamicist 
engaged in the design of the capsule because they are 
necessary in computing the degree of atmospheric 
“braking.” 

Figs. 5 and 6 include plots of the standard deviations 
of Vf and y Vf mapped to entry. An increase in time from 
injection corresponds to an increase in the number of 
observations at the rate of one observation per minute, 
which is used in the solution of the estimates. These plots 
indicate that uncertainties determined from tracking data 
(assuming no a priori ) with high frequency noise in the 



range 10 _J to 10 -3 km/sec are not significantly improved 
when the combination of a priori and tracking data is 
used to determine the uncertainties. It is interesting to 
note that a priori plus tracking data with o-£ = 10 -2 
km/sec produce statistics equal in magnitude to those 
obtained had only the a priori covariance matrix been 
mapped to entry. 

Impact statistics . Table 4 is presented to illustrate the 
dependence of the uncertainties at impact on the quality 
of the observations, the quantity of the observations, and 
the number of parameters estimated. Noise on the data 
of 10" 5 and 10~ 2 km/sec is given (statistics for 10" 4 and 



Table 4. Impact statistics 


Tracking 

data 

P 

km /sec 

77 

deg 

ff , 

7 vi 

deg 

^Timp' SeC 

(Tvi, km /sec 

5 

12 

5 

12 

5 

12 

5 

12 

Injection 

icr 2 

4.6 X 1 O' 2 

1.7 

2.8 X lO' 2 

1.3 X 10 1 

3.7 X 10‘ 2 

8.5 X 10“ 

4.5 X 10‘ 4 

1.8 X 10” 

to impact 

10" 5 

2.0 X 10' 4 

4.4 X 1 0 -1 

3.5 X 10' 4 

o 

X 

CO 

CO 

1.1 X 10' 2 

1.8 X 10” 

6.0 X 10 _s 

8.0 X 10“' 

Injection 

TO" 2 

3.5 X 10‘ 4 




6.0 X 10’ 1 


8.0 X 10‘ 3 


to entry 

10" 8 

3.3 X 10' s 



mBM 

1.3 X 10' 1 


1.0 X 10' 4 
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10~ 3 were obtained,but their inclusion here is not neces- 
sary for a general discussion). A priori uncertainties are 
included in the estimates. The numbers 5 and 12 indi- 
cate the number of parameters estimated: 5 = 4 orbital 
+ gravitational constant of Mars (GMM); 12 = 4 orbital 
+ GMM + angular velocity of planet « + 6 atmospheric 
constants. 

The increase in standard deviations when the 12 pa- 
rameters were estimated is due to (1) the correlation of 
the four orbital parameters with the eight constants, and 
(2) the inclusion of the a priori uncertainties in X/ and 
fi V! in the calculations. An increase in the noise on the 
tracking data from 10'" to 10~ 2 km/sec results in an in- 
crease of 11 min in the uncertainty of impact time and 
corresponding increases in the uncertainties in 17, y Vr 
and V/ of 1.3 deg, 0.13 deg and 0.017 km/sec, respectively. 

b. Atmospheric constants. The coefficients of variation 
(CV) on the atmospheric parameters determined from a 
solution vector that included the 12 parameters mentioned 
above are given in Table 3. The values listed under the 
heading “Occultation Experiment” were determined from 
the values of T 0 , p 0 , and H[ obtained from the Mariner IV 
occultation experiment (Ref. 3). In this study, H' b was 
not determined directly. The CV of H ' was evaluated by 
equation using the value of the CV of T b determined from 
the study. Judging from the accuracies attained from the 
Mariner IV mission, the a priori uncertainties appear to 
be conservative. As a result, the statistics on the atmos- 
pheric parameters can be considered realistic. The most 
significant fact obtained from this table is the degree of 
accuracy to which atmospheric constants can be esti- 
mated, even from a doppler system limited by high fre- 
quency noise to an inherent accuracy of 10 2 km/sec. 
The uncertainties on the estimates of some atmospheric 
constants obtained from this study are an improvement 
over those obtained from the occultation experiment. 


5. Conclusions 

This investigation indicates that an on-board tracking 
system of two-way coherent doppler between an entry 
capsule and an orbiting bus is a highly reliable observ- 
able for determining the orbit of an entry capsule and 
estimating atmospheric parameters. The uncertainties on 
the estimates of the parameters do not appear to be as 
dependent on the inherent accuracy of the system as 
with an Earth-based system. High fidelity estimates can 
be obtained even when such a system is limited to in- 
herent accuracy of 10 - km/sec due to high frequency 
noise. When this observable (with noise of 10' 2 km/sec) 


was used in a simulated orbit determination accuracy 
study, estimates of the surface temperature, stratosphere 
base temperature, and stratosphere scale height were 
obtained which were comparable to those attained by 
the present methods of occultation. 


The results presented here were computed using both 
statistical and mathematical models which were assumed 
to describe perfectly the atmosphere of Mars. Therefore, 
first-order changes in the model of the atmosphere would 
produce significant changes in the standard deviations. 


D. The JPL Matrix 
System, P. M. Muiler 


s^5g ra irm!?ng 


The interactions between various computer programs 
used in processing and analyzing DSN tracking data 
have been described in SPS 37-39 , Vol. Ill, pp. 34-36. 
To illustrate how these programs fit together, systems 
for processing non-real-time data were reviewed. In addi- 
tion, programs designed to establish the orbit determina- 
tion accuracies obtainable with the DSIF tracking data 
were discussed. In this paper, one of the programs for 
processing and analyzing tracking data is considered. 


1. Summary 

The purpose of this article is to briefly describe an 
automatic FORTRAN IV programming system which 
accomplishes all matrix and vector operations plus a 
sophisticated input-output and housekeeping scheme. 
This system is programmed in FORTRAN IV in such a 
way that it should be readily convertible to any modem 
computing system with a FORTRAN IV compiler and 
at least 16,000 words of core. Minor changes are neces- 
sary for such conversion between machines, or between 
single and double precision (the latter is nominal). These 
changes are all made to a small number of subroutines, 
and thereafter the system will operate correctly and 
automatically under the new requirements. Dimension- 
ing of arrays can range from 2 to 100, and is handled in 
one subroutine and the main (users) program only. The 
system allows efficient use of programming, computer, 
and documentation time, but is comparatively inefficient 
with respect to the use of core storage. With the ex- 
penditure of about 3000 words in core, any subset of the 
system may be used, and this eases somewhat the core 
requirements for small jobs. 

2. System Description 

At this time, the JPL matrix operations programming 
(MOP) system consists of some 85 subroutines plus over 
two dozen additional named entry points. This is one of 
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two basic methods of programming a system that ac- 
complishes the vector and matrix algebra operations and 
related manipulations. The other method, VECTRAN, 
to be found in the SHARE catalogue 1 , uses an addi- 
tion to the FORTRAN compiler which accepts matrix 
algebra expressions in addition to FORTRAN. The MOP 
method requires the user to code the FORTRAN pro- 
gram in which the MOP subroutines are called by the 
user. The VECTRAN type of compiler system presum- 
ably allows the user to code a FORTRAN type of pro- 
gram which has expanded capability. In practice, this 
compiler method has not proved to be completely versa- 
tile or free from errors. On the other hand, the reliability 
of the subsystem routine has been established, since it 
can and has been exhaustively checked out by an auto- 
matic debugging package. For this reason, plus advan- 
tages in facility of updating and modifying the system as 
well as greater ease of learning and use, the subroutine 
system was chosen in preference to a compiler system. 

3. Advantages and Costs 

If the MOP were merely a series of library subroutines, 
it would still have the advantages of coherent naming, 
calling sequences, and philosophy. However, it surpasses 
that role, as it also includes a large amount of automatic 
control. The error messages, page registration, page head- 
ings, output labels, parameter (row-column) names, and 
many of the subroutine conditions are all automatically 
controlled regardless of which subroutines or which sub- 
sets of the system are used. The calling sequences in- 
volve, for the most part, only the operator names, and 
are free of the control indices always found in library 
routines. The naming of the subroutines is highly mne- 
monic and lends itself to rapid mastery. The cost for all 
of this is about 3000 words of core, which must be used 
and initialized before any or all of the system can be 
effectively run. The programmer need have only a fair 
understanding of FORTRAN IV to use the basic sys- 
tem. The more advanced FORTRAN programmer can 
quickly and conveniently use the system to accomplish 
virtually any and all matrix algebra manipulations, in- 
cluding input and output, and can do so most efficiently. 

4. The System Operators 

The possible operators for the system are listed as 
follows : 

MATRIX: A square array (two-dimensional) of real 

numbers. 


’IBM Systems Reference Library “SHARE” Program Listing and 
Dictionary of Titles. 


TRIMAT: A square array matrix with the property of 
symmetry, stored in core as the triangular portion plus 
the diagonal (saves core). 

DIAG: Diagonal matrix, stored as a vector, but handled 
in subroutines as though it were a diagonal matrix. 

VECTOR: A vector of real numbers, either row or 
column vector, according to the subroutine called. 

CONSTANT : A real number to be used in some matrix- 
algebra operation. 

There is some provision for optimizing certain operations 
for nonsquare matrices stored in a square-matrix array 
(by augmenting with zeroes). Future development of the 
system will include the general rectangular matrix array 
(RECMAT). There are integer-vector arrays used for pro- 
gram control and for alphameric names and variables. 

5. System Operations 

The system includes several single-operator functions. 
They are inversion, correlation, normalization, duplica- 
tion (of an operator), conversion (of an operator from 
one form to another; e.g., TRIMAT to MATRIX), input, 
output print, output punch, summation (of the elements), 
test for nonzero elements, zero to all elements, perform 
any FORTRAN library function upon all elements of an 
operator, and a large number of reordering or partition- 
ing functions. The multiple-operator functions are addi- 
tion, subtraction, and multiplication (including nonsquare 
matrices and DIAGs used in pre-post multiply situations). 
The subroutines SETUP and IOU (Input-Output-Utility) 
are required to initialize and preserve the system. Future 
development will include subroutines for RECMAT oper- 
ations, eigenvalues and vectors, and name-list-format 
output of arrays. 

6. Automatic Functions and Controls 

The program has a common location for scratch vectors 
and scalars used by the subroutines themselves, which 
can occasionally be used in a program, but are normally 
left buried and untouched. There is a common location 
for a set of real-number indices (DPIX), fixed-point in- 
dices (INDEX), 10 word by 6 character page headings 
and labels (PHEADS), and 1 word, 6 character param- 
eter or row and column names (PLABS). All system 
constants are stored in DPIX and referred to by sub- 
routines as needed. The system control indices are in 
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INDEX. Both system-nominal and user-input page head- 
ings are addressable in PHEADS by the output routines, 
and are controlled by the user in blocks of 10 words 
from anywhere in the variable. The row and column 
names (which may be addressed independently) are taken 
from PLABS, beginning with any address chosen. These 
variables, once set by the user, conveniently allow label- 
ing of all output, automatic registration to any form 
length, page numbering, diagnostic outputs for debug- 
ging or other uses, and alterations to the nominal mode 
of operation of many subroutines. The presence of these 
controls results in very simple calling sequences. For 
example, matrix inversion (Call IM(A)), matrix multipli- 
cation (Call MXM (A, B, C)) wherein C may replace 
either operator, and many others, all allowed without 
recourse to the EPSILON increments, convergence tests, 
variable dimensions, or others indices always present in 
library subroutines. All these variables, when needed, are 
found in the INDEX variables. There they can be in- 
itialized once and forgotten, or changed at will in the 
users main program, depending upon the required pro- 
gram sophistication. All nominal system settings are made 
once in the subroutines SETUP and IOU, and thereafter 
will automatically serve the particular use in question. 

7. Additional System Design Features 

The MATRIX arrays are all dimensioned as vectors in 
the subroutines. In the main program, either single or 
double subscripting may be used for convenience. There- 
fore, the subroutines must compute their own addresses 
because of the presence of TRIMAT storage elements. 
MATRIX storage is by columns of length KD. with any 
matrix of that size or smaller stored by columns, each 
beginning at the top of the master columns and having 
length KORD where KORD ^ KD . TRIMAT storage is 
also by columns, with each successive master column 
beginning on the diagonal element instead of the first 
row. The matrix is dimensional with (KD + 1) columns 
of length KD each to permit overlay of two TRIMATs on 
one matrix. Thus, two TRIMATs of length (KD + l)*KD/2 
fit exactly (by EQUIVALENCE) on top of one MATRIX. 
This type of core allocation allows use of given core for 
either array type or for both sequentially. Most subrou- 
tines allow simultaneous use of these array types, even 
though overlaid at the time the answer is delivered. It 
is important to note that every operation was programmed 
to permit the maximum mathematically allowable over- 
lay and identity of operators. Therefore, it is possible to 
do very large matrix allocations and still use the system 
efficiently. For example, in a typical case, the IBM 7094 

14 


computer, with 32,000 words of core, may be allocated 
as follows: 10,000 words for the monitor system software, 
12,000 words for two MATRIX arrays overlaid with 
four TRIMAT arrays plus several vectors, and the 
remaining 10,000 words of core for MOP subroutines. 
It is possible to execute any of the MOP subroutines 
under this minimum configuration, and it is possible to 
load about one-third of the subroutines without overlay. 
A loading such as this permits the writing of a very large 
matrix program, with at least two MATRIX and four 
TRIMAT operators available for manipulation. This lay- 
out has been used for several NORMAL programs de- 
scribed here. 

If care is exercised, the system can be run with mat- 
rices of several master dimensions all existing at the same 
time, even if overlaid. The subroutines, with the excep- 
tion of duplication, are limited to having arrays of the 
same master dimension (KD) in any one calling sequence. 
The duplication routines are used, for example, to remove 
a 6 X 6 submatrix of a 50 X 50 and thereafter operate 
upon it as a 6 X 6 to save core. The two INDEX vari- 
ables, KD and KORD , are used by the main program to 
keep the logic operating correctly. 

It has been questioned why the system did not use 
the variable dimensions which FORTRAN IV permits, 
thereby eliminating the necessity of a master and sub- 
dimension. The answer is twofold: First, TRIMAT core 
allocation cannot be handled by convenient variable di- 
mensions; second, the requirement that all such variable 
dimensions must appear in the calling sequences was not 
considered advantageous. In addition, single-subscripted 
subroutines run from slightly faster to a factor of 4 faster 
than double-subscripted subroutines, depending upon 
the computer system used. Therefore, the MOP system 
was compiled in such a fashion that the user sets KD and 
KORD to the proper values, where they remain until a 
further change is required. 

8. Examples of MOP System Programs 

The author has written ten distinct programs using 
the system, ranging from a matrix partitioning plus 
punchout, to the entire NORMAL and STATISTICS 
link, to the selenodesy programming system which will 
process the astrodynamical results from the Lunar Orb iter 
spacecraft data. The results have encouraged repeated 
use of the system for JPL Systems Analysis Section pro- 
grams, and MOP is soon to be used as the basis for a 
very large and general NORMAL study program 
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(DPMMP). The system was completely revised into its 
present form, as documented here, in preparation for 
this effort. 

As an example of a simple but representative 
NORMAL program, a description of the program 
MOPATA follows. It uses card or tape (standardized 
general format) input of vectors of partial derivatives 
relating some observables and some parameters at a 
number of points. Each such input is called an A matrix, 
and has dimensions O by P , where O is the number of 
observables, and P is the number of parameters. W is a 
weighting matrix of dimension O by O, which is the set 
of weights assigned to the observables. The program ac- 
cumulates the product A T WA over the observations, the 
so-called normal, or J matrix. Its inverse is the covariance 
matrix on the parameters estimated. Such a result may 
be “mapped” by application of a matrix containing the 
necessary coordinate transformation or translation par- 
tial derivatives U, and computed by UITJ T where T is 
the covariance matrix in question. The answer covariance 
matrix has undergone the desired transformation. 

In addition to these computations, the user may wish 
to begin the J matrix accumulation at some point other 
than zero, and he may wish to study the effects of these 
various possible a priori conditions. The program 
MOPATA solves this general problem where any num- 
ber of a priori conditions up to 50 may be given, any set 
of mapping matrix conditions given, and where any sub- 
set or subsets of the parameters to be estimated (to be 
actually inverted) are specified. We call this a para- 
metric analysis program, since one run of the program 
results in answers involving all possible combinations of 
the user inputs, which are conveniently output for com- 
parison. This program (to a dimension of 50 X 50) was 
designed, programmed, checked out, and delivered by a 
single person in one week. It included useful page head- 
ings, labels, name list input of matrices, and many other 
features that are an integral part of the MOP system. 

there is good evidence that this system reduces by a 
factor of 4 to 10 the time required to do the same job by 
other standard means, even assuming that library sub- 
routines are available. The automatic features of the 
MOP system are convenient and useful, and allow very 
sophisticated output (for example) with minimal pro- 
gramming effort. In addition, the resulting programs are 
easily and reliably modified to suit short-term user needs. 
The MOPATA program, as well as several others, will 
be documented in future Space Programs Summaries. 


N67 15906 

E. Theoretical Basis for the Double Precision 
Orbit Determination Program: V. Doppler 
and Range Observables, T. D. Moyer 

1. Introduction 

This article is the fifth of a series presenting the theo- 
retical basis for the double precision orbit determina- 
tion program (DPODP), and gives the formulation for 
computing two-way doppler and round-trip range in 
light seconds. These are the primary observables taken 
by the tracking stations of the DSN. 

An electromagnetic signal is transmitted continuously 
from a tracking station on Earth, received and retrans- 
mitted by the probe, and continuously received by the 
same tracking station on Earth. The frequency of the 
received signal differs from that of the transmitted signal 
because of the doppler shift. Two-way doppler is the 
average value of the frequency shift (transmitted minus 
received) over a period of time T c (called the count time). 
An observable is generally recorded every T c seconds. 
The round-trip time of the signal in seconds of station 
time is the two-leg ranging observable. 

The formulation for computing two-way doppler and 
two-leg range is given in Sections 2 and 3, respectively. 
The formulation requires the solution to the light time 
problem (SPS 37-41, Vol. Ill, pp. 31-38); namely, the 
heliocentric position, velocity, acceleration, and jerk vec- 
tors of each participant evaluated at its epoch of partici- 
pation. The rectangular components of each vector are 
referred to the mean Earth equator and equinox of 
1950.0. The participants are the tracking station at the 
reception time t A (midpoint of the count interval T c for 
doppler), the probe at its reception/transmission time £>, 
and the tracking station at its transmission time t lm 

The effect of atmospheric refraction on two-way dop- 
pler and two-leg range is accounted for by an additive 
correction, which is not considered here. 

2. Computation of Two-Way Doppler 

a. General expression. An expression is given here for 
two-way doppler as a function of the doppler frequency 
shift and its second-time derivative, both evaluated at 
the midpoint of the count interval T c . Expressions for 
computing these quantities will be derived in the follow- 
ing two sections. 

The output from the electronic equipment at the track- 
ing station is a signal whose frequency in cycles per 
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second of station time is f: 



( 1 ) 


where 

f q (t 3 ), f q (t 1 ) = transmitter reference oscillator frequency 
at reception time t a and transmission 
time t u respectively (generally identical 
unless f q has been changed between 
and f 3 ) 


The observed two-way doppler is the average of f—fn a , 
over the count time T r - 



3 T c 


(f fbias)dt 


't 3 -*Tc 


( 5 ) 


where t 3 ~ epoch at midpoint of count interval, station 
time. 


f R /f T = ratio of received to transmitted frequency 
for unity frequency multiplication at 
spacecraft 

c.JKil)’ l - w 

i 06 (Si) 1 s ‘ band 

HP , L-band 

10 G , S-band 



The transmitted frequency is 30 or 96 times the reference 
oscillator frequency for L- or S-band operation. The 
spacecraft transponder multiplies the frequency of the 
received signal by 96/89 or 240/221 for L- or S-band 
operation before retransmitting. The expression for f may 
be written as 


The expression for f—fbuis may be expanded in a Taylor 
series with coefficients evaluated at f 3 . Eq. (5) may then 
be integrated term by term. The odd derivatives vanish, 
and the fourth and higher derivatives are ignored. The 
resulting expression for f is 

7= C4M + -24 ( x ~ ^r) ] 

(6) 

where [1 — (f*//r)] and its second derivative with respect 
to station time are evaluated at the midpoint of the count 
interval. The first term which has been truncated in 
Eq. (6) is (1/1920) (TJ ) [l-(/«//r)] ir . In order to mini- 
mize to truncation error in Eq. (6), the count time T c 
must be limited to approximately 1 sec for a near Earth 
spacecraft, and to approximately 1000 sec when the space- 
craft is in heliocentric cruise. 


f — ^1 + fbins ( 2 ) 

where 


b. Derivation of expression for [1 — (fi</fr)]- The ratio 
of received to transmitted frequency, where the fre- 
quencies are measured relative to the atomic frequency 
standard at the tracking station, is 


hias - CAUt 3 ) - /„(*,)] + c 4 (3) 

Note that the first term of f bi(l s normally is zero. 

The signal with frequency f is fed into an electronic 
counter whose register is incremented by 1 each time the 
amplitude of the signal changes from minus to plus. N 
cycles are counted during the count time T ( . The two- 
way doppler observable which the data editing program 
passes on to the orbit determination program (ODP) is: 

7 = jr - h,», (4) 


__ dn ' drr _ dr t ( 7 ) 

f T dr 3 dn dr 3 

where 

dn = infinitesimal number of cycles transmitted and 
received 

dr, = infinitesimal period (of proper time r) of trans- 
mission of dn cycles from tracking station at 
time t ! 

dr :i = infinitesimal period (of proper time r) of recep- 
tion of dn cycles at tracking station at time t :i 
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Note that proper time is time obtained from an atomic 
clock fixed in the observers frame of reference. Eq. (7) 
may be written as 

f R (dn/dtj / dh \ / dt 2 \ 
f T (dT S /dta) \dt 2 )\dt s ) ' ' 


where 

dt u dt 2 , dt 3 = differentials of transmission time t u probe 
time t 2 , and reception time t 3 (ephemeris 
time). 

The ratios dt^dtz and dt 2 /dt 3 will be obtained by differ- 
entiating the light time equations for the up and down 
legs of the light path. An interval of proper time is given 
by 



where 

c = speed of light 

ds = invariant interval between two near events in 
four dimensional space-time 

The interval ds of transmission (or reception) of dn cycles 
is related to the transmission time dt (ephemeris time) 
and to the differentials of the heliocentric space coordi- 
nates of the station during dt . Using de Sitter’s n-body 
line element (Ref. 4) (which reduces to the Schwarzschild 
isotropic line element when the mass of the Sun only is 
considered), and retaining terms of order (l/c)°, 


ds 2 



dt - — dx 2 — dy 2 — dz 2 


( 10 ) 


where 

x,y,z = rectangular components of station position in 
heliocentric nonrotating frame of reference 

t = uniform coordinate time in heliocentric frame 
of reference (ephemeris time) 

<j> = Newtonian potential at station given approxi- 
mately by 


Substituting Eq. (9) into Eq. (10) and using 
where 

s = heliocentric velocity of tracking station 

gives the ratio of proper time to ephemeris time for trans- 
mission or reception: 



From SPS 37-41 , Vol. Ill, pp. 31-38, the light time equa- 
tion for the up leg of the light path is given by 



where 

I* = gravitational constant of Sun, km 3 /sec 2 

^12 — 1 1 r 12 1 1 (15) 

*12 = r 2 — r, , r — » r, r, r (16) 


r,* 5 rj, r/, r! = heliocentric position, velocity, accelera- 
tion, and jerk of participant i at its 
epoch of participation t h with rectangu- 
lar components referred to the mean 
Earth equator and equinox of 1950.0. The 
index i — 1 (tracking station at transmis- 
sion time fj), 2 (probe at f 2 ), or 3 (track- 
ing station at reception time t 3 ). 

cos = o < < = 180° (17) 

r i ?12 

cos , rii . 0 < B* < = 180° (18) 

Formula (14) is indeterminate for light traveling radially 
to or from the Sun. For this case, 


*=+ 1 ? (id 

3 U 

M — gravitational constant of body /, kmVsec 2 
r j — distance from station to body /, km 


U 


tl = rn + ^±i n j± 
c c' r. 


(19) 


where rj and r s are the larger and smaller values of i\ 
and r 2 . Eqs. (14) through (19) apply for the down leg of 
the light path when the subscript 1 is replaced by 2 and 


JPL SPACE PROGRAMS SUMMARY 37-42, VOL. Ill 


17 


2 is replaced by 3. Differentiating Eq. (14) with respect 
to # 2 and collecting dtjdt 2 terms gives 
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Differentiating Eq. (17) with respect to # 2 and t, and 
adding gives 


JPL SPACE PROGRAMS SUMMARY 37-42, VOL. Ill 


Expanding Eq. (33) in powers of 1/c and retaining terms to order (1/c) 3 gives 


0- 


"^) + ' r2 *> + "^7 j^Pi2^2 + P^r 23 — n,r 23 + - <f> 3 + -i- (i* — 

+ -^r | P 12 ^12 + Pl 3 h 3 - KA 3 (Pl2 + P23) - (fi2 + r 23 ) ~ <*>3 + J (ij - s|) J 


— 2 /i (e 12 + £23) 


(34) 


where e 12 and e 23 are given by Eqs. (30) or 

(31) and 

of differentiation, the expression for [1 — (f R /f T )] is taken 
to be 

• r i2 * 

r 12 = • r 12 

^12 

(35) 

(l ~ -^-) ~ (f 12 + h 3 ) + 4 r (^12^12 + p 33 r 33 — >'l2»'=.l) 

• r 23 

r 23 = — • r 23 

^23 

(36) 

(42) 



The second derivatives of the omitted 1/c 2 potential <£ and 

* 1*12 
Pi* = ~r ,r ' 

(37) 

velocity i 2 terms and the 1/c 3 terms contribute less than 

•12 


10' 5 m/sec to observed average range rate. Differentiating 

- J*o 3 


Eq. (42) twice with respect to t 3 (ET) gives 

^ = 7 T‘ r2 

(38) 


* 2 = r; • r, 

(39) 

(‘ i) + c r ”>] 

The Newtonian potential at the tracking 
puted from 

station is com- 

+ c dtl [/ + c (P22 fl2) ] 

J — P* | PE 

1 | r* 1 1 station radius ’ 

i = 1,3 



(40) 



where 

p f *> pe = gravitational constants for Sun and Earth, km 3 /sec 2 

c. Derivation of expression for [ 1 — ( f R /f T ) ] “ . The sec- 
ond derivative of the doppler frequency shift with respect 
to station time, evaluated at the midpoint of the count 
interval, is approximated by the second derivative with 
respect to ephemeris time. 


+ 


2^ ( df 12 dp 12 
dt 3 du 




dr 23 dp 23 __ dr 12 c/r 23 \ 

3 dt 3 dt 3 J 


+ dt 3 dt 


(43) 




£ P 


dt; (ST) ( J u) 

* -0-fc) 


dt\ (ET) 


The quantities r 12 , f 23 , p a2 , and p 2 3 are functions of the 
state vectors at times f 3 , t 2 , and t x . In order to differentiate 
these quantities with respect to J 3 (ET), the following 
derivatives (obtained by differentiating the light time 
equations) are used: 


(41) 


dt 2 ( ET) 
* 3 (ET) 


1 23 

c 


The error in this approximation is less than the accuracy of 
the doppler formulation; namely, 10' 5 m/sec. For purposes 


d# 3 (ET) 1 c ' 12 + 23 ' 


(44) 

(45) 
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In carrying out the derivatives in Eq. (43), the following quantities appear: 


*12 




P 12 


III 

+ 

3^*12 

0#! 

r 12 . 

— . r 

*12 

1 — » 2 
2 — ► 3 

df 12 

+ 

0fl2 

_ *i2 • r 12 + r 12 • r 12 “ r\ 2 

l-»2 

a 

CO 

0#! 

*12 

2— » 3 

dr 12 

_L 

0r'i2 

ri 2 • r*i 2 + 3ri 2 • ri 2 — 3f 12 r 12 

1-h>2 

= dt 2 

1 

0#! 

r, 2 

2-> 3 

r l2 




1 -> 2 

*12 

* i* 

L ) 


2— > 3 

_ dp i2 

+ 

0Pl2 

r 12 • r’i + ri * r 12 — p 12 f 12 

l-> 2 

— 0#2 

0#1 

fl2 

2-+3 

0P 12 

+ 

0P 12 

r 12 *r*i + 2 ri 2 • r'i + r x • r 12 — 2* 12 p 12 

- Pi 2 r'i2 

0f 2 

dU 

* 12 



1- 

2 - 


2 

■3 


Also, the following relations are used: 


0r 


12 


St, 

0fi 

0f» 

0#! 


= ~P 1 


= - Pit, 


P 12 > 


1 — » 2 

2-> 3 

1 -> 2 

2— > 3 


1 

2 - 


• 2 

■3 


(46) 

(47) 

(48) 

(49) 

(50) 

(51) 

(52) 

(53) 

(54) 


Differentiating r 12 and f 23 twice with respect to t 3 using the definitions in Eqs. (47) and (48), and using Eqs. (44), 
(45), (53), and (54), gives 


d 3 f 12 

dt'l 

d 2 r 23 

dt\ 


r,'o + (2 p\J l2 — 2r" 2 r' 2: , + KtPit ~ r t ir« 3 ) 

’f‘. 3 + — (2 p S af a3 + faaPzs) 
c 


(55) 

(56) 


For the derivatives in the 1/c 3 terms of Eq. (43), the derivatives (44) and (45) are taken to be unity, and 


df 12 

dt 3 

1^2 

~ r '-’ 2 — ► 3 

(57) 

dp u 

dt :l 

l-»2 

— Pit, 2 — > 3 

(58) 

d-p 12 

dtl 

l-»2 

~ p ‘ 2 ’ 2— > 3 

(59) 

Substituting Eqs. (55) through (59) into Eq. 43 gives 


(*-£) = K 

1 + y(Pia - 3f 2 .,)J + -i V'[l 4- “■ (p 2 :i - n=)* j 


+ 4 - 

C“ 

3 

()V 12 f 12 + p^h*) + — (r, 2 pi 2 4- r 23 p23 ^12^*23) 

(60) 
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where the quantities used are computed from Eqs. (46) 2. Pioneer VII Postlaunch Activities, j.e. Ball and 

through (51). J. f. Gallagher 


3. Computation of Two-Leg Range 

The two-leg range observable is computed from 


p — t 3 — t x + 


Rc 
10 3 c 


(61) 


where 


*3 — reception time of signal at tracking station, 
station time 


t ! = transmission time of signal at tracking station, 
station time 


Pioneer VII was launched on August 17, 1966, 15 h 20 m 
17. s 286 GMT. A nominal parking orbit was achieved and 
four AFETR tracking stations observed the coast phase. 
Transfer orbit injection occurred 25 m 21 8 from launch, 
and the heliocentric orbit was nominal. Three deep space 
stations (DSS 51, DSS 42, and DSS 11) tracked the space- 
craft. After approximately 36 hr in orbit, the spacecraft 
went through a Type-II orientation maneuver to prop- 
erly align its spin axis with respect to the ecliptic plane 
for communications purposes. Table 5 summarizes the 
results of the preorientation and the postorientation orbit 
computations. 


R c = solve for ranging bias, meters (specified by 
time block for each station) 

c = speed of light, km/sec 


The light time problem is solved in ephemeris time, and 
the round-trip time is converted to station time using 
the time transformations given in SPS 37-39, Vol. Ill, 
pp. 36—38. 


N67 15907 


F. Pioneer Project Support, o. W. Curkendall, J. E. Ball, 
and J. F. Gallagher 


1. Introduction, D. W. Curkendall 


The primary objective of the Pioneer VII mission was 
to send the spacecraft through the Earths magneto- 
sphere at a distance of several hundred Earth radii from 
the geocenter. On September 22, 1966, at about 12 h GMT, 
the Pioneer VII passed through syzygy with the Sun and 
Earth. The celestial longitude of the Earth at this time 
was 359.0383 deg, and the celestial longitude and lati- 
tude of the spacecraft were 359.04009 and 0.0555337 deg, 
respectively. The distance of the spacecraft from the 
Earth was 5,146,000 km, or about 807 Earth radii. The 
distance of the probe above the ecliptic plane was 150,504 
km, or about 23.6 Earth radii. This syzgy position satisfied 
the scientific needs of the magnetosphere experiment. 


The Pioneer program consists of a series of spacecraft 
designed to gather scientific information concerning the 
deep-space environment in the region near the ecliptic 
plane and at a distance from the Sun varying from 0.8 to 
1.2 AU. The DSN is committed to track the Pioneer space- 
craft, obtain engineering and scientific telemetry, deter- 
mine the orbit of each probe, and provide facilities to 
the Project within the SFOF and the DSIF for the con- 
duct of the mission and reduction of data. Specific items 
under this task will be selected during each report period 
for further discussion. 

In this issue, the successful launch and early orbit 
determination results of Pioneer VII are reported. 


The best orbit determination to date (October 26, 1966, 
postorientation 3) shows that there will be an occupation 
of the Pioneer VII by the Moon on January 20, 1967, at 
about 5 h 30 m GMT. The right ascension and declination 
of the Moon at this time are about 42.5 and 16.5 deg, 
respectively. This occupation will be visible to the Stan- 
ford antenna where the lunar occupation experiment will 
be carried out. A gas leak was recently discovered on 
the Pioneer VII. A cursory analysis of the perturbing ef- 
fect of the gas leak on the trajectory of the spacecraft 
reveals that although its existence will considerably de- 
grade the prediction of the probo position on January 20, 
an occupation should still occur. 
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Table 5. Summary of orbits 


Items 


Preorientation 

Postorientation 3 


Epoch 

15 h 45 m 38: 626(66/08/17) 

1 3 h 00 m oo: 0(66/08/19) 


x, km 

-6317.0086 

405987.78 


y, km 

-1695.2897 

-23031.865 

Injection conditions 

z, km 

-1689.9426 

-9448.2170 

x, km/sec 

3.4921793 

1.8976926 


y, km/sec 

7.7612659 

0.0266948 


z, km/sec 

-6.8737429 

0.0774747 


ax, km 

0.063235 

2.417940 


ay, km 

0.344739 

3.065523 


crz, km 

0.315043 

6.504306 

Parameter statistics 

ax, km/sec 

0.000533 

0.000003 


ay, km/sec 

0.000345 

0.000005 


crz, km/sec 

0.000811 

0.000014 


DSS 51 (Johannesburg) 
CC3 

504 

3050 

Number of data points in orbit 

DSS 42 (Tidbinbilla) 
CC3 

224 

3860 


DSS 1 1 (Pioneer) 
CC3 

183 

2640 


Epoch 

1 5 h 44 111 53:9(66/08/17) 

15 , ’48 m 08: 5(66/08/17) 


SMA 

-238677.91 

-241866.63 


ECC 

1.0282655 

1.0277063 


C3 

1.6700384 

1.6480210 

Geocentric conic" of pericenter passage 

RCA 

6746.3487 

6701.2334 


TA 

4.1508808 

160.14153 


INC 

41.924554 

42.170418 


LAN 

358.30088 

358.22291 


APF 

197.83713 

197.87560 


Epoch 

17 h 10" 1 53: 8 (66/05/28) 

00 h 00"’ 12: 5(66/07/18) 


SMA 

118127660.0 

162858590. 


ECC 

0.42326200 

0.07928469 


C3 

-1123.4804 

-814.90397 

Heliocentric conic 1 ' of pericenter passage 

RCA 

68128711.0 

149946400. 


TA 

148.19573 

32.939440 


INC 

7.8074922 

0.11492503 


LAN 

144.28250 

326.02138 


APF 

31.805896 

327.21709 

"Osculating elements of injection epoch. 

Definition of terms: 

SMA— semi-major axis TA 

ECC— eccentricity IN 

C3— energy LA 

CC3— two-way doppler data AP 

RCA— radius of closest approach 

—true anomaly 
C— inclination 

N— longitude of ascending node 
F— argument of perifocus 
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III. Communications Research 
and Development 


N67 1590® 

A. Frequency Generation and Control: Atomic Hydrogen Frequency 
Standard, w. H. Higa 

The dual channel traveling wave maser (TWM) has been designed and is pres- 
ently being integrated into a closed-cycle refrigerator (CCR). Each section will 
conform to the following specifications: 

f = 1420 GHz 
3-db BW = 10 MHz 
Gain ^ 25 db 
Noise temperature — 10°K 

A simple test configuration will be used to evaluate the stability characteristics 
of the TWM. Subsequently, the atomic hydrogen masers will provide separate 
signals for the two channels of the dual TWM. Well known signal processing 
techniques (Ref. 1) may then be used to derive the power spectra of the maser 
oscillators. 

Fig. 1 shows the dual channel maser installed in the CCR (SPS 37-39, Vol. Ill, 
pp. 79-81). Fig. 2 shows the top view of the assembly. 




i 


i 
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WAVEGUIDE 


DUAL 

CHANNEL TWM- 


OUTPUT 

CHANNEL 


CHANNEL I 


Fig. 1 . Dual channel maser installed in the 
closed cycle refrigerator 


B. Improved RF CaDb&tfcn T^I§i$<Q fj^aily 
System Noise Temperature Measurements, 

C. T. Stelzried 

The daily noise temperature calibration technique used 
at the Mars and Venus stations has been discussed in 
detail (SPS 37-39, Vol. Ill, p. 96). Noise temperature cali- 
brations are made by switching the maser input between 
the antenna and an ambient load with a waveguide 
switch. The principal advantage of this technique is the 
long term stability and reliability of the ambient wave- 
guide termination. This technique depends on a previous 
calibration of the maser/follow-up amplifier effective 
noise temperature. A 1°K error in this calibration results 
in approximately a 0.1 °K error in the system temperature 
calibration. A refinement has been made consisting of a 
follow-up temperature contribution measurement for each 
calibration with a maser off-on Y-f actor measurement. 



PUMP 

WAVEGUIDE 


OUTPUT 
CHANNEL I 


OUTPUT 
CHANNEL 2 


INPUT \ 
CHANNEL I 


IANNSL 


Fig. 2. Top view of the maser/CCR assembly 


The follow-up contribution defined at the maser input 
measured by turning the maser off with the input con- 
nected to the ambient termination is given to good approx- 
imation by 


T f = 


T m + To 


Yon - 1 


( 1 ) 


where 


T m = maser noise temperature defined at the maser 
input, °K 

T () = ambient termination temperature, °K 

Y 00 — receiver output noise ratio between maser on 
and off 


The receiver noise temperature defined at the maser input 
is then given by 


T/f — T M 4- T f 


( 2 ) 


A revised calibration sheet, 1620 computer program 
(modified by Lois Busch), and sample print-out are shown 
in Figs. 3-5. 


JPL SPACE PROGRAMS SUMMARY 37-42 , VOL. Ill 


25 


a 

w 


Program Return a Copy to C. Stelzried 

CTS/20 B JPL 238-737, Ext. 3838 


RECEIVING SYSTEM NOISE TEMPERATURE DAILY CALIBRATION 


Operators Name 

lL 

. 



,K i A i 

6/ 

Mission 

3 

iM. 

A 


/ 

,/v.r 

/? / V 

W eather 

2T 

lL 

L 


A 



F requency 

43 

lL 

2 

9 , 

7 

, • ,6 

s 1 


63- 


Station No. | / 4- • I Month I 0 . 9 . • | Day | 0 . 8 . • | Y ear L / . ^ ^ ^ 

3 6 10 13 

Calibration. Time: hour (GMT) 

Day No. | 5 / . • | (to be recorded half way thru | O , / . • | Minute s | 0 . O . • | 

18 Y-factor measurements) 22 25 

Antenna Coordinates: AZ | £ . 7 . 6 . • . 6 . £ 1 EL | 8 „ 7 . • . 0 . 5 \ 

28 34 


l_ 

39 

Comments 

1 L 

65 


Signal Generator on Freq. lL 

I 

Pre-Cal Uj Post Cal lZj 

4 5 

Maser Gain GM (db) | 4 . 0 . • . 0 \ 
7 


Detector Bias, 4 mv LL Filter, Normal | 2 \ 

2 3 

Other | / | 

Maser Gain GM(dB) if reset or peaked 

(Do not adjust on Post Cal or if gain | . . • . I 

down less than 1 dB from nominal) 13 


Reflectometer Measurements: Ambient Load 

lL 

7 , 

,fl| 

db 


20 




Antenna 

lL 

/ 


db 


25 


Maser on-off, Yoo (db) \2 . S . • . 6 | 

30 

Ambient Temperature, T rt ~,| 2 ^5 ^ • . / . / | 

o c 


(Switch to ambient load. Measure before 
noise. Y-factors on pre-calibration and 
after noise. Y-factors on post-calibration . 

deg C. (To be recorded half way through 
Y-factor measurements) 


Maser Temperature, T^ , | 7 • 2 \ deg K. 

M 45 


Noise Temperature (AIL Readings, dB) 


AMB 

A 0 (I) 


ANT 

a a (D 


L 2 . 0 . .J 9 . , 

1 20 . 1 7 1 

1 .2.0 .18 1 

1 20,./ 6 1 

1 20 . 18 . 1 

3 

13 

23 

33 

43 

1 9.8f 1 

1 9.83 1 

1 . 9.04 1 

1 9.83 1 

1 9 . 83 . 1 


3 13 23 33 43 


Fig. 3. Daily system temperature calibration data sheet 
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C RECEIVING SYSTEM NOISE TEMPERATURE DAILY CALIBkAT IUiviS 

C PROGRAM FOR DATA ONLY (NO PLOTS) CTS/20B 

DIMENSION AODBI 5) , AADBI 5), YOADBI 5) , YOA( 5) , TSA ( 5 ) , OP ( 20 ) , MI SSI 20 ) , 

1 WTHR( 20 ) ,FRE0(7 ) ,COMM( 40) , ND( 50) , TP( 50) ,TS (50) ,TM(50) ,F (50) , 

2 AM (50 ) , AN ( 50 ) ,C,( 50 ) , D( 50) 

COMMON M 

M = 0 

IF ( SENSE SWITCH 9)40,40 

40 READ 1, OP, MISS, WTHR, FREQ 

READ 2, NO, MON, I DAY, IYR,NDAY, I HR, M I N, AZ , EL »COMM 

READ 3 , 1 SIG, I DET, IF I L, I PRE, I POST, I OTHER , GM ,GMR , AMBL , AiMT , YOO, T02 , 

1 TEM 
DAY =ND AY 
HR = IHR 
DM I N =M I N 
M=M + 1 
J = M 

T0K=T02+273. 16 

TF= (TOK + TEM ) /( 10 .-■=*( YOO / 10 . ) - 1 . 0 ) 

TR=TEM+TF 

READ 4, ( AOD B ( I ) , 1 = 1, 5 ) 

READ 4, ( A AD B ( I ) , 1 = 1, 5) 

39 PRINT 71 
PRINT 72, NO 
PRINT 73, OP 

PRINT 74,M0N,IDAY,IYR.NDAY 
PRINT 75, IHR, MIN 
PRINT 76 , FREQ 
IF ( I OTHER ) 34, 34, 36 
34 PRINT 93, IPRE, IPOST 
GO TO 38 

36 PRINT 92, IPRE, IPOST, IOTHER 
38 PRINT 77 , W THR 
IF ( AZ )5 1,52,51 

51 IF(EL)53»54,53 

52 IF(FL)56»55,56 

54 PRINT 23, AZ 
GO TO 58 

55 PRINT 25 
GO TO 58 

56 PRINT 24, EL 
GO TO 58 

53 PRINT 2 1 , AZ , E L 
58 I F ( GMR ) 42,42, 41 

41 PRINT 18 ,GM ,GMR 
PRINT 22, AMBL, ANT 
GO TO 43 

42 PRINT 19, CM 
PRINT 22, AMBL, ANT 

43 PRINT 82, YOO 
PRINT 83 , TO 2 
PRINT 84, TEM 
PRINT 85, COMM 
PRINT 11 

44 DO 50 1=1,5 
N = I - 1 

IF ( AOD B ( I ) >50,70,50 

Fig. 4. Daily system temperature IBM 1620 computer program 
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50 CONTINUE 
N = 5 

70 PETO 2 = 0 . 1 
PETR= 1. 0 
A=0 .003 
B=0. 00354 
EN = N 
SUM=0 .0 
SV0A=0.0 
DO 100 I = 1 , N 

YOADBl I ) = AO DB ( I ) -AADB ( I ) 

YO A ( I )= 10.** (YOADBl I 1/10. ) 

TSAI I ) = ( TOK + TR ) / YO A I I ) 

SUM=TSA I I l+SUM 
SYO A = YO A I I l+SYOA 

99 PRINT 1 2 , AODB ( I ) , AADB (I ) , YOADtt ( I > , TSA ( I ) 
100 CONTINUE 

STSA=SUM/EN 

YOAS=SYUA/EN 

DBYOA=LOGF I YOAS ) *4. 3429448 


S T T = 0 . 0 
DO 200 I = 1 1 N 

STT= ( STSA-TSA( I ))**2+STT 


200 CONTINUE 

PEMI=0.6745*S0RTF I STT/I EN-1. ) ) 

PEM=PEM I /SORTF I EN) 

235 PE YQ = Y0AS**2* ( A**2+B**2* DB YOA** 2 ) * • 05301 898 1 

PET SQ = P EM**2 + ( PET02**2+PETR**2)*( 1 . / YOAS ) **2+PE YO* ( (TUK + TR )**2/ 
1 YOAS* *4 ) 

PETSA=SORTF I PETSO) 

TPS= ST SA + PE TSA 
TPM=STSA-PETSA 


238 PRINT 86 , STSA 
PRINT 87, TR 
PRINT 88, TF 


PRINT 8 9 , P F M I 
PRINT 90 , PEM 
PRINT 9 1 , PE TSA 


245 IF I SENSE SWITCH 9)255,40 
255 CONTINUE 

1 FORMAT I 2X , 31 20A1 ) , 7A1 ) 

2 FORMAT (14, 1X12, 2X12, 1X14, 1X13, 2(1X12), 1XF6. 2, F5. 2, 40 Al) 

3 FORMAT (611, F 4. 1 , F 6. 1 , F 7 . 1 , 2 F 5 . 1 , F6. 2 , F 9 . 1 ) 

4 FORMAT ( 5F 10 .0 ) 


71 FORMAT ( 1H1,8X40HSYSTEM NOISE TEMPERATURE PROGRAM CTS/20B) 

72 FORMAT ( 1H ,8X1 1HSTAT I ON MO. 13) 

73 FORMATdH , 8X9H0PERAT0R 20A1) 


74 FORMAT ( 1H 

75 FORMATdH 

76 FORMAT ( 1H 

77 FORMATdH 
2 1 FORMAT ( 1H 
25 FORMATdH 
24 FORMAT ( 1H 
23 FORMATdH 
18 FORMAT! 1H 

1 F5.1) 


, 8X 5H DAT E I 2 , 1H/ I 2 , 1 H/ I 4, 10H uAY.NU.I4) 

,PX9HTIME(GMT) 13, 6H HUURS13,8H MINUTES) 

» RX 14HFR EQUENC Y( MC ) 7 A 1 ) 

, 8X7HWEATHER , 1X20A1 ) 

, 8X24H ANT E NNA COORD. (DEG.) AZ . F 7 . 2 , 2 X3HE L . E6 . 2 ) 

, RX 24HANTENNA COORD. (DEG. ) AZ . » 9X » 3HEL . ) 

, RX24HANT ENNA COORD. (DEG.) AZ. ,9X ,3HtL.F6.2 ) 

» 8X 2 4H ANTENNA COORD. (DEG.) AZ . F 7 . 2 , 2 X 3HEL . ) 

, R X 1 4HMA S ER GAI N< DB ) , F6. 2/ 1H ,8X21HPEAKED mASFR GAIN (OB ) 


Fig. 4 (contd) 
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19 FORMAT ( 1H , 8X 14HMASER GAIN( DB) , F6.2/1H ,8X23HPEAKED MASER GAIN(DB) 

___1__J 

22 FO RM AT ( 1H ,8X2 3~HR~EF L ME A~S . (DB ) AMB. LOAD , F~5 ."l , YX4HANT^F'5 . 1 ) 

§2. _ F_Q RM A_T j A H_ _ , j3 X JL 6 H. M A S £_R __0 N - 0_F F ( DB 1_,_F_5_.1A 

83 FORMAT ( 1H , 8X20HAMB I ENT f EMP. ( DEG. C ) , F6".~2 )" " “ 

84 FORMAT ( 1H , 8X18HMASER TEMP. (DEG.K ). F5. 1) 

85 FORMAT ( 1H , 8X9HC0MMENTS 3X30A1,/1H ,8X10A1) 

1 1 _ f. 0 RM AT J 1 HO ,_8X_5H A0 (I),4X5H AA < I )_»_4X6 HYQA( 1 ) , 4X6 HJ S_A_(_I 

12 F"6rMA"T(7XF7.2,2F9.2,F10.2) ' ' “ ““ 

86F0RMATJ 1HO , 8 X_1 9JH SYS JEM TEMP . ( DEG.K ) , F9. 2) 

87 F'ORMAf ( 1H , 8X 2 1HR EC ¥f V ER“T EMTi ( DEG'.'K ) , F7".'2 j 

88 FORMAT ( 1H , 8X22HFOLLOW-UP CONT. ( DEG. K ) , F6 .2 ) 

89 FORMAT ( 1H , 8X 39HPRO B. ERROR SYSTEM TEMP. I ND .MEAS . ( DEG .K ) , F7 . 2 ) 

90 FORMAT ( 1H , 8X35HPR0 B. ER ROR SYSTEM TEMP. MEA S. ( DEG.K ) , F 1 1. 2) 

91 FORMAT (lH , 8X38 HP ROB. ERROR SYSTEM TEMP70VE“RALL (DEG.K) ,F8.2) " 

92 FORMAT ( 1H ,8X25HTEST CONDITIONS PRE-CAL , I 2 , 2 X8HP0ST-C AL , I 2 , 5X 

"“l“I7HR0W-5TWWiny;T2T 

93 FORMAT ( 1H ,8X25HTEST CONDITIONS PRE-CAL , I 2 ,2X8HP0ST-CAL , I 2 , 5X 

— 1 12hN0N-STAKiDArD) 

END 


Fig. 4 (contd) 


SYSTEM NOISE TEMPERATURE PROGRAM CTS/20B 
STATION NO. 14 

OPEKATUR N M KAKU 

DATE 9/ 8/1966 DAY NO. 251 " 

TIME (GMT ) 1 HOURS _J)_ MINUTES _ 

FREUUENCY ( MC) 2297.6 

TEST CONDITIONS PRE-CAL 1 POST-CAL 2 _ NON-STANDARD 1 

WEATHER CLEAR 

ANTENNA COORD. (DEG.) AZ . 276.62 EL. 87.05 

MASER GAIN(DB) 40.00 

PEAKED^ MASER GAIN(DB) 

RE FL .ME AS. ( DB ) AMB. LOAD 39.8 ANT. 31.4 

MASE R ON-O FF ( OB ) _25^6 

AMBIENT TEMP. (DEG. C) 25. 11 

MASER TEMP. (DEG.K) 7.2 

COMMENTS 


AO ( I ) A A ( IJ YOA ( IJ TS A ( I ) 

20. 19 9.84 10.35 28.25" 

20 . 17 9.83 10.34 28.3 2_ 

20.18 9.84 10.34 28.32 

,20^18 9.83 10.35 _ 28.25 

20.18 9.83 10.35 28.25 

SYST E fTTEMP. f DEG .K) 2 8.2 8 

RECEIVER TEMP. (DEG.K) 8,0 4 

FOIL OW-lIP CONT. (DEG.K) .84 


PROB. ERROR _SYSTEM_TEMP . IND.MEAS. ( DEG. K ) _.02 
PR OB. ERROR SYSTEM TEMP. ME AS. (DEG.K) .01 
PROB. ERROR SYSTEM TEMP. OVERALL (DEG.K) .25 


Fig. 5. Daily system temperature print-out 
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Day 

No., 

1966 

Ts, 

°K 

Tf, 

°K 

Cm, 

db 

Zen 


Test conditions 

Clear 

weather 

Frequency, 

MHz 

Comments 


Pretest 

Posttest 

Other 

203 

21.3 

0.34 

41.2 

X 


X 


X 

2388 


203 

21.0 

0.37 

45.1 



X 








204 

22.3 

0.32 

43.8 

i 











205 

22.41 

5.02 

43.3 

? 





f 





206 

22.4 

5.23? 

45.2 

> 






X 




210 

22.5 

0.7 

40.2 







X 




210 

21.3 

0.6 

39.8 




f 


X 

No 




211 

25.0 

0.7 

40.4 



X 








213 

21.0 

0.5 

40.4 




? 







214 

22.1 

0.4 

39.1 







t 




261 

22.2 

0.3 

38.0 







X 




264 

22.2 

0.4 

39.9 






X 

X 




262 

21.7 

0.3 

37.0 







No 




263 

21.1 



41.0 







X 




265 

22.4 



38.4 






X 






266 

22.9 



40.8 



X 








266 

22.0 



37.2 





X 






267 

22.6 



38.2 












268 

22.0 



37.6 












269 

21.1 



40.1 



X 


1 


1 





270 

21.9 



40.5 





X 

No 




271 

21.0 

1 

f 

40.0 



X 



x 




272 

20.9 

0.4 

40.9 





X 

No 




272 

34.9 

3.0 

34.5 



X 







Rain on horn 

273 

21.6 

0.3 

40.1 




0 







274 

21.3 

0.4 

40.1 



X 




f 




275 

21.7 

0.3 

37.6 

1 

r 





X 




276 

22.0 

0.4 

37.4 

No 

1 

r 








276 

24.7 

0.7 

38.9 

X 



X 






277 

21.7 

0.4 

37.5 








1 

f 


277 

27.1 

0.6 

40.0 






1 

f 

2297.6 

Mariner IV 

279 

20.8 

0.4 

37.3 



X 



No 

2388 


279 

24.5 

0.1 

48.8 



X 




< 

2297.6 

Mariner IV 

279 

24.8 

0.1 

52.0 




X 




2297.6 

Mariner IV 

280 

21.6 

0.5 

37.5 






< 



2388 


282 

21.3 

0.3 

37.6 












283 

22.3 

0.6 

39.2 









1 

L 


284 

24.5 

0.1 

51.5 



X 






2297.6 

Mariner IV 

286 

21.9 

0.5 

39.7 



X 


1 

f 



2388 


287 

22.8 



39.8 






( 






288 

22.6 



39.7 












290 

23.7 

1 


39.0 










r 


291 

24.2 

0.1 

53.0 









2297.6 


292 

23.0 

0.5 

39.0 









2388 


293 

23.8 



39.0 



X 


1 

r 






293 

23.1 

1 

r 

39.0 






K 






294 

22.1 

0.4 

39.7 






X 






295 

20.3 

0.5 

39.7 



X 








296 

21.9 

0.4 

39.1 






X 






297 

22.1 

0.5 

39.0 

1 

r 




X 


r 


? 



Fig. 6. Tabulated results of the Venus station daily noise temperature calibrations 
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The principal results of the daily noise temperature 
calibration for the period July 12 (day No. 193) to Octo- 
ber 24 (day No. 297), 1966 are tabulated in Figs. 6 and 7 
for the Venus and Mars stations. 1 It is suspected that the 


'These data were taken by the Venus station personnel. 


high followup contributions on days 194 and 234 at the 
Mars station and days 205 and 206 at the Venus station 
were due to making a maser on-off Y-factor measurement 
with the antenna switched to the antenna instead of the 
ambient termination, as required for proper data reduc- 
tion. The average system noise temperatures and reported 
maser gains for this period are tabulated in Table 1. 


Day 

No., 

1966 

Tx, 

°K 

B 

Gm, 

db 




»ar 

ther 

Frequency, 

MHz 

Comments 



E 

193 

26.1 

0.4 

40.1 






X 

2338 


194 

27.2 

0.4 

40.0 





X 






194 

29.0 

4.1 

40.0 

1 



X 







201 

45.5 

0.6 

40.0 

No 


X 







202 

26.2 

0.7 

40.4 


X 



X 


t 




204 

24.7 

0.7 

42.5 




X 


No 




208 

27.1 

0.6 

41.0 

? 


X 


No 




209 

27.3 



40.5 


< 



X 


X 




209 

28.8 



40.5 




X 







210 

45.3 

1 

f 

40.2 





X 

1 




RF spur 

213 

46.3 

0.9 

38.6 




X 


No 



RF spur 

232 

28.3 

0.6 

40.4 




X 



X 

2297.6 


233 

27.9 

0.6 

41.3 





X 






234 

30.7 

24? 

39.8 




X 







234 

28.3 

0.6 

41.3 











244 

29.5 

0.8 

40.4 











244 

28.8 


40.1 

1 

f 


f 

X 






250 

34.4 

0.9 

40.0 

No 


X 







250 

28.1 

0.8 

40.0 

X 

X 







251 

28.3 

0.8 

40.0 





X 

i 

f 




252 

26.5 

0.9 

38.6 



X 


No 




253 

24.8 


38.3 





X 



1 

r 


270 

26.5 

0.8 

39.0 



X 




2388 


271 

27.6 

0.8 

38.6 





X 





RF spur 

278 

43.5 

0.4 

39.0 

1 

r 



X 



1 

f 


283 

31.7 

0.8 

40.6 

No 

X 




2297.6 


283 

68.0 

0.7 

40.8 

X 



X 






284 

36.3 

0.8 

40.5 

No 


X 






2-way diplex mode 

285 

34.5 

0.7 

40.5 

X 



X 

1 

f 


f 

2-way diplex mode 

285 

33.8 

0.9 

38.8 

X 

X 



< 

2388 


285 

29.7 

0.8 

40.4 

No 






2297.6 


285 

34.7 

0.7 

39.5 

No 






2388 


286 

28.6 

0.8 

40.0 

X 






2297.6 


286 

27.7 

0.7 

38.8 

) 


1 

f 

X 



2388 


287 

37.7 

0.9 

38.6 

? 

X 


1 

r 

2388 


287 

28.9 

0.8 

40.2 

) 




X 

No 



288 

28.2 

0.8 

39.6 





X 

X 

2388 


290 

28.9 

0.7 

40.0 



X 




2297.6 


291 

31.9 

0.8 

39.4 



> 








291 

28.6 

0.9 

39.6 





X 






291 

35.0 

0.9 

40.0 





X 





2-way diplex mode 

292 

27.7 

0.8 

39.9 



> 



1 





293 

27.5 

0.8 

39.9 



> 



No 




294 

27.9 

0.9 

40.0 





X 

No 

1 



294 

25.5 

0.7 

39.3 



> 



X 

2388 


294 

26.2 

1.0 

37.6 

1 

f 



X 

No 

2388 



Fig. 7. Tabulated results of the Mars station daily noise temperature calibrations 
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Table 1. Summary of the system noise temperature 
and maser gain averages 


Station 

Frequency, 

MHz 

System temperature, 
°K 

Maser gain, 
db 

Venus 

2297.6 

25.0 ± 1.2 (l-o variation) 
5 data points 

51.3 ± 1.8 (l-o variation) 
4 data points 
(excluding day No. 277) 

2388 

22.2 ±0.9 (l-o variation) 
32 data points 

39.5 ± 2.0 (1 -o variation) 
45 data points 

Mars 

2297.6 

28.9 ± 1.2 (1-a variation) 
13 data points 

40.2 ±0.5 (l-o variation) 
23 data points 

2388 

28.1 ± 2.4 (1-a variation) 
9 data points 

39.5 ±1.1 (l-o variation) 
23 data points 


This is the first reported use of the Venus station R & D 
cone at 2297.6 MHz, since installation of the modified 
maser with broadband frequency tuning capability 
(SPS 37-40, Vol. Ill, p. 90). The system temperature aver- 
ages were computed using only data which were acquired 
as follows: 

(1) Antenna at zenith. 

(2) Clear weather. 

(3) No RF spur in receiver pass band. 

(4) Not in two way diplex mode. 

The maser gain averages were computed using all data 
points, except as noted. The maser gain was optimized 
(usually with magnet current or pump frequency), if the 
gain had dropped more than 1 db below nominal during 
this period. The maser gain variations have some effect 
on system temperature variations. The next series of mea- 
surements will be made to optimize the maser gain during 
each precalibration. 

N67 15910 

C. Improved RF Calibration Techniques: Weekly 
Noise Temperature Calibrations, T Y. Otoshi 

1. Calibration Technique 

By the method described for daily calibrations in sub- 
section B of this section, it was shown that system tem- 
perature can be determined to good accuracy using a 
single thermal noise standard (ambient load), even 
though the receiver temperature might be known only 
approximately. 

To obtain an absolute calibration of receiver tempera- 
ture, however, or to separate system temperature into its 


antenna and receiver temperature components, it is re- 
quired that a second thermal noise standard (cryogenic 
load) be used. The Y-factor calibration technique involv- 
ing two thermal noise standards was discussed in detail 
in SPS 37-39, Vol. Ill, pp. 86-91, and also in SPS 37-40, 
Vol. IV, pp. 202-209. 

For reference purposes the block diagram for the basic 
measurement system is shown in Fig. 8. The double- 
primed symbols indicate temperatures defined at the 
maser input reference flange; all single-primed symbols 
refer to temperatures defined at the antenna ouput refer- 
ence flange. The double- and single-primed temperatures 
are related by the dissipative losses of the interconnecting 
transmission lines. 

The notations used are defined as follows (the primes 
have been dropped for generality): 

T c = cryogenic load effective temperature, °K 

T r = receiver temperature, °K (includes the contribu- 
tion of the receiver system following the maser) 

X A = antenna temperature, °K 

To 2 = physical and effective noise temperature of the 
ambient load, °K. (The physical temperatures 
of the transmission lines are also assumed to be 
at temperature T 02 ) 

T e = excess noise temperature of the gas tube injected 
into the maser amplifier, °K 

For this measurement system, the expression for system 
temperature, as defined at the maser input reference 
flange when switched to the antenna, is given as 


Expressing the antenna and receiver temperatures in 
terms of their respective temperature components gives 

X" - T' t + T l (2) 

and 

X" - T" 4- X'; (3) 
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c I = O.OOI352°K/mm Hg 


&P - BAROMETRIC PRESSURE 

DIFFERENCE FROM 760 mm Hg 

Fig. 8. Block diagram of basic noise temperature measurement system 


where 

T. \ = antenna temperature defined at the antenna out- 
put reference flange, °K 

T l = temperature contributed by the dissipative loss 
of the transmission line between the antenna out- 
put and maser input reference flanges, °K 

TZ = maser temperature defined at the maser input 
reference flange, °K 

T p — follow-up receiver temperature contribution de- 
fined at the maser input reference flange, °K 

Substitution of Eqs. (2) and (3) into Eq. (1) gives the 
following expression of antenna system temperature: 

T" = T' a + T l + T" + T" (4) 


The results of weekly noise temperature calibrations in 
the following section will be presented in terms of the 
components of system temperature given in Eq. (4). 

2. Calibration Results 

Periodic calibrations of the absolute receiver and an- 
tenna temperatures as well as calibrations of gas tube 
excess temperatures for the Mars Deep Space Station 
(DSS) and Venus DSS antenna systems have been con- 
tinuing on a weekly basis. 2 Earlier weekly calibration 
work on these antenna systems has been reported in 
SPS 37-39, Vol. Ill, pp. 92-96 (for Mars DSS), and 
SPS 37-41, Vol. Ill, pp. 82-86 (for Venus DSS). 

The antenna system being calibrated at Mars DSS con- 
sists of the 210-ft diameter paraboloid and a multifre- 
quency cone (formerly called the listening cone). At 
Venus DSS, calibrations are being done on the antenna 


'This work is being done by Venus DSS personnel. 
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system consisting of the 85-ft diameter paraboloid and 
the R&D cone. The type of cryogenic thermal noise stan- 
dard being used to calibrate both of these antenna sys- 
tems is a liquid nitrogen-cooled WR 430 waveguide load. 

A plot of all weekly noise temperatures calibrated to 
date may be seen in Figs. 9 and 10 for the Mars and 
Venus DSS antenna systems, respectively. Each value 
represents an average obtained from 10 sets of Y-factor 
measurements. These figures show that problems are be- 
ing encountered in these weekly calibrations. For instance, 
Fig. 9 shows that recent calibrated values have begun to 
fluctuate excessively. The Mars DSS calibrated antenna 
temperatures, as defined at the antenna flange, have been 
varying between maximum and minimum values of 12.5 
and 3.4°K, while the calibrated maser temperature values 
have been fluctuating between 17.5 and 3.6°K. Antenna 
temperatures of 12.8 to 4.5°K and maser temperatures of 
10.4 to 6.3° K have been calibrated for the Venus DSS 
antenna system. 

Although the calibrated temperature fluctuations are 
not as severe for the Venus DSS antenna system, definite 


cr 

ZD 

<t 

cr 

LU 

CL 

5 

LU 


CO 

O 



80 100 120 140 160 180 200 220 240 260 280 300 320 

1966 GMT DAY NUMBER 


Fig. 9. Mars DSS advanced antenna system 
weekly calibrated noise temperatures 


calibration problems appear to exist for both stations. 

It is suspected that the major cause of the inconsistencies 
in the calibrated data is the instability of the liquid 
nitrogen-cooled loads. This source of error is suspected 
because there does not appear to be similar types of 
fluctuations existing in the daily calibrations. A study has 
been in progress to determine several possible causes of 
errors in the weekly temperature calibrations. Some of the 
problems being investigated are: 

(1) Are the measured Y-factor ratios dependent upon 
the detector bias voltage setting? 

(2) What is the effect of not venting the cryogenic load 
dewar to the atmosphere? 

(3) Is the termination element of the liquid nitrogen- 
cooled load actually at the cryogenic liquid tem- 
perature? 

(4) Are the temperature fluctuations due to changing 
mismatch conditions of the thermal noise standards, 
antenna, or receiver? 

(5) Are the weekly calibrations being affected by sys- 
tem bias errors (such as external spurious signals 
or other noise contributions) not present in daily 
calibrations? 

Special tests have been made* to determine the effects 
of (1) changing detector bias voltage settings, (2) pres- 
surizing and not pressurizing the liquid nitrogen load 
with helium gas, and (3) venting and not venting the 
cryogenic dewar to the atmosphere. A 1-psi increase in 
pressure on the liquid nitrogen fluid will raise the liquid 
nitrogen temperature by approximately 0.6°K. The re- 
sults of the tests are summarized in Tables 2 and 3, and 
the effects of the various operating conditions on the 
calibrated temperatures can be noted. Each temperature 
value shown in the tables is the average obtained from 
10 sets of Y-factor measurements. These special tests have 
established that the detector bias should be set at approxi- 
mately 4.0 mv for linearity reasons, and the cryogenic 
dewar should always be vented to the atmosphere. The 
various problem areas will continue to be investigated. 

Even with optimum operating conditions, a liquid 
nitrogen-cooled load may not provide the resolution and 
accuracy needed for the weekly calibrations. An analysis 
of probable errors and mismatches indicates that for cali- 
bration of low-noise receiver and antenna temperatures, 
it would be more accurate to use a waveguide liquid 
helium-cooled load rather than a liquid nitrogen-cooled 
load. However, a satisfactory field-type liquid helium 
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Fig. 10. Venus DSS antenna system weekly calibrated noise temperatures: (a) antenna and system noise 
temperatures are via Arm A path only; (b) antenna and system noise temperatures via Arm B path 
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Table 2. Special noise temperature calibration tests on the Mars DSS 
advanced antenna system at 2295 MHz 


Part A. Detector bias voltage tests 

Description of test 

1966 

GMT 

day 

No. 

GMT, 

hour 

and 

minutes 

T.i via 
RCP 

path, °K 

T.J, 

°K 

Tf, 

°K 

_ " 

Ini via 
RCP 

path, °K 

T EL, 

°K 

Teh, 

°K 

Maser 

Reflectometer readings, 
db 

gain, 

db 

Antenna 

via 

RCP path 

Ambient 

load 

LN 

load 












39.6 

38.3 




10.83 








41.1 

38.2 

Less sensitive detector used, 
2-mv bias 

159 

2030 

11.37 








41.2 

38.3 








LN load pressurized with He gas 
{normal procedure) 

261 










33.9 

41.0 

LN load not pressurized 

261 

0530 

10.90 

7.35 

0.66 

28.27 

4.18 

33.48 

1 40.2 

30.7 

33.9 

43.4 


Part C. Liquid nitrogen dewar pressurization tests 





LN dewar vented through slit in 
rubber tubing 

282 

0615 

12.49 

3.02 

0.78 

25.51 

4.27 

31.93 

40.3 

30.0 

33.7 

41.8 

LN dewar vent tube open to the 
atmosphere 

282 


8.70 

9.61 

mi 

28.16 

H 

34.82 

41.7 

30.0 

33.4 

41.7 

LN dewar vented through slit in 
rubber tubing 

289 

ion 

13.34 

3.59 


27.09 


32.48 

40.0 

31.0 

40.8 

41.8 

LN dewar vent tube open to the 
atmosphere 

Q 

jjHj 

11.67 

5.48 



■ 



31.4 

40.9 

. 

40.6 

* Reflectometer was retuned on day No. 



Table 3. Special noise temperature calibration tests on the Venus DSS 
antenna system at 2388 MHz 


Liquid nitrogen dewar pressurization tests 

Description of test 

1966 

GMT 

day 

No. 

GMT, 

hour 

and 

minutes 

Ta via 
Arm A 
path, °K 

ii. 

°K 

n. 

°K 

tt 

Tsa via 
Arm A 
path, °K 

TkV, 

°K 

Teh, 

°K 

Maser 

gain, 

db 

Reflectometer 
readings, db 

Antenna 

via 

Arm A 
path 

Ambient 

load 

LN 

load 

LN dewar vented through slit in 
rubber tubing 

283 

0402 

3.35 

10.39 

0.51 

22.66 

4.24 

45.8 

38.8 

35.7 

48.4 

53.4 

LN dewar vent tube open to the 
atmosphere 

283 

0429 

4.53 

9.57 

0.51 

22.99 

4.31 

46.4 

38.5 

35.3 

47.3 

56.9 

LN dewar vented through slit in 
rubber tubing 

290 

1000 

5.11 

8.07 

0.46 

21.81 

4.23 

44.5 

38.5 

36.0 

44.8 

38.7 

LN dewar vent tube open to the 
atmosphere 

290 

1056 

6.13 

9.67 

0.46 

23.92 

4.22 

44.2 

38.8 

36.0 

41.4 

42.0 
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waveguide load has not yet been developed. It is sug- 
gested that a good operational liquid helium load be 
developed and used on a trial basis for absolute tem- 
perature calibrations now being performed on the antenna 
systems at Goldstone. 
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D. Efficient Antenna Systems: Calculated Gain 
of the Advanced Antenna System, A. c. Ludwig 

The gain of a 210-ft diameter antenna with uniform 
aperture illumination at 2295 MHz is 63.75 db. Expected 
gain of the real antenna is computed by evaluating the 
effect of deviations from the ideal antenna. Losses occur 
due to (1) nonuniform aperture illumination and spill- 
over, (2) aperture shadowing or blockage by the sub- 
reflector and its support structure (quadripod), and (3) 
deviations of the reflecting surfaces from a perfect parabo- 
loid or hyperboloid, including misalignment or defocus- 
ing of the feed and subreflector. 

1. Nonuniform Aperture Illumination and Spillover 

The primary feed for this antenna is a modified dual 
mode horn (Ref. 2 and SPS 37-36, Vol. IV, pp. 249-251). 



Fig. 1 1 . Primary feed pattern 


Experimental principal plane patterns, which are suffi- 
cient to completely determine the radiation of this type 
of feed (SPS 37-26, Vol. IV, pp. 200-208), are shown in 
Fig. 11. The phase patterns, which are not shown, are 
virtually flat over the region of interest. The subreflector, 
which consists of a hyperboloid with a conical flange 
(Ref. 3), subtends ±14° 39.3' of the feed pattern, resulting 
in 7.44% forward spillover. 

The scattered pattern from the subreflector (main re- 
flector illumination) is shown in Fig. 12. This pattern was 
computed numerically using a computer program (Ref. 4) 
that has demonstrated outstanding agreement with experi- 
mental results for similar situations.* The illumination 
efficiency of this pattern was evaluated with a second 
computer program (Ref. 4) that solves an exact integral 
for efficiency, and also divides the net loss into several 
components, some of which are approximate (SPS 37-26, 
Vol. IV, pp. 200-208). These losses 4 are shown in Fig. 13 
along with other losses, discussed shortly. The net loss 
due to nonuniform aperture illumination and spillover 
past the subreflector and main reflector is 1.15 db. 

2. Aperture Blockage 

The aperture blockage due to the subreflector and 
quadripod is illustrated in Fig. 14. (The antenna shown is 
a scale model of the 210-ft diameter antenna.) The out- 
line of the structure blocks 7 l A% of the aperture. However, 
since the quadripod is not a solid structure, one may 
assume an opacity factor to more accurately approach 
the true blockage. 

Fig. 15 shows a projection of the quadripod structure, 
as seen by a ray from the focus reflected by the inner 
portion of the reflector (inside the quadripod structure). 
Optically, the projection is 71% opaque. A similar projec- 
tion shows that for a ray coming directly from the focus 
towards the outer portion of the reflector, the structure 
is optically 47% opaque (approximately). After weighting 
by energy and area distributions, the effective average 
optical opacity is found to be approximately 60%. 

The exact relation between optical and RF opacity is 
not known. An estimate of RF opacity of 80% was taken 


n A. C. Ludwig and W. V, T. Rusch, Digital Computer Analysis and 
Design of a Subreflector of Complex Shape ( a JPL Technical Re- 
port to be published). 

4 The loss due to nonuniform amplitude illumination is fairly high 
because the feed system was designed for optimum figure-of-merit 
rather than optimum gain. Illumination efficiency was sacrificed to 
obtain extremely low spillover past the main reflector, to reduce 
antenna noise temperature (SPS 37-34, Vol. IV, pp. 239-241). 
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Fig. 14. Aperture blockage 



Fig. 15. Projection of quadripod structure onto aperture of the antenna 


as nominal. (The uncertainty in this estimate is discussed 
under error analysis, presented later in this article.) The 
subreflector was taken to be 100% opaque. On this basis, 
loss due to blockage, as shown in Fig. 13, is 0.62 db. 

3. Structural Deviations 

A structural analysis of this antenna has been consid- 
ered in detail previously."’ The resulting errors have an 
rms value of ^0.057 in., under very ideal conditions (no 

J M. S. Katow, private communication, 1966. 


wind or thermals, approximately 45-deg elevation angle). 
Although components of this error do not have a strictly 
normal distribution, direct aperture integration indicates 
close agreement with the value calculated with the 
method formulated by Ruze (Ref. 5). This value is a loss 
of 0.080 db. 

4. Error Analysis and Predicted Gain 

The ground rule adopted for this error analysis is to 
make estimates for each loss component that is believed 
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to have a 95% probability of bracketing the true value. 
For normal distributions, this would represent the 2 -<r 
deviation of these quantities. Total error is computed as 
the square root of the sum of the squares of individual 
errors. Although some errors are neither strictly normal 
nor totally independent, this is believed to be the most 
reasonable evaluation. 

Error estimates are shown in Fig. 16. By far the domi- 
nant uncertainty is due to blockage, a result of the uncer- 
tainty in the opacity factor as discussed earlier. The 
tolerance shown corresponds to opacity factors of 50 to 
100%. The total tolerance in loss factors is -f 0.17, -0.19 db. 

The sum of the loss factors in Fig. 13 is 1.85 db. Com- 
bining this with the gain for the ideal case, the resulting 
predicted gain is 61.90 +0.19, —0.17 db, under ideal en- 
vironmental conditions as mentioned earlier. Note that 
the tolerances listed apply only to the factors considered 
in the loss computations. 

Although every effort was made to make the analysis 
thorough, it is conceivable that a loss factor could have 
been omitted. Therefore, in a strict sense, the above tol- 
erance is not necessarily the overall tolerance on the pre- 


dicted value versus the true value. Ultimately, this 
uncertainty can be removed only by a good experimental 
gain measurement. 
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I. Efficient Antenna Systems: X-Band Gain 
Measurements, D. A. Bathker, and G. S. Levy 

1 . Summary 

In ideal antennas the gain increase is proportional to 
the frequency squared; however, in real antennas aper- 
ture phase errors, due to surface tolerance effects, do not 
permit the realization of this gain-frequency relationship. 
When the frequency rate of change of gain loss due to 
surface effects is equal and opposite to the f 1 law rate of 
change of gain, the antenna is said to be at gain limit. 
An experimental program has been undertaken to investi- 
gate the behavior of the DSIF antennas as they approach 
gain limit. Gain measurements have been made in the 
past at S-band, and the next step is measurements at 
X-band. 

An X-band 85-ft antenna cassegrainian cone has been 
built and will be used to measure gain at the Venus 
station. Celestial radio sources will then be observed, and 
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the flux density will be determined. The next step will 
be to measure the gain of the 210-ft advanced antenna 
system at X-band, using the same cassegrainian cone. Ulti- 
mately it is planned to repeat this procedure at C- and 
K-bands. 

2. Introduction 

In Ref. 6 the tradeoffs involved in optimizing large 
antennas for space communications are discussed. The 
concept of gain limit, which is the maximum gain a given 
antenna can achieve, is introduced. 

An ideal antenna s gain is proportional to the square 
of operating frequency. However, in real antennas the 
performance of an antenna is deteriorated by surface 
imperfections. The loss in gain is given by G/(G 0 ) = 
ex P — [(4 7r <r)/A] 2 , where A is the free space wavelength, 
and <t is the rms RF path length error due to surface 
deviations (Ref. 7). At some given frequency the deriva- 
tive of the gain loss due to surface imperfections cancels 
the f 2 gain increase; this is said to be gain limit. Fig. 1 of 
Ref. 6 shows gain versus frequency for the DSIF 85- and 
210-ft antennas for two ratios of rms surface deviation 
to paraboloid diameter. 



I 2 4 6 8 10 

FREQUENCY, Hz 


Fig. 17. Gain vs frequency, size, and surface 
tolerance for one ground antenna 

Measurements well below the gain limit have been 
made at S-band on the 85-ft dishes (Ref. 8), but there 
are very few experimental data in the region near and 
above gain limit. An experimental program has been 
started to measure the 85- and 210-ft antennas at X-band. 
There are plans to repeat these measurements at C- and 
K-band. 
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3. X-Band Measurement 

An X-band feed system has been designed to fit in a 
standard DSIF cassegrainian support cone. This feed sys- 
tem can be used on either the 85-ft Venus station antenna 
or on the 210-ft advanced antenna system. 

The first 8448 MHz tests will be conducted at the 
Venus station, where antenna gain and pattern measure- 
ments can be made using the Mount Tiefort collimation 
station. Measurements will also be made of the antenna 



Fig. 18. Diagram of feed and waveguide 
in the X-band cone 
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temperature before and after installation on the 85-ft 
paraboloid to determine the spillover contribution. Mea- 
surements of celestial source temperature will be made, 
and these data will then be used in conjunction with the 
gain data to obtain absolute source flux density. 

Fig. 18 is a diagram of the feed and waveguide system 
in the X-band cone. A dual mode horn, which is a scaled 
version of the Venus station radar feed (Ref. 7), is em- 
ployed. The feed will use circular polarization generated 
by the turnstile junction and fed through the mode gen- 
erator. A section of waveguide below the turnstile junc- 
tion is replaceable with two rotary vane attenuators to 
be used for antenna pattern calibrations. Waveguide 
switch 1 will permit either the maser and radiometric 
receiver or the crystal mixer receiver with automatic gain 
control (AGC) to look at either the antenna or the pre- 
cision terminations. The maser, radiometric receiver, and 
precision loads will be used for antenna temperature and 
radio source measurement and are discussed elsewhere in 
this report. The AGC receiver will be used for antenna 
pattern measurements and can be replaced by a thermistor 
and power bridge for gain measurements. 

The hyperboloid, flange, and vertex plate used in this 
experiment will be S-band devices. Ludwig (SPS 37-41, 
Vol. Ill, p. 90) has computed the perturbations caused 
by using the S-band devices at X-band, and these correc- 
tions will be used in determining gain. It was decided 
that construction of vertex plate and flange for X-band 
would not be desirable with regard to cost and station 
time. 

After completion of the Venus station experiment, it is 
planned to use this feed-cone at the Mars station. The 
same celestial sources will be examined, and by using 
the Venus station values of flux density, the gain can be 
measured for the 210-ft dish at X-band. 
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F. Low Noise Receivers: Microwave Maser 

Development, S. Petty, and R. C. Clauss 

1 . Summary 

A traveling wave maser (TWM) has been constructed 
for use at 8448 MHz. It provides a net gain of 41.5 db 
when installed in a closed-cycle helium refrigerator oper- 
ating at 4.4°K. The maser provides more than 30 db net 
gain over a tunable range of 140 MHz. An equivalent 
input noise temperature of approximately 18°K has been 
measured. 


2. Maser Structure Description 

The maser structure (Fig. 19) accepts both signal and 
pump waveguides at one end. Transitions from rectangu- 
lar to ridged waveguide are used to reduce the maser 
diameter. The slow-wave comb structure is folded, result- 
ing in 4 in. of active length in 2 in. of uniform magnetic 
field. The end loop, which couples the halves of the maser, 
is shown in Fig. 20. 

Loading of the comb structure is shown in Fig. 21. The 
maser material is ruby (0.05$ Cr 2 O a ). C-axis orientation 
of the ruby with respect to the DC magnetic field is 
54° 44', which allows push-pull pumping. Ten thin disks 
of polycrystalline yttrium iron garnet (YIG), spaced uni- 
formly throughout the comb structure, provide sufficient 
isolation to ensure stable amplification. Additional load- 
ing to give the desired center frequency and electrical 
length is provided by alumina strips. Beryllium copper 
shims 0.005-in. thick secure the alumina and ruby in place. 

In Fig. 22, the disassembled maser structure is shown. 
All parts are oxygen-free, high-conductivity copper, and 
pure indium sheeting is placed in all joints to prevent RF 
leakage and to improve heat transfer across these joints. 

3. Package Description 

The maser structure is shown installed in a closed-cycle 
helium refrigerator developed at JPL (SPS 37-39, Vol. Ill, 
p. 79) (Fig. 23). The signal waveguides are 0.010-in. wall 
stainless steel, copper plated with a gold flash. Commer- 
cial mica pressure windows provide the vacuum seal. 

Fig. 24 shows the complete maser package. The welded 
aluminum frame and the Alnico-5 magnet (modified to 
provide a 3800-gauss field for 8448-MHz operation) are 
the same as those used in existing S-band masers at the 
Mars and Venus stations. The pump package contains a 
Varian VA-282 klystron and provides 200 mw of pump 
power at 21.85 GHz. 

During laboratory testing, maser gain changes occurred 
when changing input terminations. These changes were 
found to be not a function of the match at the signal 
frequency (the maser being very insensitive to input or 
output VSWR changes), but a function of reflected pump 
energy which behaves in an unpredictable manner in the 
signal waveguides. Gain changes as large as 1 db oc- 
curred, even though the maser appeared well saturated 
with pump power. A low pass waveguide filter (15 GHz 
cutoff frequency) was installed at the maser input (Fig. 24). 
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INCHES 


Fig. 19. Maser structure partially disassembled 
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RUBY CLAMP 


^ \ — COMB STRUCTURE 

Fig. 20. End coupling loop in position behind 
comb structure 


REXOLITE 

SUPPORT 


RESONANT 
ELEMENT (FINGER) 


Fig. 21. End view of maser structure 
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Fig. 22. Maser structure disassembled 
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MASER 

STRUCTURE 



Fig. 23. Maser in closed-cycle refrigerator 

This filter has greater than 40-db rejection at the pump 
frequency and 0.042 db insertion loss at 8448 MHz. Gain 
changes due to input load changes are now less than can 
be measured with existing test equipment. 

Table 4. Maser performance at 8448 MHz 


Signal frequency 

— 
8448 MHz 

Net gain 

41.5 db 

Magnetic field 

3800 gauss 

Bandwidth (3 db) 

17 MHz 

Forward loss 

18.5 db 

Electronic gain 

60 db 

Measured inversion ratio 

2.8 

Pump frequency 

21.85 GHz 

Measured equivalent input 
noise temperature 

1 8°K 


4. Performance 

Table 4 summarizes the electrical characteristics of the 
TWM at 8448 MHz. The net gain of the maser (41.5 db) 
has a short term stability (minutes) better than 0.1 db, but 



Fig. 25. Bandpass characteristic of 
TWM at 8448 MHz 
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the long term stability is seriously impaired by tempera- 
ture variations affecting the magnet. The gain and band- 
width at 8448 MHz are shown in Fig. 25. The measured 
equivalent input noise temperature of 18°K is higher than 
the predicted value (based on insertion loss measure- 
ments of input waveguide components, and the calculated 
maser noise temperature). 

The maser can be tuned over a 140-MHz range with 
net gain exceeding 30 dh. This is shown in Fig. 26. Cool 
down time of the closed-cycle refrigerator with this maser 
is 4.6 hr (no liquid nitrogen precooling) and 2.4 hr with 
precooling. 



Fig. 26. Net gain as a function of frequency 




instrumentation, C T. Stelzried 


The X-band cone RF instrumentation is described. This 
cone is immediately adaptable to the Coldstone Venus 
and Mars stations antennas and will be used to provide 
operational performance data of these antennas at 
8448 MHz. A total power radiometer, using a maser am- 
plifier and precision waveguide terminations, will be used 
for the necessary noise temperature calibrations. 


The experimental X-band 8448 MHz cone can be used 
with either the JPL Coldstone 85-ft Venus station antenna 
or tiie 210-ft Mars station antenna. It will be used to 
provide information about antenna performance at this fre- 
quency. The first installation will be at the Venus station. 


Preliminary noise temperature calibrations of the re- 
ceiving system using the horn, maser, (see section III .F, 
this SPS) and monitor receiver have been made on the 
roof of Building 238 (Fig. 27). This equipment was used 
in a total power radiometer configuration with precision 
cryogenic waveguide' terminations. 


Results of the system temperature measurement will be 
presented at a future time. The difference between the 



Fig. 27. Photograph of 8448-MHz maser with liquid 
nitrogen cooled termination used in noise temper- 
ature calibrations on the roof of Bldg. 238 


system temperature measured with the horn connected 
directly to the maser input and that measured in the cone 
with the normal waveguide installation will provide a 
measure of the transmission line losses and a cross check 
on the normal waveguide insertion loss measurements. 

A block diagram of the complete instrumentation is 
shown in Fig. 28. Provisions are made for measuring the 
horn, ambient and cryogenic wave guide termination 
VSWIVs, maser gain, noise temperature and for monitor- 
ing the' receiver local oscillator frequency. Fig. 29 shows 
the waveguide system undergoing laboratory insertion loss 
and VSWR measurements. Fig. 30 is a photograph of 
the initial phase of the cone installation. Figs. 31, 32, 
and 33 show the monitor receiver, noise box, and rack 
instrumentations for the cone and control room. 
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Fig. 28. Block diagram of 8448-MHz cone and control room RF instrumentation 
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Fig. 30. X-band cone interior during initial phase of component installation 
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Fig. 31. Photograph of X-band monitor receiver 
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SWITCH N I 


SWITCH N2 


Fig. 33. Photograph of Venus station RF instru 
mentation rack X2 and the X-band cone 
instrumentation rack X3 
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Fig. 32. Photograph of X-band noise box 


RACK X3 
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H. Digital Communication and Tracking: 

Ranging Measurement, R. C. Tausworthe 

I. Ranging the Mariner Venus 1967 Spacecraft 

a. Introduction. In June 1966, NASAs Office of Space 
Science and Applications approved ranging as a part of 
the celestial mechanics experiment of the Mariner Venus 
67 mission. Inasmuch as plans had already been made to 
equip the spacecraft with a turnaround transponder, it 
remained only to provide a ground system capable of 
extending measurements out to planetary distances. 

The ranging measurements to be made can be taken 
in two parts: 

(1) The first part consists of ranging the first few million 
miles of the mission, during which time range mea- 
surements can be taken from any of the assigned 
tracking stations using standard DSIF equipment. 
This kind of ranging data is of value in improving 
orbit determination, and thus guidance, and in pro- 
viding actual system checks with the ranging tech- 
nique to be used later in the flight when the DSIF 
standard ranging system thresholds. 

(2) The second part of the ranging mission is to be 
conducted from Goldstone by R&D personnel and 
equipment to supplement the DSIF. In addition to 
improving the tracking accuracy during the latter 
part of the mission, the ranging data will be very 
valuable in two other ways: On the one hand, rang- 
ing will extend more precise and deep results to the 
celestial mechanics experiment, for while doppler 
and range are certainly related, the relationship 
depends on the details of the gravitational field — 
the area to be explored. Stating it in another way: 
doppler and range provide different ways of looking 
at the situation, and having both permits greater 
insight. At present, for instance, compatible range/ 
range-rate ephemerides of the planets are non- 
existent. 

The second way in which ranging data is helpful is as a 
probe of the intervening interplanetary medium. Ranging 
will provide additional information here because, while 
doppler measurements are related to the carrier, and thus 
the phase velocity of the signal, ranging is related to 
modulation, and thus to the group velocity of the signal. 
In a vacuum, the two are supposedly the same, and equal 
to the speed of light. In a dispersive medium, such as an 
ionosphere, the two are not the same. Therefore, having 
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both provides a sensitive way of measuring the dispersive 
effect, and compensating for it in translating time of flight 
or doppler into true range. 

The experiment will require the 210-ft Mars DSS 
antenna fitted with a 20-kw transmitter plus a supplement 
to the existing receiver. It is probable that at least 1 to 
1 ¥2 hr per pass at DSS 14 will be necessary during the pre- 
encounter period. Measurements immediately prior to 
encounter, and as near post-encounter emergence as 
feasible, are of extreme value in determining the astro- 
nomical unit. 

This report describes the ranging system in detail, and 
the ground system mechanization in particular. In sum- 
mary, however, there is a simple straightforward extension 
of present ranging techniques which will enable precise 
turnaround ranging of a spacecraft at planetary distances. 

b. Expected signal and noise parameters near encounter . 
Based on an 80 X 10° km communications distance along 
with the nominal transmitter-receiver parameters given in 
Tables 5 and 6, the spacecraft carrier tracking loop will 
be operating at about 21.9 db margin, giving about 11 or 
12 deg RF phase jitter, and the ground receiver carrier 
tracking loop will have about 25.7 db margin, or about 
7 or 8 deg of jitter. Also as for the ranging sidebands, the 
S/N into the transponder will be about — 19 db and the 
ground-received ranging power about — 172 dbm, giving 
a total-ranging-power-to-noise-density ratio S/N (1 equal to 
about 9.6 db. 

c. The ranging equipment . The R&D receiver ranging 
supplement consists of two correlation channels and a 
stored program signal processor, as depicted in Fig. 34. 
The correlation channels are constructed from modules 
much like, or the same as, modules in the GSDS receivers. 
In order to provide IF stability required for the narrow 
bandwidths involved, the usual free-running reference 
oscillator in the GSDS receiver is replaced by an atomically 
derived reference. All the reference frequencies indicated 
in Fig. 34 are coherently related to an atomic standard, 
the frequency and timing system (FTS). 

The signal processor consists of an SDS-920 general- 
purpose computer with analog/digital and digital/analog 
interfaces, supplemented by a number-controlled oscil- 
lator (NCO), various digital phases detector, counters, etc. 

The automatic feature of the programmed device makes 
it possible to effect code acquisition and range tracking in 

JPL SPACE PROGRAMS SUMMARY 37-42, VOL. Ill 



Table 5. Analysis of a turnaround ranging system 


Table 6. Carrier-tracking loop analysis 


Uplink 

Ground transmitter 

73.0 dbm 

20 kw 

Transmitting antenna gain 

59.5 db 

210-ft dish at 2115 MHz 

Ranging modulation 
suppression 

- 0.6 db 

87% of power in code 

Space loss 

~256.8 db 

80X 10 6 km at 21 15 MHz 

Receiving antenna circuit 
gain 

3.8 db 

Low gain at encounter 

Spacecraft received ranging 
power 

-121.1 dbm 


Receiver noise density 

-164.3 dbm/Hz 

2700° K noise 
temperature 

Transponder noise bandwidth 

62.2 db Hz 

3.3 Me -1.65 MHz 
baseband 

Total spacecraft receiver 
noise 

-102.1 dbm 


S/R into transponder 

- 19.0 db 


Downlink 

Spacecraft transmitter power 

40.0 dbm 

10 w 

Transmitter circuit loss 

— 1.8 db 


Transmitter antenna gain 

18.5 db 


Ranging modulation 
suppression 

- ll.Odb 


Limiter suppression 

— 21.0 db 

P(ranging)/P(limiter 

output) 

Space loss 

-258.0 db 

80X 10 6 km at 2298 MHz 

Receiving antenna circuit 
gain 

61 .3 db 

210-ft dish at 2298 MHz 

Ground received ranging 
power 

-172.0 dbm 


Receiver noise density 

-181.6 dbm 

50° K maser 

S/No of ranging signal 

9.6 db 


Effective tracking bandwidth 

- 17.0 db Hz 

0.02 Hz 

S/N in tracking bandwidth 

26.6 db 


Code tracking jitter after 
acquisition 

5.0 m 



the most effective way. For example, minute gain offsets 
in each channel can be compensated for by switching 
codes and channels back and forth periodically. 

There are seven steps to the acquisition process: (1) the 
clock is acquired and the clock-tracking loop locked, 
(2) through (6) each component is acquired, and (7) the 
entire code is tracked. 



Uplink 


Ground transmitter 

73.0 dbm 

20 kw 

Transmitter antenna gain 

59.5 db 

210-ft dish at2115MHz 

Carrier suppression 

- 8.9 db 

1 3% of power in carrier 

Space loss 

-256.8 db 

80 X 1 0 6 km 

Receiving antenna circuit 

3.8 db 

Spacecraft low gain at 

gain 


encounter 

Spacecraft received carrier 

-129.4 dbm 


power 

Receiver noise density 

-164.3 dbm/Hz 

2700° K noise 
temperature 

Spacecraft receiver design 

13.0 db Hz 

20 Hz 

bandwidth 

Receiver threshold 

-151.3 dbm 


Spacecraft carrier loop 

21.9 db 


margin 


Downlink 


Spacecraft transmitter power 

40.0 dbm 

10 w 

Transmitting circuit loss 

- 1.8 db 


Transmitter antenna gain 

1 8.5 db 


Carrier suppression 

— 5.1 db 


Space loss 

-258.0 db 

80X10° km at 2298 MHz 

Receiving antenna circuit 

61.3 db 

210-ft dish at 2298 MHz 

gain 

Ground received carrier 

— 145.1 dbm 


power 

Receiver noise density 

-181.6 dbm/Hz 

50° K maser 

Carrier loop design 

10.8 db Hz 

12 Hz 

bandwidth 

Ground receiver threshold 

-170.8 dbm 


Ground carrier loop margin 

+ 25.7 db 



There are two different tracking modes: the first is an 
early-late correlation-and-subtraction mode to produce an 
S-curve, as shown in Fig. 35, used in steps (1) and (7) only. 
The other mode is RF doppler-aided track during steps (2) 
through (6). Briefly, it operates as follows: 

Both the transmitted frequency and the clock frequency 
are accurately known, as both are developed from the 
atomic standard. With the carrier loop locked, the RF 
doppler detector provides an exact measurement of the 
incoming clock doppler (except for the slight effects of a 
dispersive medium), when multiplied by the clock-to- 
transmitter frequency ratio. It is thus possible, once the 
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Fig. 34. Manner 67 Venus ranging equipment 






































Fig. 35. Early-late correlate-and-subtract scheme 
to create tracking S-curve 


clock component is locked, to cause the receiver clock to 
vary exactly as it should by making it follow the RF 
doppler, properly scaled. Charged-particle effects and 
imprecise f c i/fr ratios cause a slight receiver clock drift 
which can be readjusted periodically as needed by revert- 
ing to tracking of the first type. 

However, one notes that steps (2) through (6) are accom- 
plished with both early and late channels free for cor- 
relation with the incoming code, speeding component 
acquisition by a factor of two over the one-correlation- 
channel Mark I system. 

The RF doppler can be used during steps (1) and (7) to 
reduce the steady-state phase error in the code tracking 


loop produced by a doppler acceleration, by merely pro- 
gramming the VCO to remove (measured) doppler effect. 

During steps (2) through (6), it is necessary to generate 
a waveform at the correct (computed) clock doppler offset, 
to which the measured clock doppler can be compared 
and thus be made to follow. Generation of the clock 
doppler offset frequency is conveniently done by a 
number-controlled oscillator (i.e., one whose frequency is 
set by a number put into it); but there is a practical limit 
on the lowest frequency it can generate. To eliminate this 
difficulty at zero and low doppler offsets, a clock doppler 
bias frequency (e in Fig. 34) is coherently inserted, both 
as an analog frequency into the clock doppler detector 
and as a known number into the NCO. 

d. Range code components . The transmitted ranging 
code is a six-component composite code of the same 
general type used in the Mark I system. The first of the 
components is a period-2 series of alternating 0’s and l’s, 
which we hereafter refer to as the clock component. The 
remaining five components are pseudonoise sequences of 
length 7, 11, 16, 19, and 23, generated by the recursions 
listed in Table 7. These are inserted into a majority- vote 
Boolean function in which tie votes are resolved to pro- 
duce a balance of l’s and Os. This balance is necessary 
to minimize the DC component of the transmitted code. 
A majority-vote logic is known to maximize the sum of 
correlations between the transmitted code and each com- 
ponent, and tie votes can then be resolved in such a way 
as to make these correlations be approximately propor- 
tional to the periods of the components, to minimize range 
acquisition time. 

Receiver component codes are duplicates of those trans- 
mitted, generated at a rate set by a VCO programmed to 
track the incoming code rate as described in part c. 

During steps (1) through (6) of the acquisition pro- 
cedure, adjacent phases of the component sequences 
themselves are inserted into the correlation channels; 
step (1) — clock acquisition — thus proceeds by early-late 
tracking as described. At the end of step (1), bit transi- 
tions of the incoming transmitted code and those of the 
clock are displaced 90 deg. During steps (2) through (6), 
the receiver clock VCO is programmed over 90 deg so 
that bit transitions are aligned to give the most effective 
detection. At step (7) the 90 deg shift is removed, and 
the system again is put into the early-late correlation 
tracking mode. 
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Table 7. Ranging code components 


e. Initial range acquisition. At the start of a run, both 
the transmitter coder and the receiver coder are run from 
the atomic-standard 1-Mc code rate, and all code com- 
ponents are synchronized together. Then the receiver 
coder is switched to its own VCO. From the time this 
switch occurs until clock acquisition (lock), the two sets 
of codes may separate in phase in an undetermined man- 
ner; however, by monitoring the clock doppler counter 
(and removing the bias) the exact amount of slip is 
recorded and thus accounted for. At the end of step (1) 
(clock acquisition), the receiver coder VCO frequency is 
the same as that of the incoming code; at the beginning 
of step (2), there is a fixed integral number of phase-steps 
separating events of the incoming code with those cor- 
responding ones in the receiver components. Step (2) is 
thus concerned with correlating the incoming stream 
with each of the seven possible phases of the second 
component. This could be accomplished serially by aver- 
aging at the first phase for a time T, recording the value, 
stepping to the second phase, and repeating the proce- 
dure until all seven correlations are obtained. That 
phase number m 2 having the largest correlation value is 
regarded as the most probably correct choice, so this 
component is shifted then m 2 steps to bring the second 
component into alignment. This could be repeated with 
all successive components to align all components. With 



ber Theory) listed. It is interesting to note that, because 
of the way the clock is acquired, is always zero in the 
computation for M. 


two channels available for correlation, however, the 
process can be speeded somewhat by obtaining two 
averages simultaneously. 


N § ZigitllSo9rli6 


cation and Tracking: Ranging 


Subsystem Acquisition Probabilities, R. Didday 


At the end of acquisition we have, by manipulating 
each component, moved the code an equivalent of M 
steps; we have stored five numbers At any 

subsequent time we have available in the computer the 
accumulated clock doppler count, which indicates the 
integral number of steps in addition to M, separating the 
transmitted and received codes, and we also have the 
fractional number of clock cycles as a vernier to the clock 
doppler counter. 

f. The computation of range. At any desired time after 
acquisition the integers = M (mod p*), the clock dop- 
pler count, and the fractional doppler cycle count (or 
vernier) are available to compute the range. The range 
in light microseconds is thus given by 

R = M + doppler count 4- vernier 

Table 7 shows how M is computed from each of the 
numbers m x : M is the linear combination of each m x with 
the coefficients (called “Chinese” numbers because they 
are a result of the Chinese Remainder Theorem of Num- 


1. Introduction 

The Mark I ranging subsystem now in use in the Deep 
Space Network has been previously described (S PS 37-20 
through 37-24 and 37-27, Vol. IV). Briefly, the system 
provides range measurements by detecting the phase shift 
between a modulating signal generated on the ground 
and the same signal after return from a vehicle being 
tracked. Detection of the phase of the received signal is 
accomplished by correlation techniques. The signals used 
are binary pseudonoise sequences chosen to have a long 
period (thus permitting unambiguous measurement of 
large distances) and yet to have characteristics which 
allow rapid acquisition (detection). 

It is, of course, essential to be able to test the system, 
and one technique for so doing is based on the behavior 
of the probability of error in acquisition as a function of 
the signal to noise ratio (S/N) and the integration time 
used in correlation. An error occurs when excessive noise 
causes a particular (wrong) phase position to appear to 
have the highest correlation of all possible positions, and 
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thus the occurrence of an error is related to the S/N. The 
test procedure consists of adjusting the S/N to a desired 
value, then attempting to acquire the signal from an 
object at a known range (a collimation tower, for ex- 
ample). A certain number of trials are made, and from 
the number of successes, the probability of correct acqui- 
sition on each trial can be placed within some range. 
If enough trials are performed, this range is small enough 
so that this technique provides an adequate measure of 
the system’s performance. 

2. Acquisition probabilities 

To use such a testing procedure, first of all we need to 
find the relationship between the probability of correct 
acquisition and the other system parameters. Using the 
result presented by Viterbi (Ref. 9, p. 120) we have the 
expression 


where C max and C m i n are the appropriate values of the 
correlation. Making the substitution x = /*,/( 2)^ we 
obtain 



X £l + erf ^ x + 

This expression may be evaluated given the value of 
ST /N Q and the characteristics of the code. 

In the Mark I subsystem there are two different codes: 
(1) the normal (or long) code and (2) the short code. Each 
of these is of the form 



P (correct acquisition) = P c = 



exp ( — p 2 / 2) 

(2rr) Va 





+ (2ST/No) 1/2 exp ( — o 2 /2) 


for an orthogonal, single component code. In this expres- 
sion, S is the signal energy, T is the integration time per 
phase position, N 0 is the noise power single-sided spectral 
density, and M is the period of the code. Since in the 
case we wish to analyze we do not have an orthogonal 
code, we must alter this expression. The nonorthogonality 
of the code we are using results in smaller changes in 
the correlation levels obtained, so the effective signal 
power is decreased by a factor of the square of the dif- 
ference in the maximum and the minimum correlation 
observed. Making this substitution, and using the defi- 
nition 



e ~ t2/ 2 dt 


the expression reduces to 


xcl + x [{ab + be -f ac) ® cl] , 

where x, a, b , and c are pseudonoise sequences with 
specified periods and cl clock is a component of 2-bit 
length. The clock component is acquired by locking the 
receiver loop, and does not enter into the correlation 
process. To discover the overall probability of correct 
acquisition it is necessary to compute the product of the 
probability of correct acquisition for each component. 
The effect of the different receiver codes generated by 
the different program states (Ref. 10) is taken into account 
by the values used for the maximum and minimum corre- 
lations. 

These calculations have been made for both the long 
and the short code, using an IBM 1620 computer. Curves 
of probability of error ( = 1 — P c ) versus S T/N 0 (where T 
is the integration time per code component phase posi- 
tion) are shown in Figs. 36 and 37. The results obtained 
are believed to be correct to four significant figures, and 
this limit is imposed only because rather large intervals 
were used in the numerical integration. Using more 
computer time, more precise answers could easily be 
obtained, since the rest of the computation is accurate to 
six significant figures. The values used for C max and C min 
were taken from Ref. 11, p. 19. 



By performing a suitable mapping, it is possible to 
obtain curves of the total acquisition time versus S/N 0 
(sideband power-to-noise-density ratio) with the proba- 
bility of error as a parameter. For example, from Fig. 36 
we note that for P (error) — 0.001, S T/N 0 — 374. Thus, 
a plot of T (time per phase position) versus the ranging 
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Fig. 36. Probability of error in acquisition for long code 



ST/N 0 — 


Fig. 37. Probability of error in acquisition for short code 


signal-to-noise-density ratio is the hyperbola (or straight 
line on a log-log plot) given by (S /N 0 ) — 374 /T. Total 
acquisition time for the normal code is calculated from 
T as 232 X T + 1.59 sec, since there are 232 (the sum 
of the component periods for the long code) phase posi- 
tions, and there is a 1.59-sec period required for the 
system to operate (Ref. 12, pp. 26-28). Plots obtained by 
this technique are shown for various probabilities of error 
for both the long code (Fig. 38) and the short code (Fig. 
39). Given the signal-to-noise ratio, from these plots it is 
possible to determine the time required for acquisition 
with a particular probability of error. 
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RANGING SIDEBAND POWER/NOISE SPECTRAL DENSITY, db— > 

Fig. 38. Long code acquisition time 



Fig. 39. Short code acquisition time 
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3. Estimation of P (correct acquisition) 

Next, we wish to be able to relate a certain number of 
successes in a given number of trials to a range of values 
which probably includes the actual probability of cor- 
rect acquisition on each trial. That is, we wish to find a 
confidence interval on P c , the probability of correct acqui- 
sition. To do this, we note that the number of successes 
in n trials is governed by the binomial distribution, that 
is, the probability of k successes in n trials given that the 
probability of success on each trial is P c , is 


Table 8. Acquisition confidence testing 


Number of 
successes observed 
in 10 trials 

95% confidence limits on 
Pfcorrect acquisition) 

Limits on Pterror 
in acquisition) 

6 

0.2650 

0.8786 

0.1214 

0.7350 

7 

0.3475 

0.9336 

0.06640 

0.6525 

8 

0.4439 

0.9752 

0.02480 

0.5561 

9 

0.5548 

0.9965 

0.00346 

0.4452 

10 

0.6914 

1.0 

0.0 

0.3086 


( r) (V a - r,)-* = (r,r (1 - r.)- ■ 

However, our problem is somewhat the inverse of this. 
What we wish to know is some estimate of the value of 
P c given that k successes were observed. Specifically, we 
wish to find some P x and P 2 so that 


P{P 1 <P C <P 2 ) =0.95. 


We can do this by forcing P± and P 2 to satisfy the re- 
quirements 


P {k < k 0 | P x ) = 0.025 
P{k>k 0 I P 2 } = 0.025 


where k 0 is the number of successes actually observed in 
the total number of trials (Ref. 13). 

We know the distribution of k given P 1? so we can 
write these two expressions explicitly as 


£(£)(?,)* a -Pi)*-* 
£(£)(P 2 )(i-p 2 ) n -* 


= 0.025 


= 0.025 . 


Solutions for various values of k 0 are accomplished by 
use of tables of the binomial distribution (Ref. 14). A 
table of the results for one case of interest (n = 10 and 
k 0 = 6, 7, 8, 9, 10) is shown in Table 8. 



SPECTRAL DENSITY, db — 

Fig. 40. Operating regions (95% confidence) for 
eight successes observed in ten trials 


The proper interpretation of this graph is that given a 
particular time switch setting, and given that 8 correct 
acquisitions were observed out of ten trials, then the 
probability is 0.95 that the system is operating within the 
interval (bar) shown. This interval represents about an 
order of magnitude change in the probability of error 
but corresponds to only 3 or 4 db in S/N 0 . Thus, the 
operating point is fairly well specified. 


These limits on P(error in acquisition) may be mapped 
onto the curves of Figs. 38 and 39 to produce a chart 
which is of use in testing the ranging system. A sample 
of such a plot, done for 8 successes is shown in Fig. 40. 


One method of using a curve such as this would be to 
select a time switch setting (thus fixing the total acqui- 
sition time), and then attempt to acquire an object at 
known range after causing S/N 0 to assume some desired 
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value. For instance, suppose we had selected time switch 
setting 10, had caused the signal-to-noise ratio to be 
30.8 db, and had obtained 8 successes. We would note 
from Fig. 40 that with 95% confidence, the operating 
point of the system is in the interval shown (which runs 
from S/N 0 — 28.83 to 32.87 db). We would conclude in 
this case that the system is functioning properly. If we 
had obtained some other number of successes, a plot 
similar to Fig. 40 would have been used to determine if 
the system seemed to be operating in the proper region. 

4. Remarks and Summary 

Briefly, a technique has been developed which allows 
quick computation of the probability of correct acquisi- 
tion in a ranging system. Information concerning the 
transmitter and receiver codes must be obtained in the 
form of length of the period of each component and 
the change in correlation before and after that code has 
been acquired. Using these probability functions, it is 
possible to determine if a particular system is operating 
properly by making a number of trials under known con- 
ditions of signal-to-noise ratio and range. Although the 
results presented here are specifically for the Mark I 
ranging subsystem, the programs used to obtain these 
results are general enough to allow similar results to be 
quickly obtained for future ranging codes. 

w6 7 1591 7 J , 

J. Digital Communication and Tracking: a Phase- 
Locked Receiver Analysis Program, 

R. C. Tausworthe 

1. Introduction 

A receiver analysis program to present a profile of a 
receivers phase noise and loop parameters under oper- 
ating conditions was reported in SPS 37-38, Vol. Ill, pp. 
29-33. In SPS 37-40, Vol. IV, pp. 219^225, a method was 
presented for the cascade analysis of Mark I ranging 
system clock loop jitter. The present article describes a 
program operable on an SDS-920 or SDS-930 for analysis 
or design of either or both of the carrier-tracking and 
clock-tracking loops of a DSIF receiver. 

2. Loop model 

The nonlinear theory of loop operation which this pro- 
gram embodies is developed in Refs. 8 and 9; the more 
interested reader is referred to these for detailed infor- 
mation. Basically, the program solves an equation 

= 2 Po m ( w* ) (~) (1) 


in which m is the receiver margin above loop design 
signal power level F 0 » p 0 is the SNR in the receiver pass- 
band w Lq Hz wide at this power level, p 0 = P 0 /kTw Lo , 
w Lo is the two-sided linear theoretic value of loop band- 
width at design point, w L(eq) is the actual equivalent loop 
bandwidth, T is a limiter performance factor, and y 2 is a 
parameter relating to nonlinear operation of the loop. 
Both w L(eq) and y 2 are functions of the parameter a 2 , so 
Eq. (1) is transcendental. The program uses Newtons 
method to find solutions for a 2 both in the carrier and 
clock loops. For the clock loop, y 2 , w L(eq ), etc. are func- 
tions of a\ arr which is thus determined first. 

Once the a 2 are found, they are related to the phase 
error variance by an approximate formula 

<T 2 = jl - exp (1 + 0.13a 2 ) . 

(2) ' 

The remainder of the loop performance parameters follow 
by mere substitution (Ref. 9). 

3. The Program 

The program was FORTRAN-compiled, using SDS 
series 900 FORTRAN II from which a paper tape dump 
was made, suitable for either the SDS 930 or SDS 920 
without further modification. It would be a simple matter 
to recompile and dump to provide an SDS 910 program. 
The tape is entered into the machine by a “standard fill.” 

Upon completion of the fill, control is automatically 
given to the program, causing the typewriter output 

TYPE OUT ANALYSIS — CARRIER LOOP, 
CLOCK LOOP, OR BOTH 

whereupon the operator responds CARRIER, CLOCK, or 
BOTH . The computer then asks for the computer con- 
figuration in carrier and clock loops, and for the type 
numerical data to be given it. The last can be either 
DESIGNED, MEASURED, or THEORETIC param- 
eters. 

4. Parameters 

Values called for by the analysis are: 

WLO = W Lq , loop linear theoretic 
two-sided bandwidth at 
design point, Hz 
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DAMPFAC(LIN) = £ 0 , the loop linear theoretic 

damping factor at design point 

T2 = numerator time constant of loop 
filter, sec 

= r 2 

T 1 = denominator time constant of 
loop filter, sec 

— T i 

BH = predetection one-sided 
bandwidth, Hz 
= w H / 2 

MARGIN = starting value of analysis, db 
= m 

A MARGIN = difference in margin between 
steps, db 

TMP = noise temperature, °K 

P/PTOT ~ ratio of power in carrier (clock) 
channel to total power, db 

PO/KTWLO ~ p 0 , the design point SNR, db 

LOOP GAIN = /* „ j! I / r“’, S ‘T le , l00 ? ) 

\K d K VC oMF (1/sec, limiter loop) 

CORR MAX = value of clock-loop correlation 
when a particular component 
is acquired (absolute) 

CORR MN = value of clock-loop correlation 
when that particular component 
is unacquired (absolute) 

5. Output 

A typical computer run is shown in Fig. 1. Most of 
the quantities are self explanatory. The phase error is 
printed in degrees, radians and radians-squared, both 
actual value, and linear-theoretic value. The limiter per- 
formance factor ( LIM.PER.FAC ), or T, is shown, and 
the limiter suppression factor a is also given. 

6. Cautions 

To eliminate as much operator failure as possible, sev- 
eral error checks were put into the program. Some things 
of note are: 

(1) False input characters recognized: 

[BOTF] 

NOT ACCEPTABLE RESPONSE , MESSAGE 
IGNORED 


(2) To restart, put machine in IDLE , press START , 
then STEP-IDLE-RUN . Also can be done by man- 
ual branch to cell 03507. 

(3) Only the first four letters of assignment words are 
recognized, although nominally these are given as: 
INCLUDED , EXCLUDED , SIGNAL , NOISE , 
DESIGN , MEASURED , THEORETIC , CARRIER 
LOOP , CLOCK LOOP , BOTH . 

(4) When analyzing the carrier loop only, at the end 
ANALYSIS COMPLETE is typed. If BKPT 1 is 
down , the entire program is reinitiated. If BKPT 1 
is up, the program only asks for a new set of pa- 
rameters. 

(5) When analyzing the clock loop, the computer asks 
for a new CORR MAX and CORR MN, to com- 
plete the ranging profile. Only the clock loop per- 
formance is printed subsequently. 

N67 15918 

K. Digital Communication and TrackingrTime 
Synchronization Experiment, W. Martin 

1. Introduction 

The following article discusses the results of time syn- 
chronization experiments conducted with several DSN 
stations. Over a two-day period the Mark I ranging sys- 
tem accurately measured time offsets between paired 
stations. The resulting clock errors were later verified 
when the Naval Observatory independently calibrated 
the stations’ clocks. 

2. Experiment 

Orbit calculations for the Lunar Orbiter spacecraft 
placed stringent requirements on the timing accuracy of 
the DSN. Whereas previously it had been adequate to 
synchronize two stations to within 10 msec, now it became 
necessary to calibrate the station clock offset to better 
than 50 fx sec. 

A new synchronization method was devised for which 
a tabulation of possible error sources showed the worst- 
case uncertainties to be substantially below 10 p sec 
(Ref. 17). The new idea utilized the existing Mark I 
ranging system and a spacecraft, thereby realizing sub- 
stantial cost savings over competing methods. 

Time synchronization tests were conducted on August 
16 and 17, 1966. Three stations participated in the experi- 
ment: DSS 12, DSS 41, and DSS 61. Each test consisted 
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SECOND-ORDER PHASE-LOCKED RECEIVER ANALYSIS 


CARRIER LOOP WITH LIMI TER, S I GNAL VARIATION 

CLOCK LOOP WITHOUT LIMITER, NOISE VARIATION 

INITIAL PARAMETER VALUES 


CARRIER LOOP- 



T1,SEC 

5.351E 

02 

T2.SEC 

1.249E- 

•01 

LOOP GAIN - 

1.000E 

06 

ALPHA0 OR A0- 

6.852E- 

■02 

TEMP, DEG K - 

5.000E 

01 

CLOCK LOOP— 



T1,SEC 

1.482E 

02 

T2.SEC 

1.482E 

00 

LOOP GAIN - 

5.084E 

12 

ALPHA0 OR A0- 

2.627E- 

-11 

MX CORR VAL - 

1.000E 

00 


WL0,CPS- 1.200E 01 
BH.CPS - 2.000E 03 
R0 * 1.999E 00 
T2/T1-Z- 2 .335E-04 


WL0 ,CPS* 1.000E 00 
BH,CPS - 2.000E 03 
R0 - 1.979E 00 
T2/T1-Z- 1.000E-02 
MN CORR- .000E 00 


KTWL0.DBM —170.821 
P0/KTWL0,DB- .000 
IF SNR,DB * -22.219 
P/PTOT ,DB - -3.000 

KT, DBM/CPS —181.612 


KTWL0,DBM —181.612 
P0/KTWL0,DB- .000 
IF SNR.DB - -33.010 
P/PTOT, DB - -3.000 


RECEIVER PERFORMANCE 


CARRIER LOOP — 10.000 


SIGNAL LEV.OBM— 160.821 
LOOP MARGI N,DB- 10.000 
P/KTWL[ACT],DB- 6.437 
P/KTWL[LIN],C9= 6.166 
IF SNR, DB - -12.219 


WLt L I N] ,CPS ■ 2 .90 IE 01 
WL T ACT ] ,CPS - 2 .726E 01 
CAMPINGtLIN]- 1.250E 00 
DAMPI NG r ACT ]- 1.206E 00 
LIM.PER.FAC.- 1.129E 00 


ERR[ACT1,DEG* 2.201E 01 
ERR[ ACT ] ,RAD- 3.841E-01 
ERR [ ACT ] ,RSQ* 1.475E-01 
ERRtUN],RSQ- 1.209E-01 
ALPHA OR A - 2.142E-01 


CLOCK LOOP — 10.000 


SIGNAL LEV, DBM— 181.612 
LOOP MARGI N,DB- 19.206 
P/KTWL[ACT],DB- 19.324 
P/KTWLtLIN),DB- 19.336 
IF SNR, DB • -13.805 


WLtLI N],CPS • 1.000E 00 
WLt ACT], CPS » 9.695E-0 1 
DAMPI NGE LIN ]- 7.070E-01 
DAMP I NGT ACT 1- 6.867E-01 
LIM.PER.FAC.- 1.000E 00 


ERRt ACT 1 ,DEG- 4.604E 00 
ERR[ ACT ] ,RAD” 8.035E-02 
ERRE ACT ] ,RSQ- 6.456E-03 
ERRt L I N ] ,RSQ- 5.825E-03 
ALPHA OR A - 2.627E-11 


CARRIER LOOP— 15.000 


SIGNAL LEV, DBM— 155.821 
LOOP MARGI N,DB- 15.000 
P/KTWL[ACT],DB- 9.186 
P/KTWLtUN],OB* 9.051 
IF SNR, DB - -7.219 


WL C L I N ] , CPS - 4.722E 01 
WLt ACT ], CPS ■ 4.577E 01 
DAMP I NGt L I NT- 1.643E 00 
DAMP I NGT ACT ]- 1.615E 00 
LIM.PER.FAC.- 1.070E 00 


ERRt ACT ],OEG- 1.505E 01 
ERR [ ACT ] , RAD “ 2.626E-01 
ERRt ACT ] ,RSQ- 6.897E-02 
ERR[LIN1,RSQ- 6.222E-02 
ALPHA OR A - 3.702E-0 1 


Fig. 41. Phase-locked receiver analysis program output 
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of measuring the time offset between two of these sta- 
tions. Three tests provided experimenters 5 with sufficient 
data to compute an “error of closure ” providing a valid- 
ity measurement on the system. 

The paired stations made two sets of ten measure- 
ments, one set with each station transmitting. A uniform 
set, based on consistency, was culled from each original 
set of twenty. Corresponding readings from each station 
were examined to determine the measured time differ- 
ence between the two stations and the spacecraft. This 
result was compared with the expected time difference 
as computed by the Orbit Determination Program. Re- 
sults are found in Figs. 42, 43, and 44 and represent the 
time offsets between the respective station clocks. 

Improper setting of the counter trigger level at DSS 
61 during the tests caused a small time measurement 


5 W. Martin and F. Bomcamp. 


bias. DSS 61 personnel set their counter to start on the 
trailing-edge of the 1-sec time tick rather than on the 
desired leading edge. The problem is illustrated in Fig. 
45. During the first half of the DSS 61 and DSS 41 tests 
the threshold level control was set on a relatively steep 
part of the trailing edge. A change in this setting prior 
to the second half caused the counter to be started earlier, 
where the change of voltage with time was not so great. 
The effect can be seen as a larger variation in the mea- 
sured time during the latter set of data (Fig. 42). This is 
understandable when one remembers that the counter is 
started when a specific voltage is reached. Thus, a lower 
rate of change of voltage with time affords more oppor- 
tunity for noise to influence the measurement. 

The effect is further illustrated in Fig. 43 (the DSS 41 
and DSS 61 tests). Here the counter's threshold was set 
on a relatively steep part of the timing pulse. The result 
was highly consistent data whose RMS variation was less 
than 0.3 /xsec. 
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Fig. 42. Time synchronization Stations 4T and 61 # August 16, 1966 
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Fig. 43. Time synchronization Stations 12 and 61, August 16, 1966 
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Fig. 44. Time synchronization Stations 12 and 41, August 17, 1966 


Figs. 44 and 45 represent more typical sets of data. 
From these it is obvious that the time synchronization 
scheme can produce results which are consistent to better 
than 0.5 pse c. Possible instrumentation biases notwith- 
standing, the data suggests that the technique is very 
accurate. 

Clock offsets together with RMS variations are shown 
in Table 9 for the three experiments. Due to frequency 


Table 9. Ranging time synchronization test results 


Day 

Time, 

GMT 

Stations 

Time offset, 
/xsec 

RMS 

variation, /xsec 

228 

07:00 

41 and 61 

9731.6“ 

1.27 

228 

18:00 

61 and 1 2 

-7637.5 

0.28 

229 

01:00 

1 2 and 41 

-2091.7 

0.39 

a Correction bi 

ased on first (trailing sdgs) meaiursmsnt. 
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STATION 61 


STATION 41 



SECOND TEST 


FIRST TEST 




TIME-*/ 


(a) STATIONS 61 AND 41 l-sec TIME TICKS 



ALL TESTS 


STATION 61 STATION 41 



(b) STATIONS 61 AND 12 l-sec TIME TICKS 

Fig. 45. Counter starting points 


drifts in the rubidium standards, these offsets are valid 
only at the stated times and cannot be used to compute 
the “error of closure” (a matter to be considered later). 
In the cases where DSS 61 participated, the time differ- 
ences have been corrected to the value that would have 
been obtained had the counter been triggered properly. 
The correction is based upon pulsewidth information 
supplied by the station. 

Considering the variation in the data (Table 9) the 
ranging time synchronization scheme should produce 
resolutions well below 1 fisec. This level should be at- 
tainable when the procedure becomes more familiar to 
the operators. Thus, the only uncertainties in the com- 
puted offsets should be due to differences in the ground 
equipments signal processing time. Even the most pessi- 
mistic estimates of the delay time difference place this 
figure well below a microsecond. 

One can test the validity of the synchronization scheme 
by summing the measured time offsets between stations 
throughout the DSN. A correction must be applied to the 
measurements in Table 9 to compensate for frequency 
drifts in the stations' master oscillators. These corrections 
are equivalent to adjusting all synchronization measure- 
ments to a common time (18:00 GMT on Day 228). Data 
relevant to oscillator drift was supplied by the Naval 
Observatory in an independent, but concurrent, time 
synchronization experiment. 


Table 10 shows the required corrections together with 
the adjusted offsets. Algebraically summing these time off- 
sets yields an “error of closure” of 0.5 /xsec. One s enthu- 
siasm for the result must be tempered by a knowledge of 
the uncertainty in the DSS 61 data. However, it seems 
probable that the offsets are accurate to within a few 
microseconds. 


Table 10. Adjusted ranging time synchronization 
test results 


Day 

Time, 

GMT 

Stations 

Correction, 

Time offset, 

/usee 

/itec 

228 

18:00 

41 and 61 

-1.6 

9730.0 

228 

18:00 

61 and 12 

0 

-7637.5 

228 

18:00 

12 and 41 

+ 1.3 

-2093.0 



Error of closure 

-0.5 


A further check on the ranging time synchronization 
experiment was provided by an independent station 
clock calibration. The Naval Observatory, acting under 
contract to the Lunar Orbiter Project, circulated a cesium 
clock throughout the DSN immediately following the 
ranging tests. Both time and frequency errors were mea- 
sured relative to the primary cesium standard. Frequency 
measurements were made over a period of several hours 
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Table 11. Station clock calibration by 
Naval Observatory 


Day 

Time, 

GMT 

Station 

Time 

synchronization 

error, 

/zsec 

Frequency 

error 

Drift 

rate, 

/i sec/day 

229 

12:00 

12 

157.4 lags 

1.74 XI 0' 11 fast 

1.50 

232 

19:40 

61 

7482.9 leads 

3.02 X 1 0~ u fast 

2.61 

237 

09:00 

41 

2189.7 lags 

7.00 X 1 0 -u fast 

6.04 


to eliminate short-term variations in the cesium clock. 
Test results are summarized in Table 11. 

To facilitate a comparison with the ranging data, the 
Naval Observatory measurements have also been ad- 
justed to 18:00 GMT on Day 228 (Table 12). The right- 
most column in Table 12 indicates the discrepancy 
between the Naval Observatory and the ranging time 
synchronization calibrations. In each case the ranging 
method produced larger offsets. There appear to be three 
reasons for the differences: (1) the basic uncertainties as 
outlined in Ref. 17; (2) the irregularity in the DSS 61 
measurements; (3) the assumption that the rubidium 
clock’s drift rates are constant, in adjusting both sets of 
data to a common time. This is unlikely, particularly in 
the case of DSS 41, where the drift rate is large and the 
adjustment interval long (almost nine days). Past expe- 
rience with rubidium vapor oscillators together with the 
literature (Ref. 18) indicates that the frequency stability 
of these units is on the order of 1 part in 10 11 . Thus, 
both the time separating the calibrations and the fre- 
quency variation could account for a large part of the 
discrepancy. 

3. Data Processing 

A computer program has been written for reducing the 
time synchronization data. This program is capable of 
accepting inputs from up to three stations simultaneously. 
In addition to computing the offset between paired sta- 
tions, the program also corrects for inaccuracies in the 
counter time base oscillator and for the nonuniform 
arrival of the code marker pulses. 

A simple calibration procedure was devised to measure 
the frequency error in the counters’ time base oscillator. 
The internal standard is counted between two 1-sec time 
ticks generated by the station’s frequency and timing 
system. Were the frequency correct, the resulting num- 
ber would be precisely 10,000,000. To the extent that it 


differs from that amount the error is known. Present plans 
call for this measurement to be made at each station just 
prior to a time-synchronization experiment. The counter 
calibration number is processed along with the data and 
is used to individually correct the time-interval measure- 
ments from each station. 


Range predicts are available at 1-min intervals from 
the Orbit Determination Program. Since these are valid 
only on the minute, a range correction must be made for 
the arrival time of the code-marker pulse. In some cases 
this delay can be as much as 2.5 sec. Failure to make an 
adjustment can cause serious errors in the predicted time 
difference between stations. The error’s magnitude will 
be determined by the code-marker pulse arrival time and 
the station-spacecraft geometry. 

Lagrange’s interpolation formula is used to compute 
the correct range (Ref. 19). 

R ( t ) = Rj + r (ARj) -!- r(f ~ 1} - (A 2 Ry) 



where 

r = T g enerall Y ran g e (1 < r < 2) 

and 

u = time interval from first predict to arrival of 
code marker pulse 

k — predict time interval 

R ( t ) = interpolated range at time t 

Rj — range corresponding to earlier of four predict 
times 

A Rj = first difference between earliest pair of predicts 
A 2 R } = second difference between first three predicts 
A 3 R ; = third difference between all four predicts 


Four range predicts are required for each interpola- 
tion. To minimize errors, the program selects a pair on 
either side of the desired time. These adjusted ranges 
form the basis for computing the difference in signal 
arrival times between the various stations. 
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Table 12. Adjusted Naval Observatory time 
synchronization results 


Day 

Time, 

Gmt 

Stations 

Time offset, 
/usee 

Difference 
ranging — Naval 
Observatory, 
fisec 

228 

18:00 

41 and 61 

9714.1 

15.9 

228 

18:00 

61 and 12 

-7630.7 

6.8 

228 

18:00 

12 and 41 

-2083.4 

9.6 


The program is presently being run on an SDS 930 
computer. Required inputs are the predicted station- 
spacecraft range, the time interval measurements, and 
the counter oscillator correction number. Entry is made 
via punched cards each of which includes the date, sta- 
tion identification number, and time. In addition to the 
clock synchronization error the printout includes the 
data, station ID numbers, time, corrected range predict 
(both range and time), and the measured time interval. 


2. Spectral Density Function Measurement 

The basic equipment used in measuring the spectral 
density function of noise n is shown in Fig. 46. The out- 
put of our equipment is the time average of n \ ; and this 
is equal to the area under the curve of S nn , the noise 
spectral density, versus f from — B f to +B/. By measur- 
ing the output of the equipment in Fig. 1 for low-pass 
filters with several different values of noise bandwidth 
B f , the curve of S nn versus f may be found. 





1 






-t 

if 

£ 

if 


Fig. 46. Equipment for measuring the spectral 
density function of a noise source 


4. Summary 

On the basis of experimental evidence, there seems 
little doubt that the ranging time synchronization scheme 
is both practical and accurate. The data confirms an im- 
provement in accuracy of approximately 1000 times over 
previous methods. Uncertainties of less than 10 ^sec 
appear to be easily obtainable, and the character of the 
data suggests that these may be reduced further to the 
1-fisec range. As station personnel become familiar with 
the procedure, the 1-jxsec accuracy should become a 
reality. 


N67 15919 


L. Frequency Generation and Control: VCO Noise 
Measurement and Optimum Tracking Filters, 

C. Cooper 

1. Introduction 


3. VCO Phase Noise Spectral Density 

Fig. 47 shows the linearized model of a phase-locked 
loop. The symbols appearing here are defined as 

0 = input phase process 

A 

0 = output phase process 

n = input noise process heterodyned to baseband 

n v = VCO phase noise process, referred to the input 
of a unity gain VCO 

A = RMS input signal voltage 

K = open-loop gain 

F(s) = tracking filter transfer function 

L(s) = closed-loop transfer function 

The effect of n and n v upon 0 may be expressed as 


This article describes a technique for measuring the 
phase noise spectral density function of a VCO and pre- 
sents some experimental results obtained through the 
use of this technique. The accuracy obtained with these 
measurements seems acceptable. Using the measured 
spectral density function of the VCO noise, the optimum 
tracking filter is determined, using the method of R. C. 
Tausworthe in Ref. 20. Experimental results are then 
presented which confirm that the tracking filter so de- 
termined is actually the optimum. 


0(t) = L(p) + - - p L(p) n.(t) (1) 



Fig. 47. Linearized model of a phase-locked loop 
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Let us now define n 0 (t) as the VCO phase noise process 
referred to the VCO output. For frequencies above the 
loop bandwidth, n 0 (t) ~ nv(t)/p. Setting the input signal 
0(t) equal to zero, the phase detector output is 

m = [1 - L(p )] - [1 - L(p)] n 0 (t) (2) 

Let us now consider L(/<o), the Fourier Transform of L(p). 
For frequencies well above the closed loop bandwidth, 
L(ja>) ~ 0, thereby causing <f> to be approximately 


mental formula for S« ono (/<o) ~ S <p<p(ju) has been found, 
we can deduce what S nv n v (/o>) should be. Set 6 and n/A 
equal to zero in Fig. 47 and observe that 

1_ 

n 0 (t) = ^ nv{t) (6) 

1 + — AK F(p) 

P 

Taking the Laplace Transform of the expression in the 
differential operator p, 


4>(t) = ~ n„(t ) (3) 

Assuming that n{t)/A and n v (t) are uncorrelated, the 
spectral density function of <f> is 

S^(/o>) = Sn n (/<o) + Sn o n o (/o>) (4) 

a a 

We will assume that 


1 

s 1 

1 + 1/s AK F(s) ~~ S + AKF(s) 

Examine the low-frequency asymptote of this transformed 
expression. 

lim 1 _ 1 

5^0 ^ + AKF(s) " AK 


S n 0 S no (/6)) > > S» n (j w ) 

a a 

S^(/w) — Sn 0 n 0 (j<») 


( 5 ) 


We will later show some evidence which justifies this 
assumption for the particular experimental configuration 
upon which tests were made. 


SM can be measured by the technique discussed 
previously. Let us suppose that this has been done for 
frequencies above the loop bandwidth. Once an experi- 


The low-frequency limit holds for frequencies appreci- 
ably less than the closed loop bandwidth. Therefore 
rto(t) = nv(t)/AK for frequencies of the range in which 
VCO noise has its greatest effect upon the phase-locked 
loop's performance. The spectral densities are related as 

( ? ) 

Usually, some of the effective gain AK appears as a phase 
detector constant and some as a VCO constant. In this 
case, AK in the above spectral density equation is re- 
placed by the VCO gain. 



Fig. 48. Experimental apparatus for measuring the spectral density function of the VCO in the 35.2-MHz loop 

of the central frequency synthesizer 
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4. Results 

The equipment used in making the experimental mea- 
surements is shown in Fig. 48. The spectral density 
function of the phase detector output <j>(t) will now be 
measured. The Hewlett Packard (HP) frequency synthe- 
sizer supplies a pure tone to the phase-locked loop which 
is equal to the VCO center frequency. Two Fairchild 
AOO-6 low-noise amplifiers are used as shown to mini- 
mize loading on the outputs of the phase detector and 
the low pass filter. They also provide a combined gain 
of 100 in order to make the signal of interest more easily 
observable. The low-pass filter consists of a passive RC 
network. No use was made of sharp cutoff active filters, 
because they introduce an excessive amount of noise. 
The filtered noise is read on a Griebach voltmeter, a 
passive instrument which indicates the true RMS value 
from DC up to 100 kHz. 

The Griebach possesses a damping time constant of 
several seconds, which is insufficient to average the very- 
low-frequency components present in the phase detector 
output. In order to obtain the average value of the square 
of the filtered noise, about 100 readings of the Griebach 
were taken at fixed time intervals. These readings were 
first squared and then averaged to obtain the mean square 
value of the filtered noise which appears on the left-hand 
side of Eq. (1). A computer program was written to 
perform this computation, and also to determine the 
standard deviation of the mean square value so obtained. 


Fig. 49 shows some recordings of the phase detector 
output before and after filtering and squaring. The 
squared signal was obtained from a square law detector 
in the Griebach meter. 

Table 13 shows two sets of data obtained in this man- 
ner. The second data set was obtained by bucking out 
the DC offset voltage before measuring with the Griebach 
and by taking a larger number of readings for each filter. 

The data in this table was analyzed in the manner 
indicated by Eq. (1). The spectral density function should 
be of the form (Ref. 20) 

S H (j2nf) = + ( 8 ) 


Table 13. Measurements of noise voltage 



Noise band- 

Mean of 

Standard devia- 

m«n» 

width of filter, 

squares. 

tion of mean of 

technique 

B f, Hz 

volt 2 


Normal 

0.191 

1 13.8991 

3.4898 

data 

1.844 

116.5746 

3.5718 

set 

20.34 

116.4412 

4.9303 


209.31 

115.2425 

3.6107 

Data 

0.407 

1.6746 

1.0058 

using 

0.787 

1.8763 

.7996 

bucking 

1.966 

2.4063 

1.1470 

voltage 

19.46 

2.5634 

1.8008 



Fig. 49. Recordings of the phase detector output <f> and <f>\ , the filtered and squared phase detector output, 

for various low-pass filter noise bandwidths 


JPL SPACE PROGRAMS SUMMARY 37-42, VOL. Hi 


69 



Analysis of the above data indicates that no N 0 term is 
present to within the accuracy limits of this experiment. 
Both sets of data obtained seem to indicate that S nn (j2*rf) 
may be somewhat more steeply sloped than 1/f , perhaps 
more like 1/f 2 . This can be partially explained, however, 
by the fact that the rather low cutoff rate of the low-pass 
filters will accentuate the rate of decrease of a decreasing 
function when this measuring technique is used. 

Table 14 shows the derivation of the best fit values 
of iV,. The average N, for the first data set is 

2V, - 40.7 X 1Q-* V“ 


*, 

we find that 


S 


Vo 


m) = 


30.5 X 10- 4 deg 2 

I/I 


Integrating S Vo (/27r/) from 1 Hz to 100 Hz (and doubling 
to get the single-sided equivalent) gives 0.028 deg 2 . 
W. F. Gillmore (Ref. 22) presents a phase noise spectral 
density* curve for the HP frequency synthesizer against 
its driver. Integrating graphically from 1 Hz to 100 Hz 
gives 0.003 deg 2 . Hence Snn(/2rr/) < < S no n o (j2arf) as 
assumed. aa 


o- = 144.0 X 10** v J 


The average N t for the second data set is 
N, = 50.2 X 10* v 2 
a = 25.4 X 10- G v 2 


Finally, Eqs. (5) and (7) yield 

c __ 0.473(Hz) 2 

*> n v n v I i I 


There is fairly good agreement between the average val- 
ues of N t for each set; but the variances show that the 
second set is much more reliable. If an appreciable N 0 
were present in N x would increase with increas- 

ing B f . Examination of Table 14 shows that this is defi- 
nitely not the case. On the basis of this data, it is safe 
to say that 


5. Verification of Optimum Filter Selection 

Let the tracking filter transfer function be 


F(s) = 


1 + t 2 S 
1 + Ti S 


N fl < 5 X 10"* v 2 


The assumption made in Eq. (5) can now be verified. 
Starting with 



Sw(j2irf) - 


50.2 X 10- ,! v 2 

— m — 


Fig. 50. Experimental apparatus for measuring 
phase noise variance 


Table 14. Derivation of best fit values of N x 


Measurement 

technique 

B Ft 

Hz 

mean squares 
A 2 
unsealed, 
volt 2 

mean squares 
A 2 

scaled, 

voU 2 

bf 2 

,n "B7 

N„ 

volt 2 

Normal data set 


1.3377 


2.266 

236.0 X 10-’ 


-0.0667 


2.400 

- 1 1 .0 X 1 0‘" 


209.31 

-0.5993 


2.331 

-102.8 X 10'* 

Data using bucking 

0.407 

0.787 

1.009 

43.6 X 10~ 8 

0.657 

63.3 X 10~ 8 

voltage 

1.966 

0.2650 

114.5 x ur* 

1.648 

69.4 X tO’ 8 


19.46 

0.0785 

33.9 X 10-" 

2.291 

14.7 X 10" 8 
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0.05 sec 0.05 sec 

Fig. 51. Phase noise for (a) present tracking filter (r — 2, B h = 0.5 Hz), 
(b) optimum tracking filter (r = 7, B L ^ 2 Hz) 



Fig. 52. Phase noise variance for tracking filters 
with various values of r and B L 


phase detector output is recorded and 6 o- is taken to be 
that length which includes 99% of the peak to peak 
jitter. This method is used to find a 2 for a number of 
tracking filters with different values of r and B L . Fig. 51 
shows recordings made in this manner for both the opti- 
mum tracking filter and for the tracking filter presently 
in the loop. Fig. 52 shows o- 2 for tracking filters with 
various values of r and B L . A definite minimum of o- 2 
occurs at r = 7, B L = 2 Hz, which is close to the pre- 
dicted optimum. The optimum filter yields approximately 
a 10%, or 1 db, improvement over the filter presently 
used in the loop, which has r = 2, B L = 0.5 Hz. 


M. Frequency Generi 





nd Control: S- and 


X-Band Central Frequency Synthesizer, 

G. U. Barbani 


The closed loop bandwidth B L and a damping factor r 
can be defined as 

D T + 1 AKr 2 

B l = — r = — 

4t 2 n 

Tausworthe has shown that in order to minimize that 
part of the output phase variance due to noise, the opti- 
mum value of r = 7 and the optimum value of B L is the 
solution of 

Bl - 0.143 Bl ” 0.774 ^ 'f-j = 0 (9) 

where k is Boltzmann’s constant. For the 35.2-Mc loop, 
A — 1.17 v mean square and T = 290°K, using the mea- 
sured values Ni v = 0.473(Hz) 2 and N ov = 0 yields 

B l = 2.7 cps 

The equipment in Fig. 50 is used to verify that the opti- 
mum values of r and B h are actually as predicted. The 


1. Introduction 

The central frequency synthesizer (CFS) has been in 
operation at DSS 13 since May 3, 1966. In SPS 37-40, 
Vol. III, pp. 57-61 a description is given of the addition 
of a X4 frequency multiplier (1- to 4-MHz) and a 4-MHz 
distribution amplifier. The CFS block diagram shown in 
Fig. 53 has been updated to include these modifications. 
A general report on CFS performance and the results of 
an investigation of the 35.075-MHz VCO static phase 
error are given in this report, 

2. Static Phase Error 

The 35.075-MHz phase-lock loop experienced a static 
phase error > 0.5 v for a period of about 4 wk. This error 
actuated the yellow (warning) condition lamp in the 
voltage vigilant system (VVS) which monitors all four 
phase-lock loops simultaneously (SPS 37-39, Vol. Ill, 
p. 66). Investigation of loop parameters such as phase 
detector drift, tracking filter components, in-lock detec- 
tor, VCO frequency adjust and bias, etc., all proved 
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Fig. 53. Central frequency synthesizer Mod I block diagram 


satisfactory. A reverse leakage current was discovered in 
the varicap circuit of the 35.075-MHz VCO which was 
causing the excessive static phase error. Replacement of 
the varicap with a new ultra-low leakage unit restored 
the circuit to correct operation. The circuit was moni- 
tored for the next 4 wk, and was found to be satisfactory 
(SPS 37-40 , Vol. Ill, p. 61). 

3. CFS Performance 

The CFS has completed 9,423 continuous hours of 
operation as of October 2, 1966. During this quarter, 
monthly preventative maintenance schedules were per- 
formed and system data was taken. Of the systems moni- 
tored, the phase noise measurement and VCO frequency 
drift are the most important. Figs. 54 and 55 show the 
four VCO aging rates, or drift, at 30-day intervals, and 
they compare very well with the crystal manufacturer’s 
tolerance for 90 days. It should be noted that all the RF 
components of the VCO within the +55°C oven are 
aging together and that the manufacturer’s aging curve 
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is for the quartz crystal only. Fig. 54 shows both 31.44- 
and 31.84-MHz VCO aging rates to be within the manu- 
facturer’s stated limits, and Fig. 55 shows that the 
35.075-MHz VCO has exceeded its limit during the July- 
August period. This relative increase could be due to the 
leakage current in the varicap circuit, as previously de- 
scribed. 

Very good correlation is obtained between the aging 
rate and leakage current of the varicap, because the 
aging rate falls below the manufacturer’s limits with the 
replacement of the faulty varicap. This data point is not 
shown on Fig. 55 because of the interruptions in oven 
stability and the replacement of the varicap in the 
35.075-MHz VCO only. Although the 35.2-MHz VCO 
aging curve exceeded the early limits (faster aging rate), 
it has returned to specification during the 90-day period 
(Fig. 55). 

The phase noise characteristics of all four phase-lock 
loops are shown in Figs. 56 and 57 which show the 
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Fig. 55. Frequency aging characteristics for 35.075- and 
35.2-MHz VCO 


standard deviation <r for a three-month period. The vari- 
ations in phase noise appear to be random, and two phase 
lock loops follow no particular trend. The 31.44-MHz 
phase-lock loop noise curve is not very meaningful, be- 
cause the data point for the month of August was not 
available; however, a trend here may be possible (Fig. 

56) . The 31.84-MHz phase-lock loop seems also to be 
following the same trend, while the 35.075- and the 
35.2-MHz phase-noise curves appear to be cyclical (Fig. 

57) . Observations during the next quarter should show 
any definite trends. 

4. Conclusions 

Installation of the CFS along with its associated tests 
and measurements has been completed, and the system 
has been officially released for operational use. DSS 13 
personnel have been schooled on operations, test pro- 
cedures, and preventative maintenance. Three sets of 



Fig. 56. Phase noise characteristics of 31.44- and 
31.84-MHz phase-lock loops 



Fig. 57. Phase noise characteristic of 35.075- and 
35.2-MHz phase-lock loops 
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complete documentation have been provided, and will Further discussions will be limited to failure modes, 
be updated when necessary. modifications, and improvements. 
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IV. Communications Engineering Development 


N67 15921 

A. S-Band Receiver-Exciter Ranging Subsystem 
(RER) r Blocks 1 1 B and MIC, H. Donnelly, C. E . Johns, 
R. W . Tappan, and M. R. Wick 

1. Introduction 

The original S-band receiver-exciter (RE) subsystem 
was designated as GSDS 1964, or Block I, (Ref. 1 and 2) 
and was designed specifically for DSN requirements. 
The RE subsystem with the ranging (R) subsystem 
Mark I (Ref. 3) can precisely determine ranges up to 
800,000 km when used in conjunction with an S-band 
“turnaround” transponder. 

The use of this S-band RE subsystem for the MSFN 
necessitated the addition of some new designs as well as 
changes in the existing design. The RE subsystem con- 
taining these design changes has been designated Block 
II (Refs. 4 and 5). 

Additional requirements for both the DSN and MSFN 
have necessitated further design improvements. These 
requirements are: 

(1) The capability of switching from one channel to 
another without the need to rephase and recali- 
brate the range receiver. 

(2) The capability of the receiver exciter ranging 
(RER) subsystem to support both DSN and MSFN 
tracking commitments. 
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The designs providing the capability of these require- 
ments are designated the RER Blocks IIB and IIIC, re- 
spectively. 

2. Receiver-Exciter Ranging 

The Block II RER subsystem uses a range code syn- 
chronizing signal (clock) that is derived from the exciter 
voltage-controlled oscillator (Fig. 1). This is mechanized 
by multiplying the VCO frequency by a factor of 5/221. 
The clock is also used as a reference signal for the rang- 
ing receiver and transfer loop. 

Generating the range clock and reference signals in 
this manner makes it necessary to recalibrate the ranging 
receiver each time a large frequency change is made, 
such as switching between channels. This is because the 
delay in the ranging receiver varies as the clock fre- 
quency is changed. The primary contributor to this delay 
change is the predeteetion bandpass filter within the 
ranging receiver IF amplifier. 

To eliminate the necessity for recalibration, the Block 
IIB design uses a fixed clock and reference frequency 
derived from the station standard. The design of the 
Block IIB ranging receiver and ranging subsystems is 
described in the following: 

a. Ranging receiver. The fixed clock and reference fre- 
quency establishes a relatively constant delay in the rang- 
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FROM BLOCK U TO II B 















ing receiver. Except for doppler offset, the detected ranging 
receiver clock frequency is constant, thereby minimizing 
the effect of the predetection filter. A fixed clock fre- 
quency of 496 kHz was selected to be compatible with 
most of the existing hardware. The clock signal is derived 
by multiplying a 1-MHz signal output from the tracking 
station standard by a factor of 496/1000. 

The subassembly that performs the 5/221 multiplica- 
tion in the Block II system to obtain the clock signal also 
generates two other coherent signals. These two signals 
drive a 57/221 frequency shifter subassembly to provide 
reference signals for the coherent translator and doppler 
extractor portions of the RF system. Since the 5/221 
output is not required when the fixed clock frequency is 
used, a new design 57/221 frequency shifter was de- 
veloped. This subassembly replaces both the Block II 
5/221 and 57/221 frequency shifter subassemblies. 

Fig. 1 is a block diagram showing both the present 
Block II ranging receiver together with the modifications 
which upgrade it to the Block IIB configuration. 

A block diagram of the 496/1000 frequency shifter, 
which generates the fixed 496-kHz clock signal, is shown 
in Fig. 2. A 1-MHz signal from the station standard is 
used as the input to the shifter. At the output of the 
1-MHz amplifier the signal is divided into two paths. 
One path connects the 1-MHz amplified signal to one 
input of a balanced mixer. The second path is used to 
synchronize an oscillator whose output is tuned to 200 


kHz. Two similar divide-by-five stages follow, generating 
an 8-kHz signal which is applied to the second input of 
the balanced mixer. The 992-kHz difference frequency at 
the mixer output is selected by a crystal bandpass filter. 
After two successive stages of amplification, the 992-kHz 
signal is used to synchronize a divide-by-two oscillator 
circuit generating the desired 496 kHz signal. 

The new design 57/221 subassembly utilizes a new 
technique for obtaining frequency division from a locked 
oscillator and is described in detail elsewhere in this 
volume. 

The 496-kHz fixed clock frequency is used in both the 
DSN and MSFN Block IIB design. To obtain a com- 
pletely common system, which is the purpose of the Block 
IIIC design, the three Block II DSN loop noise band- 
widths are being used for both the Block IIB DSN and 
MSFN receivers. This eliminates the need for including 
both the three DSN and three MSFN bandwidths and the 
capability of selecting either set of three. The DSN noise 
bandwidths (0.8, 4, 12 Hz) are adequate to meet the 
requirements of both the DSN and MSFN. 

b . Ranging subsystem , Mark l. The fixed ranging clock 
frequency of the Block IIB receiver-exciter causes a non- 
integer ratio to exist between the clock and RF doppler 
frequency. The doppler counting logic of the Mark I 
must be changed to allow the RF doppler to be tallied 
correctly. 


AMPLIFIER AMPLIFIER AMPLIFIER 



Fig. 2. 496/1000 Frequency shifter, block diagram 
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Fig. 3 illustrates the method of doppler counting used 
in the current ranging configuration. Both the clock and 
RF doppler signals pass through an input buffer to the 
range tally. The doppler selector routes the clock doppler 
to the range tally during code acquisition, and RF dop- 
pler to the range tally during the tracking portion of the 
ranging operation. The clock and RF doppler counts are 
related to each other by a fixed ratio. 72 RF doppler 
counts = 72 range units (RU) = 1 clock doppler count. 
The length of the range unit in light-time units is de- 
pendent upon the ranging clock frequency, which, in 
turn, is derived from the transmitter frequency F t , and is 


1 RU = 



light seconds, 
one-way range. 


In the Block IIB configuration, the ranging clock fre- 
quency is fixed at 496 kHz. The range unit is, therefore, 
a constant length, and is: 

1 RU = 2 X 288 X 496 kHz = 3 ' 5002 li8ht 

one-way range or 

approximately 1.049 m 


Fig. 4 illustrates the manner in which the clock and RF 
doppler are counted in the Block IIB configuration. The 
range tally retains the Chinese Number and clock dop- 
pler tallying functions. The RF doppler tallying is done 
by a new subassembly called the RF doppler counter. 
This subassembly contains the RF doppler input buffer, 
doppler counter, and range number adder. 

The range tally stores the portion of the total range 
that resulted from code shifting and clock doppler count- 
ing, and the RF doppler counter stores the range change 
that occurred since the last clock doppler count. The RF 
doppler counter counts the number of range units of RF 
doppler between clock doppler counts (Fig. 5). When the 
range to the spacecraft is increasing, a clock doppler 
count causes 72 RU to be added to the range tally, and 
the RF doppler counter to be reset to zero. RF doppler 
counts then accumulate in the RF doppler counter until 
the next clock doppler count resets it to zero. If the range 
to the spacecraft is decreasing, a clock doppler count 
causes 72 RU to be subtracted from the range tally, and 
the RF doppler counter to be set to 72. The range number 
adder sums the outputs of the range tally and the RF 
doppler counter to provide an output that is the total 
range. 


CLOCK 

DOPPLER 


RF 

DOPPLER 


Odeg 
90 deg 


Odeg 

90deg 


SHIFT CONTROL TO 



Fig. 3. Ranging subsystem Mark I doppler tally, Block II, configuration 
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Fig. 5. Timing diagram 


A small range error may exist in the Block IIB ranging 
configuration. The ratio of the ranging clock frequency 
to the S-band transmitter frequency is not constant; there- 
fore, one clock doppler count is not always equal to 


exactly 72 RF doppler counts. Since this error is not 
cumulative, it will always be less than 0.5 RU. Tech- 
niques for removing this error are known and may be 
used if the error is considered significant. 

3. Receiver-Exciter, Block IIIC 

The purpose of the Block IIIC design is to fulfill re- 
quirement 2, the capability of tracking both DSN and 
MSFN space probes. The existing Block II design does 
not provide this capability in the following areas: 

(1) Exciter tuning range 

(2) Receiver tuning range 

(3) Reference receiver predetection and loop noise 
bandwidths 

(4) Receiver AGC loop noise bandwidth 

(5) Range receiver predetection and loop noise band- 
widths 
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Fig. 6. Block II and MIC block diagram 
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Following is a more detailed discussion of these areas. 

a. Exciter tuning range. The existing X32 frequency 
multiplier and UHF buffer amplifier of Block II are 
capable of delivering the required drive level to the 
power amplifier over either the DSN or MSFN 10-MHz 
bands but not over the combined 20-MHz band. In Block 
IIIC (Fig. 6) the X32 frequency multiplier is replaced 
with a design that has adequate bandwidth to cover both 
the DSN and MSFN bands. This is accomplished in two 
subassemblies, a X 4 and a X 8 frequency multiplier. The 
-1-27.5 dbm minimum output level of the Block IIIC 
multiplier chain provides adequate drive to obtain the 
20-MHz bandwidth from the existing Block II UHF buf- 
fer amplifier. 

The X4 frequency multiplier subassembly consists of 
two cascaded low-level balanced diode doublers sepa- 
rated by an isolation amplifier and followed by four 
stages of amplification. A coaxial commercial bandpass 
filter is installed between the first two of the four cas- 
caded amplifiers to provide rejection of harmonics of the 
input frequency other than the X4. 


has been added to the subsystem ahead of the 10-MHz 
IF amplifier. This unit is mounted in the location now 
occupied by the 0-12 db attenuator subassembly (gain 
adjustment at 10 MHz now being provided only in the 
angle channel receivers). Since both the DSN and MSFN 
predetection crystal filters are located in the 10-MHz filter 
amplifier subassembly, the crystal filter in the Block II 
10-MHz IF amplifier has been removed. 

The receiver loop filter module contains the time con- 
stants for both the three DSN and the three MSFN loop 
noise bandwidths. 

d. Receiver AGC loop noise bandwidth. Three loop 
time constants (380, 34, 4 sec) which adequately cover 
the DSN and MSFN requirements have been selected. 
These bandwidths are common to both systems, and no 
selection is required. 

e. Range receiver predetection and loop noise band- 
widths. As in the AGC loop, a design common to the 
DSN and MSFN has been selected for the ranging receiver 
(see Block IIB design). 


The X8 frequency multiplier uses some of the latest 
techniques in step recovery diode frequency multiplica- 
tion, and is discussed in detail in Ref. 6. 

b. Receiver tuning range . The problem of adequate 
tuning range in the receiver local oscillator for both 
DSN and MSFN combined bands has been resolved in 
Block IIIC by replacing the local oscillator multiplier 
chain with the same design used in the exciter. The use 
of a second coaxial filter in the X4 multiplier subassem- 
bly at the output of the final power amplifier provides 
the additional attenuation of sidebands that is required. 
With this design it is no longer necessary to provide a 
filter at the S-band output of the local oscillator multiplier. 

Adequate tuning range at the receiver input is obtained 
by replacing the existing preselector with one that covers 
the combined DSN and MSFN bands. 

c. Reference receiver predetection and loop noise 
bandwidths . The Block IIIC design contains the prede- 
tection and loop noise bandwidths of both the DSN and 
MSFN Block II subsystems. Selection is made independ- 
ently for each receiver with the switches on the acquisi- 
tion control panel (Ref. 7). 

To provide both DSN and MSFN predetection filter 
bandwidths, a new 10-MHz filter amplifier subassembly 


4. Conclusion 


Blocks IIB and IIIC provide a system with both the 
DSN and MSFN parameters and the capability to select 
these parameters by switch control. The parameter dif- 
ferences have now been reduced to the reference receiver 


predetection and loop noise bandwidths. However, the 
exciter and receiver voltage-controlled oscillator sub- 
assemblies still contain only four selectable channels. 
Two sets of VCO subassemblies are provided in Block 
IIIC, one for the four MSFN channels and one for four 
DSN channels. Switching from an MSFN channel to a 
DSN channel or vice versa requires that the VCO sub- 
assemblies installed in the subsystem be replaced with 
the second set. Present plans are to upgrade all DSN 


and MSFN RER subsystems to 

N6f 15 


I, IIB and IIIC. 
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1. Introduction 


A frequency shifter for the Block IIB receiver-exciter 
has been developed to derive reference signals coherently 
related to the exciter VCO frequency by a factor of 
57/221. These reference signals are required for the 
doppler extractor and the coherent translator. 


The block diagram of the 57/221 frequency shifter 
subassembly is shown in Fig. 7. The input signal from 
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22 MHz INPUT 




FOLLOWER FOLLOWER 


Fig. 7. 57/221 Frequency shifter block diagram 


the exciter VCO is branched into two paths. One signal 
path is used to synchronize a transistor oscillator circuit 
whose output is tuned to Math of the input frequency /*. 
The output frequency of this synchronous oscillator, 
which is 17/221 of the input frequency, is doubled in 
a transistor stage and coupled to one input of a balanced 
mixer. The second signal path is used for synchronizing 
a similar oscillator whose output is tuned to With of the 
input frequency. This signal, 13/221 of the input fre- 
quency, is multiplied by seven in a transistor stage and 
applied to the second input of the balanced mixer through 
an isolation amplifier. The difference frequency at the 
mixer output, 57/221 times the input frequency to the 
subassembly, is selected by the crystal bandpass filter. 
Amplifier stages provide two isolated outputs. 

Emitter followers precede the oscillators to prevent 
signal cross-coupling between their inputs and to reduce 
the amount of feedback to the input terminal. The output 
of each oscillator is impedance-matched to its respective 
multiplier by means of another emitter follower. 

2. Synchronized Avalanche Oscillator 

Basic to the operation of the frequency shifter is the 
type of synchronized oscillators used to perform the 
division processes. Although schematically the oscillators 
appear to be similar to the conventional type, their oper- 
ation is entirely different. 

In the conventional type the natural frequency (tuned 
circuit resonant frequency) is the desired frequency. It is 
difficult to force the circuit to oscillate at frequencies 
higher or lower than its natural frequency by an external 
signal and, therefore, it can be synchronized only over a 
narrow frequency range. During the development of 
the 57/221 frequency shifter, it was found that this type 
of oscillator would not perform the high division ratios 
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over the required frequency range. Another character- 
istic of the conventional locked oscillator is its relatively 
large amount of phase instability. Since the natural fre- 
quency of the oscillator is the desired frequency, syn- 
chronization is maintained by supporting this natural 
oscillation with the incoming signal. This does not occur 
at precisely the same phase in each cycle, and phase 
instability occurs. In addition, oscillator drift with tem- 
perature change causes loss of synchronization unless 
compensation is added. 

To meet the requirements of the Block IIB 57/221 
frequency shifter design, it was necessary to investigate 
other means of synchronizing an oscillator circuit. The 
use of avalanche breakdown as a method of synchroniz- 
ing was explored. A schematic of a divider circuit using 
this technique is shown in Fig. 8. The resonant circuit is 
formed by L2 and the series combination of C2 and C3. 
Regeneration is accomplished by the connection of the 
emitter to the junction of C2 and C3. 

Fig. 9 is a composite oscilloscope photograph showing 
operational events in their correct phase and voltage 


+ 17 vdc 
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Fig. 9. Synchronized avalanche oscillator waveforms 


relationships. Between t x and t 2 the emitter voltage 
(trace A) is positive with respect to the base voltage 
(trace B), therefore the transistor is not conducting. 
During this time the collector voltage (trace C) is in- 
creasing positively at a sinusoidal rate. At £ 3 the base- 
emitter bias voltage crosses over and starts forward 
biasing of the transistor. At £ 4 the bias becomes sufficiently 
large to cause the collector current to be driven into 
an avalanche mode of operation. The type of transistor 
used was selected for this characteristic. Between £ 4 and 
£ 5 the collector-to-emitter voltage rapidly decreases to 
zero and all the available supply voltage appears across 
the collector load. Between £ 5 and U no collector current 
flows as the collector-to-emitter voltage is zero. At the 
collector voltage again starts increasing sinusoidally, and 
the next cycle begins. 

When an input signal, which is some multiple of the 
oscillator frequency, is applied to the base of the oscil- 
lator, synchronization takes place. This occurs because 
one cycle of the input frequency (£ 4 on Fig. 9) adds to 
the base to emitter bias and forces the transistor into the 
avalanche mode. Capacitor Cl reduces the amount of 
ripple on the base, which is caused by the large collector 
variations. This permits the input frequency to be more 
effective in triggering the collector current. Inductor LI 
prevents the input signal from being shorted to ground 
through capacitor Cl. Its value is selected lo give a high 
impedance at the input frequency and a relatively low 
one at the oscillator frequency. 

Collector current flows only between the period of 

and £5 which represent approximately 10% of one cycle. 
During the remaining 90% no current flows, therefore the 
transistor is acting only as a regenerative switch. The 
transistor is prevented from jumping into a destructive 
breakdown mode by limiting the maximum collector cur- 
rent flow with an emitter resistor. During the avalanche 


conditions the collector-to-emitter voltage rapidly de- 
creases to zero, simultaneously reducing the current. 
Measurements indicate a peak-to-peak current over a one- 
cycle period of approximately 6.8 ma. 

3. Lock Range 

The avalanche-type locked oscillator depends on the 
breakdown point as well as its tuned circuit to establish 
the output frequency, making it more easily controlled 
by an external signal. The waveforms of a divide-by-five 
locked oscillator are shown in Fig. 10 in more detail to 
demonstrate this. A divide-by-five circuit was arbitrarily 
selected to give a clear presentation of the base wave- 
form (trace 2). The collector waveform is trace 1, while 
the large sine wave shown in dotted lines (trace 3) is the 
frequency of the collector tuned circuit. At point A on 
the base waveform, the collector current is triggered into 
the avalanche mode and causes the output frequency 
to be shown in trace a. If the input frequency is de- 
creased, then the collector avalanches at point B (on the 
base waveform) and the output frequency is decreased 
(as shown on trace b of the collector waveform). Con- 
versely, if the input frequency is increased (point C on 
the base waveform) the output frequency is increased 
(trace c on the collector waveform). 

Control of the output frequency is relatively inde- 
pendent of the collector tuned circuit. Consequently, 
it is not necessary to pull the natural oscillator frequency 
when changing the frequency of the input synchronizing 
signal. The avalanche oscillator, therefore, inherently has 
a much larger lock range than the conventional type. 
Measurements of two divide-by-five circuits (28.35 to 
5.67 MHz), one of each type, show the lock range for 
the avalanche oscillator to be more than four times that 
of the conventional oscillator. The conventional type 
was optimized for lock range prior to performing the 
measurements. 



Fig. 10. Typical divide-by-five circuit waveforms 
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4. Phase Instability 

Controlling the synchronism of the oscillator by ava- 
lanche breakdown rather than supporting the natural 
frequency of the oscillator gives a more positive lock. 
This has been demonstrated in phase instability measure- 
ments made on both types. To accomplish this, two iden- 
tical avalanche divide-by-five circuits (28.35 to 5.67 MHz) 
were compared, using a 5.67-MHz phase detector. The 
two divide-by-five circuits were driven by the same 
signal generator to eliminate the instability of the source. 
Then an avalanche and a conventional divider were 
compared in the same manner. The measurements indi- 
cated that the phase instability contributed by the con- 
ventional locked oscillator divider circuit is approximately 
four times as great as the avalanche type. The contribu- 
tion of the avalanche oscillator in the second measure- 
ment is negligible so that this measurement represents 
approximately the instability of the conventional locked 
oscillator. 

5. Conclusion 

Because of the method of synchronization, the ava- 
lanche oscillator has demonstrated several advantages 


over the conventional locked oscillator. The advantages 
are that it: 

(1) maintains synchronism over a greater frequency 
range, 

(2) has less inherent phase instability, and 

(3) is less sensitive to temperature changes. 

Typical characteristics of the 57/221 frequency shifter 
subassembly using the synchronized avalanche oscillator 
design are: 

Input frequency: 21.9792 MHz 

Output frequency: 5.6688 MHz 

Bandwidth (controlled by crystal filter): ±33 kHz 
(~ldb) 

Lock range: > ±2% 

Phase instability: approximately 0.01 deg rms 
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V. Deep Space Station Engineering and 
Operations 


N67 15923 

A. Flight Project Support, j. Orbison 

1. Surveyor II Mission 

The Surveyor II spacecraft was launched from Cape 
Kennedy, Florida, on September 20, 1966. On the same 
day, the Pioneer DSS, with the Echo DSS as a trans- 
mitter backup, tracked the first pass. The midcourse 
maneuver commands were transmitted from the Pioneer 
DSS early in the pass. 

Following the second pass acquisition of the 
Surveyor II spacecraft on September 22 by the Pioneer 
DSS, unsuccessful attempts were made to stop the vehi- 
cle's tumbling, and less than 24 hr before the scheduled 
landing of the spacecraft, the mission was terminated. 

2. Lunar Orbiter I Mission 

Tracking of the Lunar Orbiter I was performed daily 
by the Echo DSS until pass 36, on September 16, 1966. 
During pass 34, on September 14, the Goldstone calibra- 
tion leader in the spacecraft camera equipment was read 
out in a final video sequence for a TV system degrada- 
tion determination. 

Following pass 36, tracking of the Lunar Orbiter I has 
been on an extended mission, as-scheduled basis. Rang- 


ing and time correlation experiments using the Mark I 
ranging subsystem and the Lunar Orbiter ranging trans- 
ponder have been performed between the Echo DSS and 
the Robledo DSS, and the Echo DSS and the Tidbinbilla 
DSS. Two Stanford lunar mapping experiments using 
the Lunar Orbiter I spacecraft were performed during 
passes 56 and 62. 

3. Pioneer Missions 

a. Pioneer VI mission . Currently in the tenth month of 
its mission, the Pioneer VI spacecraft is being tracked by 
the Mars DSS 210-ft antenna on a scheduled basis. The 
Echo DSS provides telemetry and command processing 
via microwave to the Mars DSS. Routine data are being 
returned from the spacecraft, and the vehicle provides 
an accurate source for antenna and system testing in 
progress at the Mars DSS. 

b. Pioneer VII mission . Because of the Echo DSS com- 
mitment to Lunar Orbiter J, the Pioneer VII spacecraft 
was tracked by the Pioneer DSS from the launching date, 
August 17, 1966, until pass 22, September 7, 1966. 
Telemetry and command processing from the Pioneer 
ground operations equipment (GOE) and the telemetry 
and command processing (TCP II) beta computer was 
provided by the Echo DSS to the Pioneer DSS (SPS 
37-41, Vol. Ill, p. 116). 


JPL SPACE PROGRAMS SUMMARY 37-42, VOL III 


87 



With pass 24, September 9, 1966, the Echo DSS as- 
sumed the Pioneer VII prime tracking on a daily basis, 
until pass 46 on October 1, 1966. An extended mission 
tracking on an as-scheduled basis is currently being per- 
formed three to four days a week. 

4. Mariner IV Mission 

The Mariner IV spacecraft, currently in the twenty- 
third month of its flight, continues to be tracked jointly 
by the Mars DSS and Venus DSS. During a scheduled 
monthly track, two-way lock-up is established, com- 
mands are transmitted, and telemetry is received and 
recorded. The Mars DSS performs the tracking, necessary 
command transmissions, and receiving, using the JPL 
R&D maser and receivers (SPS 37-39, Vol, III, p. 112). 
The data are microwaved to the Venus DSS for proces- 
sing. During October 1966, a series of low-level telemetry 
experiments were performed by the two stations (SPS 
37-42, Vol. Ill, p. 92). 


5. Atlas/Centaur Project 

The ninth in the series of Atlas/Centaur second-stage 
booster tests (A/C 9) was launched from Cape Kennedy, 
Florida, on October 26, 1966. The Pioneer DSS partici- 
pated in the test, tracking the spacecraft for 4.5 hr during 
the first pass on October 27. Recorded telemetry data 
were forwarded to JPL for processing. 

N67 15924 

B. Facility Construction and Equipment 
Installation, J. Orbison 

1. Pioneer Deep Space Station 

The MSFN equipment installation and system testing 
are progressing. The RF switching box, enabling rapid 
transfer of the antenna and antenna-associated equip- 
ment, has been installed and operationally tested. The 
tower and antenna for the microwave link between 
DSS 11 and the Apollo Station are installed and operating. 
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Trenching for the underground primary power distri- 
bution system, temporarily halted during the Surveyor II 
mission, is again in progress. All primary power cables, 
to be placed in conduit in concrete tunnels accessible 
through strategically located manholes, will lead to both 
wings of the control building and other buildings of the 
station. 

Installation of the waveguides and the RF power com- 
biner for the dual transmitter (SPS 37-41 , Vol. Ill, p. 
122-124) on the antenna is in progress. Masers 1 and 2 
have been installed on the antenna, and maser 1 is fully 
operational in the GSDS configuration. Maser 2, which 
replaces the S-band system parametric amplifier, is cur- 
rently in the final assembly stage in preparation for 
operational testing. 

Equipment of the S-band system requiring certification 
was routed through the Goldstone calibration, certifica- 
tion, and maintenance facility, following the termination 
of the Surveyor I mission. Many of the Pioneer DSS 
S-band receiver modules are the original pilot modules, 
which were calibrated and certified to establish oper- 
ational levels comparable with the latest S-band systems. 
Concurrently, to provide an area for the new equipment, 
relocation of racks at the multiple mission support area 
(MMSA) continued. At present, the two racks for the 
telemetry and command data (TCD) handling are in- 
stalled. Fig. 1 illustrates a typical control room config- 
uration, including the MMSA. 

2. Echo Deep Space Station 

The M MSA/mission support recording (MSR) equip- 
ment has been installed and operationally tested. The 
rack rearrangement of the microwave, recording, and 
analog instrumentation subsystems to provide area for 
the installation of the MMSA/MSR also necessitated 
moving the communications teletype equipment to the 
alcove of the control room. 

Construction of the commercial power line from the 
substation north of the Fort Irwin Army Air Field to the 
Tiefort Peak collimation site began in September 1966. 
Approximately 2 mi of poles have been set across the 
relatively flat part of the desert, and drilling for the 
pole-line up the ridge to the collimation site is in prog- 
ress. Because of the extremely rugged and steep terrain 
(Fig. 2), men and equipment are transported to the work 
area by helicopter. The ridge poles will be set and the 
cable strung in the same manner. 



Fig. 2. Aerial view of the Tiefort Peak collimation site 
and power line area 


3. Mars Deep Space Station 

Except for the maser and associated equipment not 
presently in use, installation and operational testing of 
the S-band system are nearing completion. Operational 
testing of the subsystems is taking place during the 
Pioneer VI tracking. Station personnel are assisting in 
the hydrostatic bearing alignment, the master equatorial 
alignment and calibration, and the installation and align- 
ment of the auto collimater in the independent reference 
structure (IRS) which is associated with the master 
equatorial. 

The transfer from station power to commercial power 
at 10:10 PDT on October 19, 1966, marked the completion 
of the commercial power installation. Construction of the 
underground cabling (Fig. 3) from the pole-line terminal 
to the substation had been completed earlier. 
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Fig. 3. Mars DSS underground commercial power line 
construction 


The backup 1126-kw diesel generator brought to the 
Mars DSS during the period that the three 500-kw diesel 
generators were being rebuilt (SPS 37-41 , Vol. Ill, p. 121) 
was shut down and prepared for shipment. Investigation 
and testing of the torque problem in the 500-kw units 
continue; however, two 500-kw units are currently avail- 
able as emergency backup to the commercial power. 

4. Venus Deep Space Station 

Construction of the laboratory and office building 
(G-60) is essentially completed (Fig. 4) and occupancy 
took place late in October. 

N67 15925 

C. ' Mariner Venus 67 Project: Ground 

Instrumentation for Occupation, G. Levy 

Although the ground instrumentation for the Mariner 
Venus 1967 project will be similar to that used for 
Mariner IV, there will be marked differences due to the 
very dense nature of the Cytherean atmosphere. The 
prime station for the occultation measurements will be 
the DSS 14 210-ft advanced antenna system (AAS). 



Fig. 4. Completed laboratory and office building at Venus DSS 
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LOOP FILTER 



Fig. 5. Mariner Venus 67 D$$ 14 occupation configuration 































DSS 11 will receive doppler data in the normal config- 
uration, and DSS 13 will have an open-loop receiver. 

Because of the high atmospheric density of Venus, 
doppler perturbations of the order of 30 kHz are ex- 
pected (in contrast to 3 Hz, one way, in the case of the 
Mariner IV). In addition, the atmosphere also acts as a 
defocusing lens which produces refractive attenuation 
estimated to be as high as 26 db. To ensure the most 
favorable possible signal-to-noise ratio, the 210-ft AAS 
site has been chosen as the prime station. 

The preliminary Mariner Venus 67 DSS 14 occultation 
configuration (Fig. 5) shows three separate receivers being 
employed. Receiver one will be the standard phase lock 
loop DSIF configuration and will be used for telemetry 
and doppler data acquisition. Receiver two will be a 
slightly modified standard receiver with fixed tuned local 
oscillator derived from the frequency standard by using 
a synthesizer. The 10-MHz IF of this open loop receiver 
will be frequency-translated to the 10- to 200-kHz region 
and recorded. This bandwidth is necessary because of 
the uncertainties in the expected signal as well as antici- 
pated frequency shifts due to the planetary atmosphere 
and the orbit. In addition to the two DSIF receivers, a 
third receiver, the alidade instrumentation area receiver, 
is illustrated. This receiver has a programmable local 
oscillator which may be used for the occultation experi- 
ment. 

At present a simulator is being constructed to produce 
the expected signal frequency and power fluctuations. In 
addition, DSS 11 will obtain telemetry and doppler data 
and DSS 13 will have an open loop receiver and a special 
low noise “Ultra-cone.” The Ultra-cone will be a listen- 
only configuration and will have a zenith system tempera- 
ture of 20° K or less. 
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D. Venus DSS Operations, M. A. Gregg, E. B. Jackson, 
and A. I . Price 

1. Experimental Activities 

During the period from August 16 through October 
14, 1966, the major activities at the Venus DSS (sum- 
marized in Table 1) were bistatic (with the Mars DSS) 
planetary radar experiments, radio source evaluation 
(using the Venus DSS 85-ft AZ-EL antenna), and recep- 
tion from and transmission to the Mariner spacecraft. 
Additionally, certain improvements were effected in the 
station antennas, buildings, and equipment. 
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Table 1. Summary of Venus DSS activities 
(August 1 6— October 14, 1966) 


Activity 

Hours 

Per cent 

Primary experiments 

Planetary radar (Venus bistatic) 

49 

3.40 

Secondary experiments 

Radio source temperature calibration 

608 

42.22 

Mariner experiments 

Transmission, reception, and testing 
(including time required to change 
Cassegrain feed cone) 

50 

3.47 

Testing, calibration, construction, and 
maintenance 

709 

49.24 

Holidays and scheduled nonoperating time 

24 

1.67 

Total 

1440 

100.00 


The cooperative, with Mars, tracks of the Mariner 
spacecraft continued to return telemetry data, spectro- 
grams, and some doppler measurements. These tracks 
usually involve reception at the Mars DSS, with any 
necessary transmissions being accomplished by the Venus 
DSS. However, a low-level, low-speed telemetry simula- 
tion experiment with the Mariner spacecraft required 
weaker signals than were available from the Mars an- 
tenna. To accomplish this experiment, the Venus R&D 
Cassegrain feed cone, with the maser retuned to 2297.5 
MHz, was used for reception of the Mariner spacecraft in 
both the one- and three-way modes. Signal strength mea- 
sured at the Venus DSS in this mode was —177.6 dbm. 

In addition to the signal strength measurements made 
on the Mariner spacecraft during the reception at Venus, 
some data were obtained on spacecraft-transmitted fre- 
quency drift rate in the one-way mode. These data were 
a series of 1-min integrations taken at 4-min intervals 
and plotted against frequency using a 29.3-Hz sampling 
bandwidth. The results (shown in Fig. 6) indicate that 
the frequency received by the Venus antenna drifted 
downward approximately 1.7 Hz in 20 min. Fig. 7(a) 
illustrates a spectrum of the spacecraft-received fre- 
quency and shows the first pair of sidebands as well as 
the carrier. Fig. 7(b) illustrates the carrier alone, inte- 
grated over a 31-min period. 

The planetary radar experiments, with Venus as the 
target, were the first attempted by JPL with Venus ap- 
proaching superior conjunction at a one-way range of 
approximately 1.632 AU. Total spectrum and ranging 
measurements were made using a range gate depth of 
1000 /isec and disclosed that the planet Venus was 
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Fig. 6. Manner Mars 1964 spacecraft, carrier only, showing frequency drift in one-way mode 


approximately 9000 /*sec (round trip) closer than pre- 
dicted by the JPL master ephemeris. However, because 
of insufficient signal strength, accurate measurements 
could not be made. 

2. Subsystem Performance 

a. 100-kw transmitter (operation). During this period, 
the R&D transmitter operated a total of 239.7 filament 
hours and 110.4 beam hours for the Mercury and Venus 
bistatic radar experiment. A few weeks were devoted to 
intensive training of personnel in operations and mainte- 


nance of the transmitter systems and the associated 
equipment, both S-band and X-band. 

b. Mariner 100-kw transmitter. The Mariner transmit- 
ter has operated a total of 33.8 filament hours and 25.8 
beam hours with no inter-lock operations or loss of 
tracking time. 

c. Receiving systems. Although the X-band receiver 
(8448 MHz) was not used during this period, the 2388- 
MHz receiver was used in planetary experiments and 
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Fig. 7. Mariner Mars 1964 spacecraft, one-way mode, carrier plus sidebands and carrier alone 


radio source evaluation. The receiver was used in its 
monostatic mode at the Venus DSS and in its bistatic 
mode in conjunction with the Mars DSS 2388-MHz R&D 
receiver. In the latter, the 455-KHz signal from the Mars 
DSS receiver is fed, via the microwave link, into the AM 
channel and the main loop channel of the Venus DSS 
2388-MHz receiver. 

Reception of signals from the Mariner spacecraft was 
accomplished using both the Venus DSS Mariner receiver 
(for one-way and three-way reception) and the Mars DSS 
2295-MHz receiver (for one-way and two-way recep- 
tion). As before, the 455-kHz signal from the Mars re- 
ceiver was fed, via the microwave link, into the 455-kHz 
portion of the Venus Mod IV receiver. 

With the exception of one failure which involved the 
loss of output from the bias-up synthesizer, operation of 
the programmed local oscillator remained satisfactory 
throughout this period. The failure was repaired on-site, 
and normal operation was restored. 

Two failures occurred in the central frequency syn- 
thesizer (CFS). One was concerned with the + 30-v power 
supply voltage vigilant unit which senses large changes 


in the power supply output voltage. Although trouble- 
shooting is now underway, it should be noted that the 
failure of this unit does not affect the power supply 
output. The other failure resulted in excessive static 
phase error in the 35.075-MHz reference loop. A faulty 
voltage-controlled oscillator was repaired, and operation 
was returned to normal. 

3. System Improvements 

a. Antennas. To improve appearance and to provide 
additional storage space an enclosure to contain moving 
portions of the cable wrap-up has been added to the 
30-ft AZ-EL antenna. Also, the complete antenna struc- 
ture, with the exception of the surface of the dish, has 
been painted white to improve appearance and to reduce 
structural deformation due to thermal gradients. 

During the month of November, operation with an 
X-band cone at a frequency of 8488 MHz will commence 
on the 85-ft antenna. To improve performance at this 
frequency, thermal gradients have been minimized by 
complete repainting of the antenna. 1 

According to JPL Specification 1006B (modified). 
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b. 100-kiv transmitter . The solid-state rectifiers for the 
100-kw S-band transmitter were removed and replaced 
with the tube rectifiers to allow the testing of the indi- 
vidual diodes in the stacks and replacement of defective 
diodes before shipping the stacks for installation in the 
new solid-state rectifier equipment to be used for the 
400-kw system. Some tube rectifiers were lost through 
improper cooling and were replaced with existing spares. 
The inter-lock circuitry was modified to prevent this from 
re-occurring. 

c. Receiving systems . To provide dual frequency re- 
ceiving capability in one antenna cone, the 2295- to 30- 
MHz converter and power supply were removed from the 
Mariner cone and installed in the 2388-MHz R&D cone. 
Appropriate cabling was also installed, thus providing 
the two converters independent control and signal func- 
tions. A 400-cycle power cable change is the only modi- 
fication required to switch reception from one frequency 
to another. 

Because of continuing tape-reader problems in the 
digital portion of the programmed local oscillator, a com- 
plete mechanical and electrical overhaul of the tape 
reader and tape spooler was performed during this 
period. Modification of the tape-reader lamp power sup- 
ply greatly improved operation of the reader. 
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E. DSIF FR-900 Video Recorders, W. c. Allen 

Three DSIF stations— DSS 12, DSS 61, and DSS 41- 
are equipped with Ampex FR-900 rotary-head magnetic 
recorders capable of acquiring video data from the Lunar 
Orbiter. The machine presently installed at each site has 
a single wideband channel and two auxiliary tracks. The 
wideband can record either one channel of dc to 5.5-MHz 
data in an FM mode or one channel of data within a 
1- to 10-MHz bandpass in a direct-record mode without 
bias. The Lunar Orbiter video data are acquired from an 
unrestricted 10-MHz output of the S-band receiver at 
approximately —10 dbm, amplified to +3 dbm and 
predetection-recorded directly on the FR-900. Each of 
the two auxiliary tracks has a 30-kHz bandwidth, on 
which are recorded a station voice-communication chan- 
nel and serially coded GMT. A recording of approxi- 
mately 1 hr can be made on 3800 ft of the 2-in. -wide 
video tape. 

The FR-900 is a complex electromechanical device in- 
corporating three servos that control (1) the head drum 
rotation, (2) the longitudinal tape speed, and (3) two 
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electrically variable delay lines. During record, the 2-in. 
video tape is moved at a constant speed of 12.5 in./sec 
from the supply reel to the take-up reel by means of a 
capstan drive assembly, while wideband data are being 
recorded transversely across the width of the tape by 
either a rotating scanning assembly or a head drum. The 
capstan and head drum motion are separately controlled 
by hysteresis synchronous motors driven by signals de- 
rived from a crystal-controlled frequency standard. The 
wideband response is achieved primarily by a head-to- 
tape speed of over 1500 in./sec as a result of the rotary- 
head speed of 14,640 rpm. 

When a recorded tape is being reproduced, it is essen- 
tial that the longitudinal motion of the tape and the 
transverse motion of the rotating heads be tightly con- 
trolled so that the recorded tracks are accurately scanned. 
To accomplish this, a control track is recorded along one 
edge of the tape to mark each revolution of the head 
drum during the recording process. During reproduction, 
the head drum can then be rotated at constant velocity 
and the capstan speed can be servo-controlled so that 
one control-track pulse is reproduced for each revolu- 
tion of the head drum. This reestablishes the ratio of 
longitudinal-to-transverse motion that existed during re- 
cording and assures accurate registration of the repro- 
duce heads on the recorded tracks. 

An additional time-base reference or pilot signal is 
recorded transversely along with wideband data. When 
reproduced, this pilot is used as a time-base reference 
for control of head-drum rotation and permits additional 
electronic time-base correction of the reproduced data 
through electrically variable delay lines. The pilot is re- 
corded below the information spectrum at 500 kHz to 
avoid interference with the video signal. Precision time- 
base control is a principal feature of the FR-900 recorder/ 
reproducer system that makes possible the recording and 
reproducing of time-stable wideband information. When 
data are reproduced, the peak-to-peak time-displacement 
error is less than ±25 nsec with respect to the system- 
timing reference. The two auxiliary channels are re- 
corded in the conventional longitudinal manner along 
one edge of the tape, using a direct-record-with-bias 
technique. 

The Lunar Orbiter video data recorded at 10-MHz 
center frequency has a peak phase modulation index of 
±4 rad. Information sidebands extend more than ±1.24 
MHz about the 10-MHz carrier frequency and include a 
310-kHz video subcarrier modulated with video data 
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confined to a 230-kHz passband, a 38.75-kHz demodu- 
lator pilot reference, and a 30-kHz engineering data sub- 
carrier. 


During the first Lunar Orbiter mission, approximately 
300 magnetic tapes containing the above data were suc- 
cessfully recorded by the DSIF stations using the 10-MHz 
predetection record mode. Tapes recorded by the DSIF 
are processed at Langley Research Center, Hampton, 
Virginia, where film records are made for routine anal- 
ysis; and magnetic tape dubs of the detected video are 
made for digital analysis at the Manned Spacecraft Cen- 


ter, Houston, Texas. 
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F. DSIF Station Control and Data Equipment, 


R. Flanders, G. Jenkins , A. Burke, W. Frey, and H. Baugh 


The DSIF station control and data equipment (SCDE) 
was introduced and described in SPS 37-41 , Vol. Ill, pp. 
124-127. The Space Programs Summaries will continue 
to be used as the reporting medium for the status of 
SCDE implementation; this article presents the status as 
of October 1966. 


The SCDE consists of the following specific subsys- 
tems: (1) antenna pointing subsystem, (2) digital instru- 
mentation subsystem, (3) station control and monitor 
console, (4) telemetry and command processor, and (5) 
frequency and timing subsystem. The functional objec- 
tive of the equipment is to provide the capability for 
processing spacecraft telemetry and command data in 
direct mission support, to perform station monitoring, 
time tag data, supply reference frequencies, and to per- 
form antenna pointing as mission-independent support. 

1. Antenna Pointing Subsystem 


original program with the following additions and 
changes. 

Drive-tape reasonable-data-check option . The new 
program performs a check on the data on the drive-tape 
input. This new feature is a program option which may 
or may not be used, depending on the setting of com- 
puter breakpoint switch 2. If it is turned to RESET, the 
program performs the reasonable data check; if it is 
turned to SET, the program does not perform the rea- 
sonable data check. 

The reasonable data check looks at the angular data 
on the drive-tape samples and determines if it is accept- 
able for use in pointing the antenna. The program per- 
forms at straight-line extrapolation on existing data points 
to predict a new data point. When the new data sample 
is read into the computer, it is compared with the pre- 
dicted value. The new angle data must agree with 
predicted data with an allowable tolerance of the antenna 
pointing beamwidth or it is rejected by the computer 
program. If a sample is rejected, the program will search 
the drive tape for the next valid sample before it will 
process the angular data contained on the tape. This 
modification to the program will eliminate the possibility 
of pointing the antenna to erroneous angles that are in- 
troduced on the drive tape as a result of teletype trans- 
mission errors when the antenna-pointing prediction 
information is sent to the DSIF stations. 

Drive-tape input-format option. The new program has 
the option of reading either a drive tape generated by 
the teletype predict reformatting program (DOI-5004-OP) 
or the teletype tape as it is received by the station with- 
out using the reformatting program. The type of tape 
that will be accepted by the IAPS program is controlled 
by computer breakpoint switch 1 as follows: 


The antenna pointing subsystem (APS) has been in- 
stalled and is operational in its interim configuration at 
DSS 12, 14, 41, 42, and 62. Installation of interim APS 
(IAPS) equipment at DSS 62 and 72 will be completed 
during the fourth quarter of 1966. The hardware and 
software capabilities of the IAPS are described in detail 
in SPS 37-38 , Vol. Ill, pp. 74-76. The following addi- 
tional tasks have now been accomplished. 

a. Update of IAPS computer program . An updated 
version of the IAPS program (DOI-5107-OP-A) has been 
submitted to the DSIF computer-program library. This 
new program will replace the existing IAPS program 
(DOI-5107-IN), performing the same functions as the 


If breakpoint switch 1 is at RESET, the IAPS accepts 
a drive tape that has been generated by the teletype 
reformatting program; if it is at SET, the IAPS accepts a 
drive tape as it is received by the station without using 
the teletype reformatting program. 

This option allows the IAPS operator to bypass the 
time-consuming process of using the teletype reformat- 
ting program to generate drive tapes. 

Pointing after one real-time sample pull The new 
program generates commands to point the antenna after 
one real-time tape sample pull instead of two as was 
required in the original program. This method reduces 
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Except for the countdown clock, SMC Phase I is 
essentially complete. Problems have been encountered 
in the procurement of components to be used in Hi-Rel 
circuit modules from which the clock is to be made and 
are causing this portion of SMC Phase I to be approxi- 
mately six months behind schedule. 

The prototype program alarm and control panel has 
been delivered to DSS 11, where it is being evaluated. 
Follow-on procurement for five additional panels should 
be complete by March 1967. The graphical recorder 
should be delivered by December 1966 and evaluation 
completed in early 1967. Follow-on procurement should 
be initiated early in fiscal year 1968. 

4. Telemetry and Command Processors, Phase II 


which will show details of design in the digital divider 
and clock areas, as well as an expansion of the functional 
blocks in the simulation and frequency standard sub- 
assemblies. The functional specification should be ready 
for release by the end of the year. 

Procurement of components for the digital modules 
is almost complete, with approximately 30% of the 
required parts now on hand. The Varian R20 rubidium 
frequency standard has been delivered and is currently 
under test. The most recent results indicate that per- 
formance is within specifications. However, there are 
still some problems involving the microphonics and the 
physical characteristics of the equipment, and these are 
undergoing thorough investigation at the standards lab- 
oratory at Goldstone. 


The telemetry and command processors, Phase II 
(TCP II) are the major assemblies in the telemetry and 
command data handling subsystem, Phase II (TCD II). 
A description of the TCP II has been presented in 
SPS 37-38 , Vol. Ill, pp. 76-77. The major changes in 
TCP II during this period are as follows: 

(1) The TCP II-B configuration now operational at 
DSS 51 will be changed to the standard TCP II. 
Procurement has been initiated for the required 
computers, and it is estimated that the changeover 
will take place in mid-1967. 

(2) Procurement of analog to digital converters 
(A/DCs) for installation at DSS 11 and DSS 12 
has been initiated. The addition of these units to 
the TCP II will provide a needed capability for 
additional flight project support. A/DCs will be 
added to the remaining Deep Space Stations at 
a later date. 

(3) A fourth teletype channel has been scheduled for 
addition to the TCP II communications buffer. The 
addition of this channel, scheduled for March 1967, 
will provide the capability for data transfer between 
the TCP II and the DIS. 

(4) TCP II implementation has been proceeding on 
schedule and is in agreement with Table 2, p. 126, 
SPS 37-41 , Vol. III. The only change to this table 
is the configuration change (from TCP II-B to the 
standard TCP II) at DSS 51. 


FTS II progress is within scheduled guidelines, as 
indicated at the design review; the specific problem areas 
discussed at that time are in the process of being resolved. 
Of particular interest at that time were the configuration 
of the timing portion of the simulation equipment and 
the definition of voltage and current levels at the FTS II 
simulation interface. 


Prototype fabrication is scheduled for the first quarter 
of calendar year 1967. 




lets of the Operational Voice 
Communications Subsystems, F. E. Bond, Jr., and 
J. W. Capps 


This paper provides an overall description of the mam 
agement, configuration, and control of the operational 
voice communications subsystem (OVCS), which is a 
communications capability of the SFOF, fulfilling the 
requirement of end terminations and internal distribution 
of audio intelligence. Since this is a leased facility, JPL 
and associated contractor personnel may not relocate, 
move, disconnect,nor alter in any way, any of the equip- 
ment associated with this cuiimiunications oubsystem * 


The OVCS provides (1) touch-tone signaling telephones, 
(2) voice conference nets (for talk-listen headsets/handsets 
and speakers), and (3) point-to-point intercommunications. 


5. Frequency and Timing Subsystem, Phase II 

The design of frequency and timing subsystem, Phase II 
(FTS II) is now in its final phase and is estimated as 80% 
complete. An expanded block diagram is being prepared 


Any person directly involved in space flight operations 
in the SFOF has three voice communications interfaces: 


2 Configuration changes may be made only under the procedures 
outlined in SFOF/SOP 10-201A. 
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the initialization time required in the drive tape mode by 
one sample interval time. 

Station visibility check in the star-track mode. The new 
program performs a check in the star-track mode to de- 
termine if the star to be tracked is visible to the local 
station. If the star is not visible, the program zeroes all 
command registers and waits for the subsystem operator 
to enter a new set of coordinates. The criteria to deter- 
mine station visibility are as follows: 

(1) Local hour angle is between ±90 deg. 

(2) Declination is between ±90 deg. 

(3) Local elevation angle is > 0 deg. 

b. Ascension Island (DSS 72) antenna pointing pro- 
gram. A new antenna pointing program (DOI-5244-OP) 
has been written for DSS 72. The DSS 72 program is 
essentially identical to the antenna pointing program 
(DOI-5107-OP-A) that will be used at DSS 11, 12, 41, 
42, 51, 61, and 62. The major difference between the two 
programs is in the star-track mode. The Ascension Island 
program generates commands to point the antenna in the 
azimuth/elevation coordinate system, while the program 
for the other DSIF stations generates commands in the 
polar coordinate system. This change is necessary to point 
the 30-ft azimuth/elevation mounted antenna at Ascen- 
sion Island. 

The final configuration antenna pointing subsystem 
(APS I), described in SPS 37-38, Vol. Ill, will be installed 
at DSS 14, beginning November 1, 1966. Hardware and 
software checkout will continue for approximately two 
months. The APS I is scheduled to be operational at 
DSS 14 by January 1, 1967. 

2. Digital Instrumentation Subsystem, Phase II 

The digital instrumentation subsystem, Phase II (DIS 
II) prototype design is currently being evaluated, prior 
to processing the procurement for subsequent units. The 
DIS II configuration is an extension of the DIS I, and 
provides the capability for performing the DSIF system 
monitoring functions under control of the system monitor 
program. The DIS II is composed of the following major 
functional elements: (1) SDS 920 digital computer, in- 
cluding peripheral equipment, (2) computer ancillary 
input/output devices, (3) magnetic tape recording ap- 
paratus, (4) data gathering equipment, (5) communica- 
tions buffer, and (6) alarm, control, and display elements. 


The primary objectives of the subsystem are grouped 
into three general categories: (1) performance- and 
status-data presentation, (2) alarm monitoring for failure 
detection, and (3) preparation of permanent historical 
record of station performance. 

The present effort involves definition of the subsystem 
functions, determination of the hardware configuration 
to accommodate these functions, establishing the monitor 
system program criteria, definition of the subsystem 
interface characteristics, and preparation of supporting 
documentation requirements prior to initiation of a pro- 
curement. These activities entail the selection of critical 
station parameters for alarm monitoring and display. 
Formats for periodic reports and alarm messages are 
being established to present information in an optimum 
manner for interpretation by station operating personnel. 
The principal functional elements of the system monitor 
program are being established to permit work to proceed 
on the actual program implementation. Formats for tele- 
type and high-speed data messages transmitted or re- 
ceived external to the station are being coordinated with 
cognizant personnel to assure compatibility with the 
ground communication system requirements, the terminal 
equipment, and the ultimate recipient of the messages. 
Similarly, compatibility in the form and format of station 
performance data recorded on magnetic tape is to be 
established to permit post-mission processing for analysis 
at the SFOF. The present monitoring capability of the 
DIS I, under control of the interim monitor program, 
has no provision for communication of station perform- 
ance and alarm messages to the SFOF. The definition 
of these functions is now being undertaken, along with 
a major expansion in the system functions to be moni- 
tored and displayed for DSIF station operating personnel. 

The procurement for subsystem supporting documen- 
tation is being processed simultaneously with the sub- 
system hardware procurement to alleviate installation 
and testing difficulties and to provide complete operation 
and maintenance procedures at the time of installation. 

3. Station Control and Monitor Console 

The station control and monitor console (SMC) sub- 
system provides a centralized location for the station 
manager to monitor critical station functions and to con- 
trol station operations by directive. The SMC has five 
17 Vi- by 19-in. display-panel surfaces. The SMC displays 
mission-independent data for the most part; however, 
it provides for certain mission-dependent display areas. 
No control functions are available at the SMC. 
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(1) To the higher level of authority (first order). 

(2) To lateral levels of authority (second order). 

(3) To subordinate levels of authority (third order). 

The first level represents the requirement for the 
space flight operations director (SFOD), and assistants, to 
properly direct the supporting elements of his mission. 
To fulfill this requirement, the “command net” of any 
given project must be provided to each supporting 
agency — e.g., Comm, Flight, or Track. It should be noted 
that this prime net is furnished by the project supported. 
Business transacted on this net should be restricted to 
that necessary to control, direct, and determine mission 
status for the particular project involved. 

The second level represents the requirement for the 
directors of the major supporting agencies, themselves, 
to transact business. Each supporting agency has been 
provided with a “home net” which is accessible to the 
remaining agencies. The concept in this case is that each 
agency will monitor and receive calls from any other 
lateral agency only on his home net, a mandatory require- 


ment that eliminates any one person from being respon- 
sible for monitoring (and responding) on a multitude of 
nets. Business transacted on these second-order nets 
should be restricted to matters of immediate concern to 
the two parties conversing. 

The third level represents the requirement for a super- 
visor at any specific location to control those personnel 
directly under his supervision. For the most part, these 
“internal nets,” should be restricted to the pertinent area, 
room, or activity concerned. 

Net custodians have been assigned to all existing 
OVCS voice conference nets. Custodial control is defined 
as the authority to determine fixed net configurations. 
All requested changes to a specific voice net will be 
approved by the assigned net custodian prior to imple- 
mentation. Assigned custodians are listed in Table 2, 
which also reflects the presently existing OVCS voice 
conference nets. 

Documentation of the current configuration of the 
OVCS is provided, and updated for all net custodians. 


Table 2. OVCS voice conference nets 


No. 

Name 

Category 

Type 

Custodian 

3 

ETR 

Ml 

OCC 

i 

D. G. Tustin 


Thi* net is normally interconnected to a scheduled voice line to AFETR and 
is used to coordinate flight data between the FPAC at JPL and the RTCF at 
AFETR. It may also be used for discussion between the SFOF and JPL opera- 
tions at AFETR. When not connected to AFETR, this net may be used as a 
backup net for facility or project use as scheduled. 


4 


STATUS 


Ml 


OCC 


D. G. Tustin 


This net is normally interconnected to a scheduled voice line to AFETR and 
is used to report spacecraft range, countdown, and launch status to the SFOF 
from AFETR and to report DSN and SFO status to the AFETR from the SFOF. 
When not connected to AFETR, this net may be used as an internal distribution 
net, under the cognizance of the operations control chief (OCC), either for backup 
purposes or for wide drstriburion of uudio date. 


1 5 DSN 


Ml 


OCC 


D. G. Tustin 


This is the home net where the DSN-OCC (OPSCON) may be reached when 
on duty. It will be used for overall coordination of the DSN by the OCC and 
the operating system chiefs, and may be used by the latter for interface 
coordination. 


20 


COMM 1 


Ml 


COMM 


C. P. Le Veau 


This net is used for coordination purposes to establish voice contact to 
various elements of the DSN communications center from both DSN and project 
personnel. 


No. 

Name 

Category 



16 

COMM 2 

Ml 

COMM 

C. P. Le Veau 

This net is used as the prime internal net of the DSN communications center 
to permit coordination of operational internal matters within this agency. 

47 

COMM 3 

Ml 

COMM 

C. P. Le Veau 

This net is used as a maintenance/engineering net to which external activities 
are connected during those periods when communication-to-communications con- 
tact is required. 

36 

COMM 

CONF 

Ml 

COMM 

C. P. Le Veau 

The purpose of this net is threefold: (1) to provide a facility whereby groups 
may effect a conference with the DSS stations or other outside agencies without 
interfering with project-oriented nets; (2) to provide a facility whereby dictation 
may be received by the TRArCC.V. teletype operator for teletype transmission 
during periods of voice circuit failure to the DSS stations in a teletype confer- 
ence mode; and (3) to provide a listen-only facility to be bridged to other 
operational nets for the use of visiting personnel. 

30 

SIM 1 

Ml 

SUPPORT 

D. A. Nelson 

This net is used to provide contact between various elements of the SDCC 
and the remainder of the SFOF. 

51 

SIM 2 

Ml 

SUPPORT 

D. A. Nelson 

This net is used as an internal net to provide contact between the various 
elements of the SDCC. 
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Table 2 (contd) 


No. 

Name 

Category 

Type 

Custodian 

2 

DATA 1 

Ml 

DATA 

G. R. Emery 


The purpose of this net is to provide a communications interface between 
the project data controller (DaCon) and the operational components of the 
SFOF data system. For consistency, a given project will usually be supported 
on the same data net; however, this arrangement will be reassigned in real 
time by the data chief as required. 

This net shall not be used at any time for communications between users 
outside the SFOF data system without the concurrence of the data chief. 


9 


DATA 2 


Ml 


DATA 


G. R. Emery 


This net provides a communications interface between DaCon and the oper- 
ational components of the SFOF data system. For consistency, a given project 
will usually be supported on the same data net; however, this arrangement will 
be reassigned in real time by the data chief as required. 

This net shall not be used at any time for communications between users 
outside the SFOF data system without the concurrence of the data chief. 


55 


DATA 3 


Ml 


DATA 


G. R. Emery 


This net provides communications between DaCon and certain operational 
components of the SFOF data system. This net will be used as assigned by the 
data chief to supplement coverage provided by the Data 1 and Data 2 nets. 
Principal usage should be restricted to communications between DaCon and 
TPS or, in the case of tests conducted entirely within the data system, between 
DaCon and any operational component which would normally be contacted on 
the Data 1 or Data 2 net. 

With the exception of communications interface between DaCon and the data 
distribution system (DaFlo), this net will not be used at any time for comuni- 
cations between uses outside the SFOF data system without the concurrence of 
the data chief. 


56 


DATA 4 


Ml 


DATA 


G. R. Emery 


The purpose of this net is to provide communications between DaCon and 
certain operational components of the SFOF data system. This net will be used 
as assigned by the data chief to supplement coverage provided by the Data 1 
and Data 2 nets. Principal usage should be restricted to communications between 
DaCon and TPS or, in the case of tests conducted entirely within the data 
system, between DaCon and any operational component which would normally 
be contacted on the Data 1 or Data 2 net. 

With the exception of communications interface between DaCon and DaFlo, 
this net will not be used at any time for communications between users outside 
the SFOF data system without the concurrence of the data chief. 


5 DPA 1 


Ml 


DATA 


G. R. Emery 


This is the home net of the data chief, and shall be used by mission- 
independent personnel to initiate and carry on communications with the data 
chief. 

This net will also be used for communications between the various operations 
portions of the SFOF data system at the discretion of the data chief. 


54 


DPA 2 


Ml 


DATA 


G. R. Emery 


This net will be used for communications within the SFOF data system only. 
This net shall not be bridged, or otherwise extended, under any circumstances, 
into areas of the SFOF that are not in the complete control of the data system. 

Usage of the DPA 2 net within the ground rules indicated above will be at 
the discretion of the data chief. 


No. 

Name 

Category 

Type 

Custodian 

18 

DSS 1 1 

Ml 

DSIF 

C. A. Holritz 

This net is normally interconnected to a scheduled voice line to DSS 11 and 
is used to coordinate flight data between the SFOF and DSS 11. 

19 

DSS 12-1 

Ml 

DSIF 

C. A. Holritz 

This net 
is used to c 

s normally interconnected to a scheduled voice line to DSS 12 and 
oordinate flight data between the SFOF and DSS 12. 

25 

DSS 13 

Ml 

DSIF 

C. A. Holritz 

This net is normally interconnected to a scheduled voice line to DSS 13 and 
is used to coordinate flight data between the SFOF and DSS 13. 

24 

DSS 14 

Ml 

DSIF 

C. A. Holritz 

This net is normally interconnected to a scheduled voice line to DSS 14 and 
is used to coordinate flight data between the SFOF and DSS 14. 

10 


Ml 

DSIF 



C. A. Holritz 

This net is normally interconnected to a scheduled voice line to DSS 41 and 
is used to coordinate flight data between the SFOF and DSS 41. 

11 

DSS 42 

Ml 

DSIF 

C. A. Holritz 

This net is normally interconnected to a scheduled voice line to DSS 42 and 
is used to coordinate flight data between the SFOF and DSS 42. 

13 

DSS 51 

Ml ' 

DSIF 

C. A. Holritz 

This net is normally interconnected to a scheduled voice line to DSS 51 and 
is used to coordinate flight data between the SFOF and DSS 51. 

33 

DSS 61 

Ml 

DSIF 

C. A. Holritz 

This net is normally interconnected to a scheduled voice line to DSS 61 and 
is used to coordinate flight data between the SFOF and DSS 61. 

49 

DSS 62 

Ml 

DSIF 

C. A. Holritz 

This net is normally interconnected to a scheduled voice line to DSS 62 and 
is used to coordinate flight data between the SFOF and DSS 62. 

38 

DSS 71 

Ml 

DSIF 

C. A. Holritz 

This net is normally interconnected to a scheduled voice line to DSS 71 and 
is used to coordinate flight data between the SFOF and DSS 71. 

50 

DSS 72 

Ml 

DSIF 

C. A. Holritz 

This net 
is used to 

is normally interconnected to a scheduled voice line to DSS 72 and 
coordinate flight data between the SFOF and DSS 72. 

53 

SDA 

Ml 

DSIF 

D. Chaney 

This net 
Goldstone 
scientific d 
minations. 

is normally interconnected to a scheduled circuit to Building 33 at 
and is used for voice liaison between Goldstone, FPAC 1, and 
ata analysts (SDA) in the net control area relative to orbit deter- 
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Table 2 (contd) 


No. 


Name 


Category 


Type 


Custodian 


34 


NETCON 1 


Ml 


DSIF 


C. A. Holritz 


This net is used by DSIF control-room personnel to coordinate their activities 
on a particular mission. Extensions of this net to other areas are made when 
DSIF personnel are located in areas outside the DSIF control room. 

DSIF operations chief assigns use of the net and approves extensions. 





NETCON 2 

Ml 

DSIF 

C. A. Holritz 

Same as Netcon 1 . 

17 

NETCON 3 

Ml 

DSIF 

C. A. Holritz 

Same as Netcon 1. 

1 

PIONEER 

COMMAND 

MD 

PIONEER 

J. W. Davis 

— 


This is the major Pioneer net, which is used for the control of mission oper- 
ations within the SFOF, interconnecting the chiefs of all the prime operational 
and support areas concerned with the Pioneer mission. 

This net will be used to report the sequence of events as they are accom- 
plished during the launch operations. This will provide the SFOD with a 
continuing picture of the progress of the sequence of events in support of the 
Pioneer flight mission before launch, during the boost phase, and during the 
critical phases of the flight mission. This net will also provide the SFOD with 
operational readiness reports, establishing areas of responsibility and readiness 
to support the Pioneer launch operations. 


21 


AMES 


MD 


PIONEER 


J. W. Davis 


This net is normally interconnected to ARC and is used primarily to provide 
mission control for the Pioneer flight mission and, secondarily, to provide 
status reporting for NASA/ARC during the launch operations and the early 
phases of the Pioneer flight mission. 


22 


TRW 


MD 


PIONEER 


J. W. Davis 


This net is normally connected to TRW Systems, providing voice communica- 
tions between the SFOD console at JRL/SFOF and TRW Systems, Redondo Beach, 
California. This voice net will be made available to the SFOD to discuss 
areas concerned with the Pioneer spacecraft contractor. 


28 


PION 

COORD 


MD 


PIONEER 


J. W. Davis 


The purpose of this voice net is to provide the spacecraft performance 
analysis team with an internal JPL/SFOF communications system for technical 
evaluation and reporting of rne Pioneer Vtl spacecraft performance Nominal 
and anomalous detailed discussions will be held on this net to provide the 
SPAC director with the proper information relating to his responsibility. 


52 


PION S/C 


MD 


PIONEER 


J. W. Davis 


The purpose of this voice circuit is to provide the SFO team with the proper 
coordination between the flight operations control at JPL/SFOF and the launch 
operations control at Cape Kennedy Air Force Station (CKAFS). Discussions 
between TRE 1 at JPL/SFOF and the spacecraft representative at Hangar AM, 
CKAFS, will be concerned with detailed coordination of launch operations and 
spacecraft analysis being performed as part of the fight mission operations. 
This net will also provide the spacecraft analysis team at JPL/SFOF and, in 
turn, the JPL/DSN with spacecraft frequency report measurements being taken 
at Hangar AM. 


No. 

Name 

Category 

Type 

Custodian 

39 

FLITE 

CONT 

MD 

LUNAR 

ORBITER 

R. B. Rung 

This mission control net is utilized by SFOD, area directors, track chief, etc., 
and carries S/C status, SPAC and FPAC milestone reports, computer configura- 
tion, command information, and approval. 

40 

FPAA 3 

MD 

LUNAR 

ORBITER 

R. B. Rung 

The FPAC set carries all FPAC business, both inter- and intra-area, using 
the Boeing Company and JPL personnel. Orbit determination, trajectory, and 
maneuver information are the prime data carried. 

41 

DSS 1 

MD 

LUNAR 

ORBITER 

R. B. Rung 

This net is used to patch tracking station voice nets to all areas. SNOMAN 
provides patches allowing tracking stations in operation to be heard. DSS 
41, 61, and 12 may appear. The communications status display board or 
SNOMAN correlates a given DSS to a specific net. These nets are used by 
video engineers for photo data evaluation and photo acquisition information. 



MD 


R. B. Rung 

This net is used to patch tracking station voice nets to all areas. SNOMAN 
provides patches allowing tracking stations in operation to be heard. DSS 41, 
61, and 21 may appear. The communications status display board or SNOMAN 
correlates a given DSS to a specific net. These nets are used by video engi- 
neers for photo data evaluation and photo acquisition information. 

43 

SPAA 3 

MD 

LUNAR 

ORBITER 

R. B. Rung 

This is the SPAC director’s net, interfacing with other operational areas 
and carrying only milestone data of SPAC functions, and is used for command 
conferences, major troubles, and analysts’ reports. 

44 

SPAA 4 

MD 

LUNAR 

ORBITER 

R. B. Rung 

This is the SPAC analysts’ net, permitting discussion of telemetry data, 
data formats, minor troubles, command preparation. 

45 

DATA 

CONTROL 

MD 

LUNAR 

ORBITER 

R. B. Rung 

This com 
operation, 
and proces 
modes. 

puter control net carries computer status, telemetry edit tables in 
data being processed, identifies user programs which have been run, 
ses troubles. It also carries information about computer operational 

46 

MAA 

MD 

LUNAR 

ORBITER 

R. B. Rung 

This is the mission advisors net, an internal conference net for contacting 
area directors in emergencies. 

48 

MAA-INT 

MD 

LUNAR 

ORBITER 

R. B. Rung 

This is an internal conference net used primarily for coordination purposes 
for mission advisory personnel. 
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Table 2 (contd) 


% 


No. 

Name 

Category 

Type 

Custodian 

26 

TPAC 

MD 

SURVEYOR 

R. T. Hayes 

This is the SSAC internal net for Telpac group coordination. 

27 

SUR FPAC 

MD 

SURVEYOR 

R. T. Hayes 

This is the FPAC internal net. 

29 

SURV 

CONTROL 

MD 

SURVEYOR 

R. T. Hayes 

This is the primary SFOF internal mission control net, the interface between 
ACE 1 and area directors. 

31 

— 

TV COORD 

MD 

SURVEYOR 

R. T. Hayes 

This is the TVGDHS internal net. 

32 

SURV 

COORD 

MD 

SURVEYOR 

R. T. Hayes 

This is the backup to the Surveyor control net; it also takes care of unplanned 
interfaces and discussions too lengthy for the Surveyor control net. 

37 

SURV SPAC 

MD 

SURVEYOR 

R. T. Hayes 

This is the SPAC internal net. 


No. 

Name 

Category 

Type 

Custodian 

6 

SPAC/AO 

MD 

SURVEYOR 

R. T. Hayes 

This is the SPAC internal net, connected to Building AO voice line from 
prelaunch to about launch plus 1 hr. 

7 

TSAC 

MD 

SURVEYOR 

R. T. Hayes 

This is the SSAC internal net for TSAC group coordination. 

8 

TV DATA 

MD 

SURVEYOR 

R. T. Hayes 

This TV-1 /TV- 1 1 /Telpac technical coordination net is used for (Goldstone) pre- 
visibility countdown and system calibration; it is disconnected on countdown 
completion. 

12 

FPAC/RTCC 

MD 

SURVEYOR 

R. T. Hayes 

This is the FPAC internal net, connected to the RTCC voice line from prelaunch 
to about launch plus 3 hr. 

14 

SUR 

COMMAND 

1 

MD 

SURVEYOR 

R. T. Hayes 

This SFOF/DSS interface for command activity is bridged with the DSS net 
and connected to station voice line during station transmitting period. 
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